UNCLASSIFIED 


_ AD  NUMBER _ 

AD460593 

LIMITATION  CHANGES 
TO: 

Approved  for  public  release;  distribution  is 
unlimited. 


FROM: 

Distribution  authorized  to  U.S.  Gov't,  agencies 
and  their  contractors ; 

Administrative /Operational  Use;  JAN  1965.  Other 
requests  shall  be  referred  to  Office  of  Naval 
Research^  Washington,  DC  20360. 


_ AUTHORITY 

onr  Itr,  28  jul  1977 


THIS  PAGE  IS  UNCLASSIFIED 


THIS  REPORT  HAS  BEEN  DELIMITED 
AND  CLEARED  FOR  PUBLIC  RELEASE 
UNDER  DOP  DIRECTIVE  5200,20  AND 
NO  RESTRICTIONS  ARE  IMPOSED  UPON 
ITS  USE  AND  DISCLOSURE, 

DISTRIBUTION  STATEMENT  A 

APPROVED  FOR  PUBLIC  RELEASE; 
DISTRIBUTION  UNLIMITED, 


UNCLASSIFIED 


AD- 4  6  0  5  9  3 

/ 

I 

DEFENSE  DOCUMENTATION  CENTER 


SCltNlKIC  AND  TECHNICAL  INFORMATION 

CiWi’CH  suites  ilEUSORIA  VIRGINIA 


UNCLASSIFIED 


NOTICT:  When  government  or  other  drawings,  specl- 
fUcaClons  or  other  data  are  used  for  any  purpose 
otclier  than  in  connection  with  a  definitely  related 
gCQTenunent  procurement  operation,  the  U.  S. 
Gcaremment  thereby  incurs  no  responsibility,  nor  any 
otalLgatlon  ;Aat  soever;  and  the  fact  that  the  Govern - 
meaut  may  have  formilated,  furnished,  or  in  any  way 
suiyplled  the  said  drawings,  specifications,  or  other 
dm1;a  la  not  to  be  regarded  by  implication  or  other- 
wUbc  as  in  any  manner  licensing  the  holder  or  any 
oCher  person  or  corporation,  or  conveying  any  rl^ts 
orr  permission  to  manufactxire,  use  or  sell  any 
pwtented  invention  that  may  in  any  way  be  related 
tfae  ret  o . 


CATALOGED  BY;  DDC 


^  The  Synthesis  of  Sweep-Frequency 
^  Ground  Backscatter  by  Digital  Computer 

by 

T.  A.  Croft 

L  460593 

January  1965 

Technical  Report  No.  84 

Prepared  under 

Office  of  Naval  Research  Contract 
Nonr-225(64),  NR  088  019,  and 

Advanced  Research  Projects  Agency  ARPA  Order  196-65 

RnOIOSCIEnCE  lrborrtorv 

STnnFORD  ELEITROnilS  IRBORnTORIES 

STnnFORO  UniUERSITV  •  STSHFORD,  CnilFORIlin 


DDC  AVAILABILITY  NOTICE 


Qualified  requesters  may  obtain 
copies  of  this  report  from  DDC. 
Foreign  announcement  and  dissem 
ination  of  this  report  by  DDC  i 
not  authorized. 


THE  SYNTHESIS  OF  SWEEP-FREQUENCY 
GROUND  BACKSCATTER  BY  DIGITAL  COMPUTER 


SEL-65-002 


by 

T.  A.  Croft 


January  1965 


Hep  reduction 
IS  permitted 
the  Tn I t  ed 


I  li  whole  or  in  part 
fo  r  any  pu  rpo  s  e  of 
States  Government. 


Technical  Report  No.  84 
Prepared  under 

Office  of  Naval  Research  Contract 
Nonr-225(64) ,  NR  088  019,  and 
Advanced  Research  Projects  Agency  ARPA  Order  196-65 


Radioscience  Laboratory 
Stanford  Electronics  Laboratories 
Stanford  University  Stanford,  California 


f 


f 


ABSTRACT 


A  high-frequency  pulse  radar  often  receives  energy  which  has  propa¬ 
gated  through  the  ionosphere  to  the  ground,  scattered,  and  then  returned 
to  the  radar  again  via  the  ionosphere.  When  this  energy  is  displayed  as 
a  function  of  both  time  delay  and  operating  frequency,  the  record  is 
called  sweep -f requency  backscatter.  In  this  report  a  digital  computer 
technique  is  devised  to  synthesize  such  records  in  a  manner  which  is 
quantitative  in  all  parameters. 

The  synthesis  is  designed  so  that  very  few  approximations  are  needed. 
Full  allowance  is  made  for  the  following;  a  spherical  earth  and  iono¬ 
sphere;  a  wide  variety  of  electron-density  variations  in  height  and  range; 
ionospheric  absorption  which  varies  with  frequency,  range,  and  the  ray 
angle;  radar  antenna  gain  as  a  function  of  frequency,  azimuth,  and  ele¬ 
vation;  ground-scattering  characteristics  which  vary  with  frequency,  range, 
and  the  angle  of  incidence;  ground-reflection  coefficients  whicli  vary  with 
frequency  and  range;  the  transmitted  pulse  shape,  duration,  and  power 
level;  receiver  bandwidth  and  total  received  noise  power.  Only  three 
approximations  are  made;  the  discrete  number  system  is  substituted  for 
continuous  numbers;  ray  theory  is  used  with  the  assumption  that  energy 
is  distributed  uniformly  between  adjacent  closely  spaced  i-aypatlis;  the 
efioct  of  the  geomagnetic  field  on  raypath  range  and  group  time  delay  is 
ncg lec  tod . 

The  digital  syntlicsis  lia.s  two  major  uses.  First,  it  is  possible  to 
vary  a  single  parameter  in  a  known  way  and  then  see  the  effect  tliat  i  l)o 
variation  sliould  haxe  on  ground  backscatter.  For  example,  a  rijjple  of 
known  cha  lac  t  or  i  s  t  ics  can  be  insoi-ted  in  the  ionosphere  at  a  known  range 
and  the  resulting  perturbation  in  the  grounci  Ijackscattor  can  Ijo  calculated. 

The  second  major  use  of  the  program  is  the  syntlicsis  of  a  close 
likeness  of  a  specific  experimental  record.  This  synthetic  record  serves 
to  confirm  knowledge  of  tlio  mechanisms  which  play  a  part  in  ground  back- 
scatter.  For  example,  tlio  ground-scat  t  oi-  coefficient  varies  with  fre¬ 
quency,  teri'ain  typo,  and  incidence  angle  in  a  manner'  which  is  not 
presently  well  understood.  The  process  of  making  a  synthetic  record  to 
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match  an  experimental 
cal  scatter  functions 
behavior.  Since  most 
predictable  behavior, 
pattern  and  study  it 


record  will  Involve  the  trial  of  several  hypotheti- 
until  one  is  found  which  duplicates  observed 
of  the  other  contributing  mechanisms  exhibit  a 
it  will  thus  be  possible  to  isolate  the  scatter 
ndependently , 
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I.  INTRODUCTION 


In  the  study  and  use  of  radio  frequencies  in  the  high-frequency 
band,  a  commonly  encountered  phenomenon  is  the  reception  of  energy  which 
is  scattered  off  the  ground  and  propagated  to  the  receiver  after  iono¬ 
spheric  refraction.  This  "ground  backscatter"  arises  in  an  easily 
recognized  form  when  the  source  of  energy  is  a  pulse  transmitter  located 
near  the  receiver.  Echoes  of  ragged  and  fading  structure  are  often 
found  after  a  time  delay  of  a  few  milliseconds. 

When  backscatter  was  first  discovered,  it  was  attributed  to  the 
action  of  reflectors  of  some  soi-t  located  in  the  ionosphere  [  Ref .  l]. 

Only  in  recent  years  was  it  firmly  established  that  the  predominant 
source  of  such  echoes  is  scatter  from  the  ground  itself,  after  propaga¬ 
tion  through  the  ionosphere  [Refs.  2  and  s].  Several  historical  accounts 
of  l)ackscattcr  studios  are  available  [Refs.  4,  5,  and  6],  so  another  will 
not  Ijo  given  here. 

A.  DEFINITION  OF  THE  PROBLEM 

The  term  "backscatter"  will  be  used  lici-e  in  a  restricted  sense  to 
moan  i  onosijher  ical  ly  propagated  energy  scattered  by  tlie  ground.  Tlie 
phenomenon  wliich  will  be  anal\zo(i  is  tlie  eclio  whicli  would  he  received  t)y 
a  pulse  radar  in  the  high-frequency  banci  wlien  the  solo  source  of  return¬ 
ing  energy  is  ground  scatter.  This  source  accounts  for  most  of  the 
echoes  which  are  received  by  lif  radai’,  anci  oilier  echoing  sources  have 
been  analyzed  in  considerable  fletail  elsewhere.  Examples  are  the  fol¬ 
lowing:  f i old -al igned  ionization  in  the  vicinity  of  the  auroral  zone 

[itef.  7],  in  the  temper  t  ure  latitudes  [liefs.  8  and  o],  and  near  the 
equator  [Ref.  lO];  direct  echoes  from  sporadic  E  patches  [itef.  1 1  ] ; 
apparent  echo  surfaces  duo  to  local  dense  regions  or  turhulence  in  the 
E  and  F  i-egions  TRof .  I2];  man-made  regions  of  dense  ionization  [Ref.  is]; 
incoherent  electron  scatter  [uof .  14].  The  analysis  of  these  phenomena 

is  in  most  cases  less  complex  than  it  is  for  the  case  of  ground  back- 
scatter  because  they  invtilve  the  interaction  of  fewer  mechanisms.  A 
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great  deal  has  been  written  about  ground  backscatter,  and  it  is  well 
summarized  in  the  four  lists  of  abstracts  which  are  available  [Refs.  Iv5, 
16,  17 ,  and  18  ]  . 

Below  the  maximum  plasma  frequency  of  the  ionosphere  f  the 

c  ’ 

strongest  source  of  echo  is  the  multiple  return  from  vertical  incidence. 
This  return  is  usually  so  strong  that  it  masks  the  backscatter  echo. 
Consequently  this  analysis  will  be  limited  to  frequencies  above  f 

c 

(usually  3  to  6  Me)  whore  backscatter  can  commonly  bo  seen. 

These  frequencies  are  well  above  the  electron  gyrof requency ,  which 
remains  below  1.5  Me  in  all  the  regions  traversed  by  energy  which  plays 
a  pajt  in  backseat tering  [Ref.  19].  Other  workers  have  investigated  this 
matter  and  concluded  that  the  geomagnetic  field  causes  changes  in  the 
raypath  parameters  which  are  so  small  that  the  field  may  safely  be 
neglected  [Ref.  20 ] .  This  is  an  extremely  fortunate  circumstance  because 
it  leads  to  a  considerable  simplification  in  the  mathematics  with  an 
attendant  decrease  in  the  cost  of  computation.  The  upper  frequency 
limit  of  interest  is  established  by  the  maximum  frequency  that  can  be 
icfractod  back  to  the  ground  by  the  ionosphere  under  the  average  con¬ 
ditions,  roughly  25-35  Me. 

With  the  tinders  landing  of  the  nature  of  ground  backscatter  came  a 
realization  of  its  potential  value  as  a  means  for  study  of  the  refracting 
characteristics  of  tlie  ionosphere.  Using  a  pulse  radar  in  the  higii- 
f requency  band,  one  can  acquire  a  consideral)le  amount  of  information 
aljout  tlie  nature  of  the  ionosphere  at  a  distance  from  the  equipment  in 
a  desired  direction  [kef.  21  ] .  This  technique  show's  its  greatest  promise 
wiien  the  transmitter  and  receiver  are  made  to  sweep  through  the  high- 
frequency  band  so  that  the  backscatter  can  be  studied  as  a  function  of 
botii  time  delay  and  the  frequency  of  the  radio  waves  [Refs.  22  and  23], 
These  "sweep-frequency  ground -backsc a  I  ter "  records  have  been  gathered 
for  a  nutnijcr  of  years;  the  details  of  their  structure  are  complex  and 
interesting  although  not  now  well  understood. 

If  sweep-frequency  ground -backsca t ter  records  could  be  properly 
interpreted  to  shed  light  on  the  nature  of  the  ionosphere  in  ti-e  region 
of  refraction,  tiiey  would  be  of  groat  use  to  the  many  organizations  which 
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use  the  high-frequency  band  as  a  means  of  communication.  These  organi¬ 
zations  form  an  important  and  large  segment  of  the  communications  indus¬ 
try,  primarily  by  virtue  of  the  large  bandwidth  available  and  the  ability 
of  these  waves  to  refract  off  the  ionosphere  to  reach  around  the  curva¬ 
ture  of  the  earth. 

Also,  it  should  be  noted  that  the  proper  interpretation  of  back- 
scatter  records  would  be  valuable  in  the  study  of  the  upper  atmosphere, 
a  region  of  increasing  importance  with  the  advent  of  nuclear  weapons  and 
rocket  travel.  The  ability  of  a  backscatter  sounder  to  examine  the 
ionosphere  at  considerable  lateral  distance  permits  study  of  areas  which 
are  otherwise  inaccessible  due  to  geographical  or  political  limitations. 
Backscatter  techniques  offer  the  possibility  of  examining  the  ionosphere 
at  a  distance  from  a  single  station  which  operates  at  a  fixed  location 
without  remote  assistance  or  use  of  remote  equipment  (Fig.  l).  This 
single  fact  holds  out  so  much  promise  that  considerable  effort  is  justi¬ 
fied  in  an  attempt  to  understand  the  phenomenon.  It  may  prove  feasible 
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FIG.  1.  THE  PRIMARY  ADVANTAGE  OF  BACKSCATTER. 


to  study  the  detailed  behavior  of  the  ionosphere  at  a  distance  over 
areas  that  are  inaccessible;  oven  over  accessible  areas  the  backscatter 
technique  offers  the  advantage  of  technical  simplicity  when  compared  to 
the  various  alternate  methods-  that  use  multiple  cooperating  stations. 
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The  theoretical  study  of  the  phenomenon  has  led  to  a  rather  dis¬ 
appointing  lack  of  success.  In  fact,  the  reader  will  find  that  this 
paper  gives  the  first  analysis  of  such  records  which  is  realistically 
quantitative  regarding  signal  strength.  Since  a  backscatter  record  is 
merely  a  time  history  of  signal  strength,  it  can  be  seen  that  this  is 
a  rather  basic  quantity  which  must  be  understood  if  the  data  are  to  be 
properly  interpreted. 

Most  of  the  information  successfully  derived  from  backscatter  records 
has  been  taken  from  a  study  of  the  frequency  dependence  of  the  leading 
edge  of  the  echo.  This  parameter  can  be  studied  without  analysis  of  the 
signal  strength  because  a  phenomenon  called  "time-delay  focusing"  usually 
causes  the  backscatter  echoes  to  begin  abruptly  so  that  the  skip -di stance 
portion  of  the  return  can  be  identified  with  some  confidence  [kef.  24  . 

One  useful  application  of  this  has  been  the  determination  of  the  MUF  of 
a  propagation  path  by  means  of  the  equipment  located  solely  at  one  end 
of  the  path  [kef.  25]. 

B.  METHOD  OF  STUDY 

Tlie  backscatter  phenomenon  is  a  complex  interaction  of  so  many  dif¬ 
ferent  mechanisms  that  the  most  practical  approach  to  its  study  is  the 
use  of  a  modern  high-speed  digital  computer.  The  present  method  taxes 
tlie  capability  of  an  IB.M  7090--a  fairly  modern  machine,  which  may  be  one 
of  the  primary-  reasons  why  it  has  taken  so  long  to  acquire  an  ability  to 
study  backscattei-  signal  strength.  Analytic  and  graphical  methods  are 
useful  only  if  the  inv'estigator  is  willing  to  make  some  extreme  simplify¬ 
ing  assumptions  (i.e.,  a  parabolic  ionosphere)  or  if  he  is  interested  in 
just  those  aspects  of  the  phenomenon  which  can  be  attacked  by  some  non¬ 
computer  technique  (i.e.,  the  time  delay  of  the  leading  edge  of  the  echo). 
Some  authors  have  attempted  to  calculate  the  signal  strength,  but  their 
methods  have  eitlier  been  too  difficult  for  practical  use  [Ref.  26  i  or 
else  they  have  made  the  calculation  short  by  means  of  unrealistic  simplify¬ 
ing  approximations  [kefs.  20  and  27  , 

Work  has  been  done  in  determining  skip-distance  time  delay  of  simple 
ionospheres  by  analog  computer  [kef.  28],  but  the  method  did  not  give  any 
information  about  signal  amplitude  and  so  its  usefulness  was  limited. 


SEL-65-O02 


4 


The  approach  to  be  described  begins  with  a  rather  detailed  analysis 
of  the  backscatter  phenomenon.  Subsequently,  the  results  of  this  analy¬ 
sis  are  used  to  devise  a  computer  program  which  will  simulate  the  various 
mechanisms  that  comprise  backscatter,  leading  to  a  synthetic  radio  record 
which  might  be  acquired  by  a  high-frequency  pulse  radar.  The  program  is 
run  at  a  number  of  frequencies,  and  finally  a  simple  manual  operation  is 
used  to  combine  many  such  records  into  a  single  "sweep-frequency  ground- 
backscatter  record." 

Both  analysis  and  synthesis  make  use  of  raypaths  that  are  calculated 
through  a  known  ionosphere  by  computei’.  The  synthesis  process  makes  use 
of  "raysets,  "  which  are  groups  of  numbers  stored  on  punched  cards  in  such 
a  manner  as  to  permit  the  subsequent  reconstruction  of  the  characteristics 
of  a  family  of  raypaths  at  low  cost. 

Since  raypaths  may  be  considered  to  be  streamlines  in  a  flux  of 
energy,  the  basic  assumption  is  made  that  energy  which  starts  out  between 
closely  spaced  raypaths  will  stay  between  them,  uniformly  spaced  in  both 
range  and  time,  and  thus  the  computer  can  calculate  the  direction  of 
travel  and  the  distribution  of  energy  in  space  and  time.  This  approxi¬ 
mation  approaches  truth  so  long  as  the  raNS  are  closely  spaced.  A  con¬ 
siderable  cost  reduction  is  achieved  by  neglecting  the  magnetic  field, 
an  approximation  which  is  good  in  this  application.  These  approximations, 
together  witli  the  use  of  a  discrete  lather  than  a  continuous  number  sys¬ 
tem,  ai'c  the  only  ones  needed;  thus  the  "rayset  method"  is  shown  to  Ije 
an  economical  approach  for  the  realistic  calculation  of  the  amplitude  of 
a  ground -backseat  ter  signal.  Each  run  costs  about  two  dollars,  and  the 
result  is  au t oma t ica ] 1 y  plotted  for  an  extra  dollar. 

C.  APPLICATIOXS  OF  THE  TECHNIQUE 

The  backscatter  synthesis  program  can  be  used  as  a  means  of  studying 
the  phenomenon  it.self.  Because  the  synthesis  technique  is  so  lealistic, 
calculations  can  .show  the  effec'  on  giound  backscatter  of  changes  in  the 
following: 

1.  The  shape  and  elevation  angle  of  the  lobes  in  antenna  gain  patterns. 

2.  Secondary  lobes  in  the  gain  patterns. 
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3.  The  clegi'ce  of  ionospheric  absorption. 

4.  Changes  in  the  absorption  as  a  function  of  j-ange,  such  as  when  the 
radar  looks  across  a  twilight  zone. 

5.  The  duration  of  the  transmitted  pulse. 

6.  The  shape  of  the  transmitted  pulse,  including  the  effect  of  limited 
receiver  bandwidth  which  tends  to  smooth  the  pulse. 

7.  The  absolute  magnitude  of  the  scattering  coefficient  of  the  ground. 

8.  The  dependence  of  this  scattering  coefficient  with  ground  range  as, 
for  example,  when  a  radar  looks  out  ovei-  flat  land,  mountains,  and 
sea . 

9.  The  dopondence  of  the  scattering  coefficient  on  the  angle  at  wliich 
enoj'gy  strikes  the  ground. 

10.  Variation  of  scattering  coeffieient  with  frequency. 

11.  Tlie  reflecting  coefficient  of  the  ground,  which  is  a  measui-e  of  the 
amount  of  energy  specularly  reflected.  Tliis  affects  tlie  relative 
power  of  multiple-hop  backscatter  echoes. 

12.  Variation  of  the  reflecting  coefficient  as  a  function  of  ground 
range  duo  to  the  presence  of  the  various  types  of  toi-rain  which 
also  affect  the  scattering  coefficient. 

The  effect  of  changes  in  the  nature  of  the  ionosplioro  can  also  bo  studied 
by  this  technique.  The  kinds  of  changes  that  can  bo  studied  are  limited 
only  by  the  versatility  of  the  raytracing  program  used  to  generate  the 
raysots.  Practically  any  raytracing  pi-ogram  can  be  adapted  to  this  pur¬ 
pose;  t  tie  particular  one  used  to  support  those  calculations  would  permit 
investigations  of  tlio  following  cliangos  in  the  ionosphere; 

13.  Chapman  layers,  parabolic  layers,  E,  ,  and  F.;>  layers  in 
various  combinations  and  relative  proportions;  in  fact,  any  iono- 
splierc  whicli  can  lie  adec|uatcly  described  1)>  3000  sample  points. 
(Usually,  300  serve  t  tie  purpose.) 

14.  tiorizontal  gradients  of  electron  density  (ionospheric  tilts)  which 
are  gentle  but  widespread,  similar  to  those  which  are  normally 
present  in  ttie  ionosptiere  near  dawn  or  twilight. 

15.  An  isolated  irregularity  in  an  otherwise  uiut  i  s  tu  rtied  ionosphere. 

The  i  ri-egula  ri  ty  can  be  some  kinct  of  localizect  ctense  or  rarefiect 
region.  B>-  carrying  out  ttie  synthesis  with  ttie  irregularity  at  a 
succession  of  different  ranges,  it  is  possible  to  simulate  the  time 
tiistoi-y  of  ttie  effect  of  a  traveling  irregularity  on  sweep- f  roc|uoncy 
gi'ound  backscatter. 
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FOUNDATIONS  OF  THE  BACKSCATTER  ANALYSIS 


II  . 


In  preparation  for  the  analysis  of  backscatter,  three  pertinent 
subjects  will  be  discussed: 

1.  The  concept  of  a  rayset  as  a  means  of  storing  raypath  characteristics. 

2.  Experimental  measurements  of  the  ground-backscatter  coefficient. 

3.  Nondeviative  absoiption  formulas,  manipulated  into  a  useful  form. 


A.  THE  METHOD  OF  CALCULATION  BY  MEANS  OF  RAYSETS 

Before  proceeding  with  analysis  of  backscatter,  it  is  necessary  to 
introduce  a  mathematical  concept  which  has  been  called  the  "rayset.” 
Calculations  will  be  keyed  to  the  utilization  of  this  concept. 

It  is  assumed  that  there  exists  a  digital  computer  program  which 
traces  the  paths  of  radio  waves  in  the  ionosphere.  Such  programs  are 
in  operation  at  several  institutions  [Ref,  29]  and  they  form  the  basis 
of  a  rapidly  advancing  technology  in  the  study  of  the  ionosphere.  The 
Itarticular  ray  tracing  program  is  largely  immaterial,  so  long  as  it  has 
the  requisite  versatility,  accuracy,  and  sufficient  speed  to  bo  run  at 
a  reasonable  price.  Tlie  program  used  to  support  this  work  has  Ijoen 
dcscriljod  in  a  previous  publication  [Ref.  30  ] ;  basically,  it  is  a  suc¬ 
cessive  series  of  solutions  of  Snell's  Law  in  a  concent  j-ically  strati¬ 
fied  ionosplu'i-o .  Tlie  index  of  j-efraction  is  a  scalar  which  vai-ios 
voi'tically  and  ho I'i /.on tal ly  as  indicated  on  Fig.  2. 

This  piogram  is  aljlc  to  present  its  results  in  any  comliination  of 
tiii-ee  forms:  the  raypaths  will  be  drawn  in  curved  earth  coordinates; 

selected  numerical  data  can  l>e  printed  out  to  sliow  the  ray  trajectory 
paiainet  ers;  or  tlie  coiiipulor  may  punch  the  data  into  the  rayset  cards. 
Figure  3  shows  the  automatically  drawn  raypaths  which  were  calculated 
through  a  Chapman  ionosphere  whoso  vert  ical  -  incidence  critical  frequenc:, 
was  -1  Me.  These  raypaths  wore  originally  drawn  by  a  small  x-y  plotter 
controlled  by  a  tape  from  the  IBM  7090.  The  true  shapes  of  the  rays  are 
shown  over  a  curved  earth.  To  the  author's  knowledge,  this  is  the  only 
ra\t racing  program  whicii  directly  plots  its  resulting  raypaths  in  an 
undistorted  coordinate  sy.stem. 
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ELECTRON  DENSITY,  N,  IS  A  FUNCTION 
ETC.  OF  RANGE  WITHIN  EACH  SHELL. 

REFRACTIVE  INDEX,  n  =y/|  - 


A  rayset  is  a  set  of  eleven  numbers  which  characterize  a  particular 
hop  of  a  particular  calculated  raypath.  This  set  is  stored  in  one 
punched  card;  the  system  is  so  devised  that  a  deck  of  appropriately 
related  punched  cards  will  contain  sufficient  oi-ganized  information  that 
the  user  can  extract  practically  anything  he  wishes  to  know  concerning 
the  interrelationships  of  the  various  raypaths  that  have  been  computed. 
Thtis,  docks  of  raysets  can  bo  used  as  input  data  in  subsequent  computer 
progj'ams.  In  this  way.  the  information  which  was  acquired  by  raytracing 
can  be  used  over  and  over  at  practically  no  cost. 

New  programs  can  bo  written  based  on  the  assumption  that  the  relation- 
sliips  between  raypath.  parameters  are  known  quantities.  For  example,  the 
dei'iv'at i\o  which  relates  the  change  in  I’aypath  ground  range  that  arises 
from  an  a.ssociatod  change  in  the  raypath  takeoff  angle  can  be  considered 
as  a  known  number;  this  derivative  appears  in  most  expressions  for  energy 
density  wliich  reaches  the  ground.  Consequently,  the  analysis  which  follows 
may  stop  when  it  encounters  such  a  quantity. 

Tliis  is  a  tremendous  advantage  because  it  is  usually  the  unknown 
\ariations  of  such  parameters  which  make  other  forms  of  analysis  either 
difficult  or  impossible.  The  rayset  method  was  previously  used  by  the 
author  for  the  synthesis  of  oblique  ionograms  [lief.  31,.  In  that  case 
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FIG.  3.  EXAMPLES  OF  CALCULATED  RAYPATHS  AUTOMATICALLY  DRAWN. 
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the  analysis  was  very  simple  and  the  subsequent  synthesis  was  made  almost 
trivial  by  the  use  of  raysets. 

The  natui-o  of  a  rayset  can  best  be  understood  by  consideration  of  the 
process  which  goes  on  inside  a  digital  computer  when  ray  tracing  is  taking 
place.  Each  ray  is  calculated  individually;  the  computer  is  instructed 
to  go  to  a  designated  starting  point  and  calculate  the  path  originating 
at  that  point  which  propagates  at  a  given  initial  takeoff  angle  above  the 
surface  of  the  earth  as  indicated  in  Fig.  4.  Using  circular  coordinates 
and  appropriate  equations,  the  computer  calculates  the  curved  path  of  the 
ray  as  it  traverses  a  medium  of  vai-ying  index  of  refraction. 


Manx  things  can  liappen  to  the  ray;  it  can  refract  off  the  ionosphere 

and  cxt'cine  normal  hops.  It  may  go  so  far  that  it  exceeds  a  specified 

maximum  range.  H  on  Fig.  4,  in  which  event  the  computer  will  bo 

max 

instructed  to  shut  off.  Similarly,  there  might  be  a  specified  maximum 

altitude  11  at  which  ravpath  calculations  would  be  terminated, 

'  max 

If  the  index  of  refraction  varies  in  the  horizontal  direction,  the  ray 
may  climb  into  the  ionosphere  only  to  come  down  on  lop  of  a  densely 
ionized  area;  if  the  ray  then  refracts  upward  off  this  area,  it  is  said 
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lo  liavo  oxociitod  a  "porigoe.  "  This  somotimes  happens  when  rays  are 
calculated  in  an  ionosphere  with  an  E  layer  w'hose  electron  density 
increases  with  increasing  range. 

Usually,  a  ray  trapped  between  the  E  and  F  layers  will  either 
escape  or  roach  the  maximum  range.  However,  there  may  bo  a  buildup  in 
the  electron  density  between  the  E  and  F  layers  of  sufficient  mag¬ 
nitude  that  the  index  of  refraction  becomes  imaginary  in  all  foi-ward 
dii'octions.  This  electron-density  distribution  is  called  a  "closed  trap" 
by  the  author,  and  rays  are  terminated  when  they  oncoiintei'  this  circum¬ 
stance.  However,  in  the  actual  situation,  the  enoi’gy  must  reflect  and 
come  back  on  t . 

Each  time  a  ray  encounters  the  oartli,  a  perigee,  oi'  a  raypatli  termina¬ 
tion  instruction,  the  computer  stops  tompoi'arily  to  store  on  magnetic 
tape  a  sot  of  numbers  which  will  later  be  punched  into  a  card  to  t>e 
called  a  raysot.  Those  parameters  are  illustrated  on  Fig.  5.  Consider 
the  first  hop  of  the  ray  shown.  When  the  computer  finds  that  tlie  ray 
strikes  the  earth,  it  stores  the  niimboi-s  shown  in  the  top  row  on  tlio 
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FIG.  5.  P.lRAMETEaS  IN'CLUDEI)  IX  A  RAY.SET. 
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flgiiro.  The  paramotor  IID  (lonosphoro  idon  t  If  Icn  t.l  on)  Is  a  (hroo-dlRll 
numljor  (hnL  idoiitlflos  the  particular  Ionosphere  used  for  calculation. 

The  frequency  F  Is  given  In  kilocycles.  The  parameter  B  (bounce)  Is 
given  the  value  1  because  this  Is  the  first  time  the  ray  lias  hit  the 
earth.  The  parameter  p  Is  the  takeoff  angle,  and  Is  I  he  landing 

angle  as  shown  on  the  figure.  The  ground  range  from  the  transmitter  to 
the  point  of  Impact  of  the  ray,  II^ ,  Is  given  In  kllomelors.  The  group 
and  phase  time  delays,  T^^^  and  T^^,  are  given  In  milliseconds.  Tlie 
height  Is  the  height  of  the  apogee  of  the  ray  as  shown.  The 

paramotor  (indicator)  Is  a  number  which  is  coded  0  through  (>  and 

servos  to  Identify  the  general  nature  of  the  rnypath.  Notice  that  I  -  0 
when  the  ray  hits  the  cartli;  1  =  1  if  the  ray  exceeds  maximum  height 
and  1=2  or  3  if  the  ray  exceeds  maximum  range.  Similarly,  I  =  d 
or  if  the  ray  executes  perigees,  and  I  =  G  if  the  ray  encounters 

a  closed  trap. 

After  the  computer  has  stored  I  lie  raysel  described  above.  It  con¬ 
tinues  calculating  the  I'ay  until  It  again  encounters  the  earth,  a  perigee, 
or  a  raypath  termination  instruction.  In  the  example  shown  on  Fig.  .T , 
the  ray  hit  the  ground  for  a  second  time  and  the  computer  stored  the 
rayset  shown  as  the  second  row  In  the  figui'c.  Here  the  Information  Is 
similar  to  the  first  rayset  except  that  D  =  2,  because  this  Is  the 
second  encounter  with  the  ground.  Also,  the  i-ange,  group  time  delay, 
and  phase  time  delay  are  all  given  as  accumulated  totals  counted  from 
the  t  rnnsmi  t  t.er . 

As  a  mutter  of  Interest,  the  definitions  of  the  various  parameters 
as  shown  on  Fig.  .G  apply  only  when  1  =  0.  When  I  Is  nonzero,  the 
ranges,  times,  and  angles  are  defined  differently  so  as  to  maximize  the 
ability  of  the  user  to  reconstruct  the  ray  t  ra.j  oc  tor  I  es .  For  example, 
if  the  third  hop  of  the  ray  cxccedeil  maximum  range  going  upward,  ns  It 
appears  to  do  In  I'ig.  .G ,  the  time  aiul  range  would  be  measured  to  the 
point  of  ray  termination.  The  height  would  be  the  hi'lght  of  termina¬ 
tion  and  (r  would  lie  the  terminal  angle  of  the  ray  relative  to  hori¬ 
zontal.  The  indicator  would  lie  I  2.  If  the  ray  had  i-xcieeded  maximum 
range  coming  down,  tin;  same  definitions  would  apply  (>xcepl  that  the 
indicator  would  lie  I  3.  Actually,  tin*  change  in  the  indicator  numlier 
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sei'vos  only  to  establish  the  algebraic  sip:n  of  the  angle  This 

procedure  conserves  one  column  on  the  punched  card  and  thus  maximizes 
the  amount  of  information  it  can  convey. 

Figure  6  shows  the  appearance  of  a  family  of  rayset  cards  when  they 
are  printed  out.  The  family  shown  has  been  calculated  through  a  Chapman 
ionosphere  with  a  100-ktn  scale  height  and  a  4-Mc  critical  frequency. 
Thus,  the  operating  frequency  is  twice  the  critical  frequency.  The 
height  of  the  maximum  occurs  at  300  km  and  there  were  no  horizontal 
gradients  of  electron  density  in  this  example.  The  corresponding  rays 
are  drawn  in  the  sixth  plot  on  Fig.  3. 
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Using  the  information  on  those  rayscts,  one  can  reconstruct  most  of 
the  useful  Information  that  describes  all  of  the  possible  one-hop  rays 
through  this  ionosphere.  Raysots  computed  through  a  very  simple  iono¬ 
sphere  such  as  this  will  prove  useful  in  later  synthesis  for  study  of 
the  effect  on  backseat  ter  of  equipment  changes  without  the  added  com¬ 
plication  in  the  record  which  appears  when  the  ionosphere  has  layers  or 

tills. 

n.  THE  GROUND -SCATTER  PATTERN 

In  ID.^O,  Neilson  et  al  [Ref.  32]  and  Hagn  [Ref.  33]  of  the  Stanford 
Research  Institute  (SRT.)  flow  an  in  st  rumen  tod  aircraft  o-'cr  various 
terrain  typos,  investigating  the  behavior  of  backseat  I erod  power  at 
.32.8  Me.  Although  the  tlata  are  restricted  to  one  frequency,  they  appear 
to  be  the  only  comprehensive  data  availatile.  Therefore,  as  a  starting 
point,  the  SRI  data  will  be  used  at  all  frequencies;  however,  it  should 
be  noted  that  the  technique  of  backseat  tor  synthesis  l>ei ng  invented  here 
may  open  an  avenue  for  the  investigation  of  the  actual  behavior  of  the 
ground-b.acUscatter  coefficient  as  a  function  of  frequency,  since  all  the 
other  necessary  parameters  will  be  .seen  to  exhibit  a  predictable  behavior 
during  changes  of  frequency. 

The  SRI  results  were  expressed  in  terms  of  a  pai-amoter,  p.  In  this 
paper,  h  denotes  a  reflection  cootficient,  so  scattering  will  be 
symbolized  K.  On  page  51  of  Ref,  33.  K  is  explained  (in  paraphrase) 
as  lollows: 

"When  K  is  viewed  as  a  loss  parameter,  one  simple  interpretation 
is  to  consider  it  the  ratio  of  the  received  power  actually  mea¬ 
sured  to  the  power  that  would  have  been  measured  with  K  equal 
to  unity.  The  case  of  K  being  unity  corresponds  to  the  earth 
acting  as  an  antenna  (with  aperture  equal  to  the  area  actually 
Illuminated  1)V  the  pulse)  which  receives  the  incident  power  and 
isotropically' retransmits  it  without  depolarization  or  loss  into 
a  hemisphere." 

As  a  consequence  of  the  definition  of  K,  the  theoretical  treatment 
of  backseat  ter  will  begin  with  the  assumption  that  the  eartli  acts  as  a 
completely  reradiating,  hemispherically  isotropic  antenna.  The  power 
received  (according  to  the  resulting  equations)  will  then  bo  multiplied 
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by  K  in  order  to  bring  the  theoretical  predictions  into  line  with  the 

SRI  measurements  or  other  known  measures. 

On  Fig.  7  are  shown  Hagn ' s  curves  depicting  the  behavior  of  the 

parameter  K  as  a  function  of  the  angle  of  incidence  ilf  for  various 

2 

terrain  types.  Also  shown  on  the  figure  are  curves  of  sin  ■>!;,  sin  \lr, 

4 

and  sin  i!'.  The  labels  H  and  V  refer  to  horizontal  and  vertical 

poi:i"i2a‘ ion  of  the  radiation  being  measured.  It  is  somewhat  remarkable 

thai  in  all  cases  the  horizontal  polarization  seems  to  exhibit  a  shape 

4 

which  is  quite  close  to  that  of  the  cui-ve  of  sin  Since  the  vertical 

axis  of  Fig.  7  is  logarithmic,  curves  can  be  shifted  up  and  dow'n  on  the 
plot  simply  by  use  of  a  scalar  multiplioj’. 


FIG.  7.  32-Mc  BACKSCATTER  COEFFICIENTS  MEASURED 

BY  AIRCRAFT. 
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For  vertical  polarization,  the  data  are  not  so  easily  matched;  in 

2 

some  cases  it  appears  that  the  data  will  fit  a  sin  ijr  curve,  whereas  in 
other  cases  it  appears  that  a  constant  would  be  more  appropriate. 

Figure  8  shows  the  same  curves  separated  by  terrain  typo,  together  with 
analytic  functions  which  have  been  selected  as  the  closest  reasonable 
match  to  each.  It  is  seen  that  in  some  cases  the  agreement  is  quite 
close.  The  reader  should  not  bo  excessively  critical  of  small  differ¬ 
ences  between  the  analytic  function  and  the  SRI  measured  curves,  for  the 
latter  are  actually  smoothed  representations  of  a  wide  variety  of  data 
and,  in  a  sense,  they  reflect  the  judgment  of  Hagn. 


ANGLE  OE  INCIDENCE  *  lAegl 


FIG.  8.  ANALYTIC  FUNCTIONS  FITTED  TO  THE  SCATTER  DATA. 
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As  a  consequence  of  this  curve-matching  exercise,  it  was  decided  to 


represent  the  five  types  of  terrain  by  the  function 

4  2 

K(\lf)  =  K  sin  \1;  +  K  sin  +  K 

X  ^  * 

where  the  values  of  ,  K^,  and  are  as  shown  in  Table  1,  which 

gives  the  data  used  to  produce  Fig.  8. 


TABLE  1.  GROUND -REFLECT I ON  COEFFICIENT  FACTORS 


Material 

Horizontal 

Polarizat ion 

Vertical  Polarization 

_ 

'<2 

*^3 

New  ice 

lO'^ 

7.1 

X  10'“^ 

0 

Old  ice 

2.5  X 

10 

2.5 

.  10-^ 

0 

Sea 

3.3  N 

10-^ 

0 

8.3  X  10 

Smooth  moist 

1  and 

1  X 

lo"^ 

0 

2.5  \  lo'^ 

Smooth  dry 
land 

1.8  - 

1 

O 

0 

-6 

2.1  X  10 

Although  the  SRI  data  are  the  onlv  complete  set,  some  other  authors 
have  mentioned  isolated  characteristics  of  ground  scattering  derived 
cither  from  theory  or  experiment.  The  theoretical  study  of  this  suljject 
is  extremely  difficult.  It  appears  that  those  who  try  to  analyze  this 
phenomenon  make  approximations  so  severe  that  their  results  are  question- 
aljle  or  else  they  fire  unable  to  solve  the  equations  that  they  set  up. 
Tlierc  is  a  wide  disagreement  among  the  results  which  have  been  derived 
and  so  they  will  not  l)o  used  hei-e.  However,  the  following  items  taken 
from  the  literature  deserve  mention: 

1.  Sheaiinan  [Ref.  20  ]  states  on  theoretical  grounds  that  the  Ijack- 
scatter  should  vary  as  sin'll]/.  It  is  also  interesting  to  note 
that  ho  made  calculations  which  sliow  tiiat  the  magnetic  field  can 
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be  neglected  in  backscalter  calculations  without  introducing 
perceptible  error  in  the  location  of  the  scatterers. 

2.  McCue  [Ref.  6]  acquired  experimental  data  which  indicated  that 

land  is  a  better  scatterer  than  the  sea  when  the  angle  is 

small . 

3.  Ingalls  and  Stone  [Ref.  34  [  made  measurements  at  18  and  24  Me 
which  seemed  to  indicate  that  the  sea  had  a  backscatter ing  cross 
section  of  10'"^  m^  for  every  square  meter  of  sea. 

4.  Ranzi  [Ref.  35,  found  that  the  sea  was  a  better  scatterer  than 

the  land  by  an  amount  which  varied  between  5  and  16  db.  This  same 
author  [Ref.  36 j  found  that  the  scattering  coefficient  of  a 
mountain  range  depended  very  sti’ongly  on  the  particular  charac¬ 
teristics  of  the  troposphere.  It  appears  that  the  refracting 
power  of  the  troposphere  varies  slightly  from  day  to  day  and  makes 
small  changes  in  the  value  of  at  which  the  raypaths  encounter 

the  mountain  range.  If  such  small  angular  changes  acttially  cause 
measurable  changes  of  scattered  power,  this  might  be  explained  in 
terms  of  a  scattering  coefficient  which  varies  quite  strongly  in 
the  vicinity  of  t  =  0  deg  where  Ranzi 's  measurements  were  made. 

Such  a  charac tei'istic  is  exhibited  bv  the  functions  sin^y  and 
.  4 

sin  , . 

5.  .Ament  [Ref.  37  calculates  theoretically  that  the  scattering 

coefficient  of  a  rough  surface  should  vary  in  proportion  to  sin'll 
and  that  is,  inversely  proportional  to  the  wavelength  of 

tl'.e  radio  frequency. 

The  mechanism  for  scattering  of  signals  by  terrain  has  been  intensely 
sttidied,  and  tlie  above  references  form  only  a  partial  view  of  the  whole 
subject.  The  state  of  knowledge  in  this  field,  as  of  .June  1961,  is  well 
summarized  by  a  series  of  77  abstracts  wliicli  are  given  by  Hagn,  N'eilson, 
and  Smitli  [Ref.  15,.  The  most  recent  theoretical  study  of  scattering  at 
the  time  of  this  writing  was  made  by  Hagfors  [Ref.  38  .  liis  analysis 
applies  only  for  values  of  near  90  deg,  so  the  results  are  of  little 

use  liere;  but  ho  gives  a  lucid  explanation  of  the  common  approach  to 
this  problem,  togetlier  with  a  number  of  references  more  recent  than 
those  wliicli  wore  availaltle  in  1961  to  Hagn  otal. 


C.  THE  METHOD  OF  TRE.ATING  ABSORFflOX 

One  of  the  advantages  of  tlie  technique  developed  here  is  that  it 
permits  use  of  the  most  basic  theory  as  a  source  of  information  about  the 
physical  world.  The  processing  of  the  data  often  involves  a  cons ideralil e 
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amount  of  arithmetic,  but  the  application  of  the  basic  theories  involved 
is  seldom  masked  by  complex  mathematics  or  obscure  approximations.  The 
speed  of  the  digital  computer  permits  solution  directly  from  the  basic 
laws.  Developments  such  as  the  Breit-Tuve  theoi-em  or  Martyn's  equivalence 
theorem  are  unnecessary  approximations  in  this  context. 

The  method  by  which  absorption  is  taken  into  account  should  provide 
a  clear  example  of  the  above.  .-Ibsorption  losses  lead  to  a  single, 
isolated  factor  in  the  final  result,  and  so  the  analysis  of  absorption 
can  bo  separated  from  the  general  backscalter  problem  and  will  bo  treated 
independently  here. 

The  physical  principle  involved  is  that  power  is  lost  in  proportion 
to  that  which  remains,  so  that  the  amount  which  manages  to  penetrate  an 
absorbing  region  varies  as  the  factor 


exp 

where  k  is  called  the  absorption 
traversed.  The  quantity  which  is 
the  absorption  index  a,  which  is 
equat ion 


coefficient  and  z  is  the  distance 
usually  measured  experimentally  is 
related  to  the  coefficient  by  the 


O  ^  J  k  dh 

whore  h  is  the  vertical  component  of  z.  ..\n  assumption  will  be  made 
that  the  absorbing  region  is  centered  about  a  height  of  80  km,  since  this 
is  ilio  approximate  height  of  the  D  layer,  and  only  the  nondoviative 
absorption  will  be  considered  here. 

The  above  relations  are  basic  and  can  l)e  found,  foi‘  example,  in  Ref, 
39.  There  it  is  also  shown  that  the  coefficient  of  absorption  should 
theoretically  bo  inversely  proportional  to  the  square  of  the  radio  fro- 
quenct  in  the  absence  of  a  magnetic  field.  If  the  magnetic  field  is 
taken  into  account,  this  proportion  is  approximately  true  provided  the 
operating  frequency  exceeds  the  electron  gy rof roquency  (see,  for  example. 
Ref.  40),  which  will  be  a  good  approximation  at  the  frequencies  used 
throughout  this  study. 
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In  actual  practice,  one  of  the  best  measures  of  absorption  is 
obtained  by  observing  successive  multiple  reflections  of  a  pulse  travel¬ 
ing  vertically  between  the  earth  and  the  ionosphere.  The  technique  is 
described  in  Ref.  39  where  data  from  such  an  experiment  are  given.  The 
technique  yields  good  results  but  has  two  unfortunate  limitations: 

(l)  frequency  is  restricted  to  less  than  the  vertical -incidence  critical 
frequency  of  the  most  dense  layer  because  total  reflection  must  be 
obtained;  (2)  when  frequency  is  near  the  critical  frequency  of  any  layer 
below  the  maximum,  then  the  ray  suffers  a  great  deal  of  deviative 
absorption . 

Nevertheless,  such  experimental  data  can  be  extrapolated  so  as  to  be 
useful  in  the  prediction  of  the  bel\avior  of  oblique  propagated  waves. 

In  Fig.  9  the  author  has  taken  data  from  Ref.  39  and  has  attempted  to 
make  such  an  extrapolation  from  the  cusped  curve  into  a  straight  line 
where  a  is  proportional  to  the  inverse  square  of  frequency.  These 
data,  it  should  bo  noted,  applied  in  winter  at  noontime  in  a  temperate 
zone  at  a  time  when  the  sunspot  activity  was  low.  On  the  oxti-apolatod 
■straight  line,  a  equals  0.12  at  a  frequency  of  10  Me. 

Unfortunately,  there  appear  to  be  two  different  definitions  in  cur¬ 
rent  use  for  the  moaning  of  the  index  of  absorption.  In  fact,  both 
definitions  are  used  in  Ref.  39,  The  first  definition,  as  given  above, 
leads  to  the  result  that  power  is  decreased  by  the  factor  o"^.  However, 
field  measurements  of  the  absorption  index  are  sometimes  reduced  through 
tlie  use  of  a  definition  based  on  a  power  loss  proportional  to  10 
which  facilitates  use  of  the  decibel  as  a  unit.  This  latter  definition 
is  used  in  Fig,  9, 

Lot  us  call  the  former  definition  a  and  the  latter  a  Thus 

e  10 


a 


a 


10 


10 


Comparison  of  the  value  of  from  Fig.  9  with  that  which  would 

be  obtained  if  /  k  dh  wore  measured,  would  reveal  that  two  compensating 
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errors  nearly  cancel  one  another.  That  is,  f  ^  tJh  corresponds  to 
(up  and  down)  and  since  2a^  is  only  15  percent  larger  than 


2q 

e 


e 


10 


Figure  10(a)  shows  the  relation  between  an  element  of  distance  along 
an  oblique  raypath,  ^^z,  and  an  element  of  distance  along  a  vertical 
raypath,  i'll,  such  as  was  used  in  the  determination  of  the  curve  on  Fig. 
9.  Here  it  is  seen  that 


i'iz  =  (/ll)  CSC  C 


FIG.  9.  FHEQUFNCY  VARIATION  OF  THE 
ABSORPTION  INDEX. 
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rom  which  it  follows  that  absorption  loss  is 

-a  CSC  r 
e  ^ 

Now,  since  Q  is  proportional  to  the  inverse  of  frequency  squaied, 
it  will  be  more  convenient  to  define  a  new  parameter,  E,  which  obeys 
the  relation 

2  7 

E  =  F  Q  =  1.38  X  10  for  =0.12  at  10  Me 

This  new  factor  is  not  frequency  dependent  except  near  the  gyrof rcquency , 
and  the  work  whicli  follows  is  resti'icted  to  frequencies  of  several  mega¬ 
cycles,  well  above  the  gyrof requency. 

It  is  desirable  to  program  tlie  computer  so  that  the  severity  of 
absorption  may  change  with  horizontal  range.  This  will  happen,  for 
example,  if  a  backseat  tor  radar  looks  toward  or  away  from  the  sun  across 
a  twilight,  zone.  Thus  the  factor  will  be  a  function  of  range.  A 

cuijic  equation  for  E^^  piovidcs  enough  flexibility,  so  absorption  is 
described  by  tlic  expression 

F(U)  =  E^(l  I  c^R  +  e^H^  t- 

where  E  ,  e  ,  c  and  e  can  be  specified  independently  each  time 
Ci  1  2  .1 

tile  piogram  is  run.  This  equation  has  proven  to  lie  more  than  adequate. 

In  the  equations  above,  it  is  necessary  to  specify  t^,  the  angle  at 
which  the  raypath  penetrates  the  absorbing  region.  It  may  be  assumed 
tliat  tile  raypath  from  earth  to  the  D  region  is  a  straight  lino  as  shown 
in  Fig.  10(l)).  In  the  figure,  P  is  a  known  takeoff  angle,  r^  is  the 
radius  of  the  oartli  (6370  km),  and  h^^  is  the  height  of  the  absorbing 
region  (80  km).  From  this  figure  it  is  seen  that 
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Since  E  varies  with  range,  it  will  be  necessary  to  locate  the  range 
at  which  eve.-y  rav  penetrates  the  absorbing  layer  going  up  and  coming  clown 
for  each  hop.  For  example,  if  an  element  of  backseat tering  ground  is 
illuminatcKl  by  a  two-hop  ray  and  the  energy  returns  to  the  receiver  by 
a  one-hop  ray,  then  the  calculation  of  the  recenved  power  will  include 
six  facto.-s,  one  going  up  and  one  coming  down  for  each  of  the  three  hops. 
This  sounds  difficult,  but  it  is  actually  little  more  than  a  bookkeeping 
operation  to  the  computer  which  can  simply  add  c^xponents,  thanks  to  the 
re  1  a  I  ion : 

abed  a+b+c+d 
0*C  *0  *0—0  t'lC. 

Since  E  will  be  a  slowly  varying  function  of  range,  it  is  assumed 
that  a  ray  which  reaches  a  given  range  travels  by  hops  which  are  ecpial 
and  symmetrical.  This  approximation  is  made  only  for  the  purpose  of 
evaluating  E ;  in  other  stages  of  the  analysis,  it  would  be  altogethei 
unacceptable.  Here,  its  effect  is  quite  minor  and  the  approximation 
does  lead  to  an  increase  of  computer  speed. 

It  is  necessary  to  know,  for  a  ray  which  leaves  the  ground  at  angle 
p,  how  much  ground  range  it  traverses  before  reaching  the  D  region. 
This  calculation  is  straightforward  as  shown  on  Fig.  10(b).  The  ground 
range  in  question  is  symbolized  In  this  figure  it  can  bo  seen  that 

%  =  '  e 
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It  can  also  be  seen  that 


^  =  C  -  P 

Therefore  the  instruction  given  to  the  computer  is 

%  =  ■■»<?  -  « 

Using  this  value  of  the  computer  evaluates  E(r)  at  the  two 

locations  shown  in  Fig.  10(c)  for  the  case  of  a  one-hop  ray.  Here,  the 
symbol  R  stands  for  the  total  range  reached  by  the  one-hop  ray.  For 
the  two-hop  ray,  the  computer  is  instructed  to  evaluate  E  at  the  four 
locations  shown  in  Fig.  10(d).  Similarly,  for  a  three-hop  ray,  the  com¬ 
puter  evaluates  E  at  six  places  as  shown  in  Fig.  10(e).  No  more  than 
three  hops  were  included  in  any  portion  of  the  entire  analysis,  primarily 
because  the  author  has  seldom  seen  on  any  experimental  record  the  fourth 
hop  of  ground  backscatter.  (The  fourth  or  fifth  hop  could  easily  be 
included,  if  desired.) 

Let  us  recapitulate  what  has  been  done  herer 

1.  The  absorption  has  been  programmed  so  that  it  can  be  made  to  vary 
v/ith  range  in  any  manner  matched  by  a  cubic  equation. 

2.  Absorption  has  been  corrected  so  that  the  index  is  inversely 
proportional  to  the  square  of  the  radio  frequency. 

3.  Absorption  is  calculated  to  allow  for  the  angle  at  which  the 
raypaths  penetrate  the  D  ref^ion. 

4.  Absorption  is  calculated  once  each  time  a  ray  penetrates  the  region 
going  out  to  the  backscattering  ground  and  once  each  time  the  ray 
penetrates  the  region  coming  back  from  the  same  ground,  both  when 
the  ray  goes  up  and  when  it  comes  down. 

The  above  is  simply  an  expression  of  the  fact  that  the  energy  content 
of  a  ray  decreases  exponentially  as  it  passes  through  an  absorbing  region. 
Further  complexity  could  have  been  built  into  the  program,  but  tlii  s  set  of 
choices  is  optimum  for  the  objectives  at  hand  and  refinements  can  be  added 
if  they  become  desirable.  In  the  following  analysis,  the  factor  written 

a  . , 

e  will  signify  an  instruction  to  the  computer  to  calculate  absorption 
according  to  the  criteria  outlined  here. 
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III.  ANALYSIS  OF  BACKSCATTER 


The  factors  which  must  be  taken  into  account  in  this  discussion  of 
the  ground -backseat  ter  phenomenon  are  sketched  in  Fig.  11.  It  is 
assumed  that  there  are  on  the  earth+s  surface  a  transmitter,  receiver, 
and  an  associated  antenna  which  are  to  be  simulated.  The  transmitter 
operation  may  be  described  by  specification  of  the  transmitted  power 
level  as  a  function  of  time.  In  general,  the  pulse  shape  will  not  be 
rectangular  and,  when  scattered  and  received,  it  will  be  further  modified 
by  the  smoothing  effect  of  limited  receiver  bandwidth  and  ionospheric 
dispersion. 


- ELECTRON-DENSITY  VARIATION — 

I  (H,R) 


32690 

FIG.  11.  FACTORS  WHICH  INFLUENCE  THE  AMPLITUDE  OF  GROUND 
BACKSCATTER. 


The  electron  density  can  vary  with  both  altitude  and  range.  How¬ 
ever,  the  raytracing  program  used  here  does  not  permit  electron-density 
gradients  in  the  direction  transverse  to  the  propagation  path,  since  it 
is  two-dimensional.  Such  gradients  can  bo  simulated  by  carrying  out  the 
entire  synthesis  process  a  number  of  times  for  different  azimuths 
(assuming  no  transverse  variation  within  any  single  azimuth  sector),  and 
the  final  results  can  be  added. 


SEL-65-002 


26 


At.  the  lower  edge  of  the  ionosphere,  a  certain  amount  of  the  energy 
will  be  absorbed,  and  this  absorption  will  depend  on  the  radio  frequency. 
Also,  the  amount  absorbed  depends  on  P  because  shallow  raypaths  stay 
in  the  absorbing  region  for  a  longer  distance.  Since  calculations  will 
cover  many  thousands  of  kilometers,  provisions  are  made  for  the  variation 
of  absorption  with  horizontal  range. 

Antenna  gain  will  generally  be  a  fairly  complex  function  of  frequency, 
azimuth,  and  elevation  angle.  The  exact  nature  of  its  variation  exerts 
quite  a  strong  control  over  the  backscatter  amplitude  pattern  and  so  a 
good  deal  of  effort  is  justified  in  the  simulation  of  its  effect  [kef. 

25..  The  ground  scatter  pattern  and  reflection  coefficient  have  been 
discussed,  and  the  provision  will  be  made  so  that  they  can  be  changed  as 
a  function  of  range  to  simulate  the  changing  surface  of  the  earth. 

Since  the  radio  wavelength  is  comparable  to  the  scale  of  irregularities 
on  the  earth's  surface,  these  parameters  may  vary  with  frequency.  Also, 
the  data  previously  discussed  give  a  clear  indication  that  the  scatter 
pattern  varies  as  a  function  of  the  angle  !'. 

A.  THE  FLL'X  OF  ENERGY 

Figure  12  contains  sketches  showing  raypaths,  areas,  and  dimensions 
of  interest  over  the  same  "block  "  of  the  earth.  In  part  (a)  are  shown 
the  raypaths  and  measurements  which  apply  to  the  energy  which  travels 
outward  f  i-om  the  radar  to  tlio  backsca  t  tori  ng  ground.  In  part  (b),  atten¬ 
tion  is  focused  on  the  I'ays  whicli  return  from  the  backseat  tering  ground 
to  the  radar. 

Notice  in  particular  that  it  is  not  being  assumed  that  energy  which 
travels  out  by  a  certain  raypath  must  return  by  the  same  ravpath.  To 
emphasize  this  fact,  part  (a)  sliow's  one  hop  and  part  (b)  shows  two  hops. 
However,  it  could  equally  well  be  that  both  the  transmission  and  recep¬ 
tion  propagation  paths  consist  of  the  same  number  of  hops.  For  instance, 
the  energy  could  go  out  by  one-hop-E  and  return  by  one-hop-F  or  the  same 
one-hop-E.  As  the  reader  will  see,  it  is  not  necessary  to  make  any 
restriction  and  thus,  in  this  regard,  the  program  will  be  able  to  synthe¬ 
size  a  portion  of  backscatter  which  has  previously  been  completely 
neglec  ted . 
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RAYPATH  1^ 


RADAR  IN  AREA  A^ 


FIG.  12.  ILLUSTRATION  OF  THE  BACKSCATTER  POWER-DENSITY 
CALCULATION. 

In  order  to  lend  clarity  to  the  next  few  pages,  a  convention  is 
adopted  whereby  the  subscript  T  means  that  the  subscripted  parameter 
is  associated  with  the  transmitted  rays,  and  the  subscript  R  implies 
association  with  the  received  ( backseat tered )  rays.  In  other  words, 
parameters  associated  with  Fig.  12a  will  have  a  T  subscript  and  those 
associated  with  Fig.  12b  will  have  an  R  subscript. 
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Now  consider  the  Iransmittod  rays.  Since  raypaths  are  streamlines 
in  a  flux  of  energy,  it  is  assumed  that  all  of  the  energy  which  leaves 
the  radar  within  elevation  angle  increment  within  azimuth 

angle  increment  will  travel  up  into  the  ionosphere  confined  between 

raypaths  1.^,  and  2^  and  the  two  great  circles  defined  by  This 

energy  will  refract  downward  and  strike  the  ground  within  area  .  The 
angle  increments  are  considered  to  be  small,  so  that  the  energy  striking 

is  uniformly  distributed.  Then,  since  raypaths  1,^,  and  2^  are 
calculated  by  the  raytracing  program,  it  is  possible  to  calculate  the 
energy  density  which  strikes  area  A^^  directly  from  a  sot  of  known 
pai-amoters,  including  the  antenna  gain  specified  at  elevation  angle  (3,^,. 

Area  A^  may  be  illuminated  by  other  raypaths  which  leave  the  trans¬ 
mitter  at  different  values  of  Tliis  eventuality  will  t)o  covered 

under  the  assumption  of  linearity  so  that  superposition  may  bo  used.  In 
other  words,  the  percentage  of  energy  backsca t t orod  by  the  ground  will 
tto  assumed  to  be  independent  of  the  amotint  of  enoi-gy  striking  tlio  ground 
in  that  particular  region.  This  appears  to  bo  an  oxti-omoly  safe 
assumption,  ljut  it  should  bo  noted  that  this  is  an  assumption.  Notice 
also  that  use  of  the  SRI  data  implicitly  assumes  that  the  Isacksca t  toiod 
energy  is  incoherent.  This  is  generally  accepted;  see,  for  example, 

Ref.  dl.  A  detailed  theoretical  exposition  on  this  sul),]oct  is  given  in 
Itof.  5. 

If  the  earth  were  flat,  the  width  of  ai-oa  A^^  would  simply  bo 

(it)  (A^.^),  w'hore  It  is  the  mean  distance  f i-om  the  radar  to  the  area. 

However,  the  earth's  transverse  curvature  requires  a  correction  that 

cannot  bo  ignored.  Consider  Fig.  13a  wliicli  shows  the  earth  in  cross 

section  witli  the  radar  situated  to  give  an  elevation  view.  The  distance 

from  the  radar  to  the  area  A,  is  R.  The  earth's  radius  is  r  ,  and 

1  o 

the  gooconti-ic  angle  subtended  by  R  is  a.  From  the  figure,  it  can  bo 
seen  that  R  =  r^  a  and  that  sin  a  =  Now  in  part  (b)  of  Fig.  13, 

the  same  radar  is  shown  in  plan  view.  The  distance  from  the  radar  to 
the  ground  area  appears  to  be  p,  not  R.  The  width  of  area  A^  appears 
undistortod  in  this  view  and  is  seen  to  bo  pAO .  Combining  all  those 
results,  the  width  of  the  area  is  to  be: 
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width  of  A,  =.  r  sin  — 

1  T  e  r 

o 

riio  oartli  is  not  a  sphere,  of  course,  and  so  the  selection  of 
value  for  is  somewhat  subjective.  At  the  poles  it  is  6357  km;  .n  t 

the  equator  it  is  6378;  the  average  of  these  two  is  6368;  a  sphere  whose 
circumference  is  d0,000  km  would  have  radius  6366;  the  selection  of 
6370  is  used  here  in  order  to  reflect  the  uncertainty  in  the  last  digit 
(values  from  Hof.  42). 

Looking  back  at  area  A^ ,  it  Is  seen  that  can  he  computed 

dii-ectly  from  tlio  parameters  associated  with  the  raypatlis  1^  and  2 
The  equation  for  is 


A  Ait  Alt'  r  sin  — 

1  T  T  e  r 

c 

All  of  these  parameters  can  be  obtained  from  tlie  raysots.  Asseciatod 
with  tlieso  particular  raypaths.  there  will  bo  a  certain  amount  of  energy 
injected  into  area  A^^  within  the  small  solid  angle  between  lusd 

Ait'.j,,  The  actual  solid  angle  subtended  by  A^  as  seen  from  the  radar 
via  raypaths  1,^  and  2.^  is  then  Ait'.j,  cos 


a.  Elevation  view 


AREA  A, 


b.  Plan  view 


FIG.  13,  GEOMETRY  U.SED  IN  CALCULATING  A  CORRECTION  FOR  'HIE 
SPHERICAL  EARTH. 
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A1 Lhough  actual  interest  will  be  centered  on  power  levels,  this 
analysis  will  be  simplified  if  attention  is  initially  restricted  to  the 
accounting  of  energy.  Later  a  roflnement  will  produce  the  time  dis¬ 
tribution  of  a  known  amount  of  energy,  thereby  giving  the  desired  answers 
in  Lorms  of  power.  Thus,  the  pulse  duration  will  Initially  bo  set  at 
C,  which  is  a  time  Interval  so  short  that  errors  of  the  order  of  g 
are  acceptable. 

Let  the  amount  of  energy  which  impinges  on  area  and  which 

travels  between  raypaths  1,^,  and  be  symbolized  .  Let  the  gain 

of  the  transmitting  antenna  bo  symbolized  when  it  is  evaluated  In 

the  direction  of  rayy.  L.s  i  ^  and  2,^.  Tlien  the  total  energy  trans¬ 
mitted  tov.-.'.id  A^  by  those  raypaths  is 


E,  =  GP,„G,„ 
1  T  T 


f'os 

'la 


where  P.^.  is  the  transmitted  power.  All  of  this  energy  will  strike  A^, 
except  that  reductions  must  be  made  for  the  power  lost  to  absorption  and, 
if  more  than  one  hop  is  involved,  for  the  power  lost  duo  to  the  incom¬ 
pleteness  of  i-eflection  off  the  surface  of  Die  earth.  The  latter  factor 
will  be  taken  into  account  as  follows.  It  will  be  assumed  lliat  the  earth, 
in  addition  to  backseat terliig,  specularly  reflects  (lOOp)  percent  of  the 
energy  Incident  upon  it.  The  parameter  p  is  thus  a  power  ratio  and 
will  always  bo  less  tlian  unity.  Since  p  is  piobalily  strongly  dependent 
on  the  typo  of  terrain,  the  final  progi'am  will  incorporate  an  abtlttv  to 
vary  p  with  distance  from  the  radar. 

When  is  corrected  for  alisorption  and  incomplete  reflection.  It 

represents  the  amount  of  energy  whlcli  impinges  on  ai-ea  due  to 

propavai'on  between  raypaths  l.^,  and  2^.  Many  different  propagation 
modes  may  illuminate  tlio  same  area  on  the  ground  and  ttiese  must  lie 
accounted  for  in  the  final  i-osull. 

Despite  tlio  appearance  of  Fig.  12a,  recall  tlini  area  A^  is  con¬ 
sidered  to  ho  Infinitesimal.  In  fact,  it  is  shown  ns  only  a  point  of 
backscattorlng  ground  in  Fig.  12b.  Hero  arc  shown  raypaths  1^^  and  2^^ 
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wnich  i-cprcsotiL  one  of  the  possible  propagation  modes  by  which  the 
cner-gy  wliich  strikes  area  can  return  to  the  radar. 

Suppose,  for  the  moment,  that  raypaths  1^^  and  2^^  are  actually 
known.  Then  it  can  bo  soon  by  analogy  to  A^^  that  the  area  A^  is 

A„  =  ZVIT,  A*,,  r  sin  — • 

2  II  R  o  r 

e 


However,  this  equation  is  not  directly  useful  since  the  quantity 
cannot  bo  found  in  the  raysots.  It  is  associated  with  rays  which  return 
from  the  scattering  ground,  'nd  the  raysots  describe  only  the  rays  which 
go  out  from  the  radar.  It  might  bo  possible  to  perform  a  raytracing 
which  originates  at  the  scattering  ground  area,  but  this  would  be  cumber¬ 
some  and  expensive.  Fortunately  a  more  logical  solution  exists,  which 
is  as  follows.  First,  note  that 


‘^'n  ST  **1, 


This  relation  shows  promise  because  an  opportunity  will  arise  for  the 
cancellation  of  the  factor  A'^.  Thus,  if  the  derivative  can  be  deter¬ 
mined,  the  ai'oa  A^  will  be  knovm. 

Although  this  analysis  concerns  simulation  of  a  monostatic  radar,  it 
will  be  convenient  to  digress  temporarily  and  consider  the  bistatic  sys¬ 
tem  shown  on  Fig.  14a  for  the  study  of  the  needed  derivative.  Assume 
Lliat  a  raypath  exists  which  joins  the  transmitter  and  i-eceiver,  and  that 
this  raypath  strikes  the  ground  at  angles  P  and  v  rts  shown.  Con¬ 
structed  al  the  Li-ansmitter  and  receiver  arc  an<ennas  whose  gains  arc 
equal  and  ai 'ulrartly  large.  These  theoretical  antennas  will  be  con¬ 
sidered  to  have  a  single  loljo  at  tlte  appropriate  elevation  angle,  \i 
oi-  y,  so  that  tliey  have  maximum  gain  along  the  raypath.  Facli  of  these 
antennas  has  an  effective  receiving  area,  A,  related  to  gain  by  the 
equa  t ion 
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The  areas  of  these  two  imaginary  antennas  are  normal  to  the  raypath,  by 
definition. 


.^  ANTENNA  AREA=WZ 

/ 

/ 


ANTENNA 

AREA=VZ 


FIG.  14.  RELATIO.VS  U.SED  TO  PERMIT  APPLICATION  OF  REC IPttOCITY . 


Suppose  one  has  calculated  rays  which  originate  at  the  transmitter 
end  of  this  path;  then  from  the  raysets  one  has  a  measure  of  the  spreading 
of  energy  as  it  travels  from  transmitter  to  receiver.  The  measure  most 
accessible  is  the  change  in  the  range  at  which  the  rays  return  to  earth 
when  a  small  change  in  takeoff  angle  occurs.  For  example,  suppose  a 
ray  I  lacing  at  p  =  26°  and  27°  produces  two  rays  which  land  at  the 
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earth  at  R  -  2500  and  2600  km,  respectively.  Assuming  that  these  rays 
ai’e  flux  sti'eamlines,  then  the  energy  which  leaves  the  transmitter  in 
the  1-deg  increment  of  takeoff  angle  hits  the  ground  in  a  100-km  incre¬ 
ment  of  ground  range.  The  measure  of  spreading  can  then  be  written: 

Ar/Z!:<3  =  100  km  per  degree 

If  this  parameter  is  large,  the  signal  is  weak,  and  vice  versa. 

Suppose  now  that  the  transmitter  and  receiver  are  switched  without 
any  change  in  the  antennas  or  the  ionosphere.  Then,  the  measure  of 
spreading  would  bo  AR/A'ir. 

Notice  in  particular  that  for  small  A^lt,  Ap/A'lf  —  the 

desired  cierlvativo.  We  wish  to  calculate  its  value  without  raytracing 
so  wo  turn  to  tlie  Itayloigh-Carson  reciprocity  theorem,  which  was  stated 
by  Torman  [Ref.  43]  as  follows: 

"Tf  an  electromotive  force,  K,  Inserted  In  antenna  1  causes  a 
current,  I,  to  flow  at  a  certain  point  in  a  second  antenna,  2, 
then  the  voltage,  E,  applied  at  this  point  in  the  second  antenna 
will  produce  the  same  current,  I  (both  in  magnitude  and  phase), 
in  a  short  circuit  at  a  point  in  antenna  1  where  the  voltage,  K, 
was  originally  applied.  (This)  fails  to  bo  true  only  when  the 
propaga  t  ion .  .  .  is  appreciably  affected  by  tlio  presence  of  the  earth's 
magnetic  field." 

In  other  words  the  hypothetical  switch  of  transmitter  and  receiver 
sliould  not  cause  a  change  in  the  received  signal  level  under  tlio  condi¬ 
tions  assumed  thi'oughout  this  report.  Now,  for  convenience  of  notation, 
assume  tliat  and  Atjt  are  selected  to  make  Alt  constant.  I’urthor, 

suppose  that  the  antenna  on  the  left  side  of  the  figure  has  an  eciuivalent 
area  WZ  which  is  inclined  at  angle  (i  so  that  it  casts  a  shadow  of 
length  R^  as  shown  in  part  (b)  of  the  figure.  Similarly,  let  the 
antenna  at  the  right  end  have  area  VZ  inclined  at  angle  '|/.  When  the 
transmitter  is  at  the  left  end  of  this  propagation  path,  lot  the  received 
power  (at  the  right)  bo  P  and  when  transmitter  and  receiver  are 

switched,  let  the  received  power  bo  P  .  If  the  transmitted  power  in 
’  w 

both  cases  is  P.^,,  then  all  of  the  following  relations  hold  true; 
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V 


and 


=  R  sin  \lr 

V 


W  =  R  sin  P 
w 


cos  P 

^  K  ^2 

p  - - ,  ■  ' 

V  4^  RAR  sin  \lr 

P_GA\1/  cos  \lr 

p  _  T  _ _  ^  WZ 

w  ~  4:t  RZ^iR  sin  P 

Reciprocity  indicates  that  P  =  P  .  Since  WZ  =  VZ, 

V  w 


Ap  _  sin  cos  \1<' 
A\Sr  sin  p  cos  p 


from  which 


AR  _  AR  sin  cos  ^|/ 
A’i'  ~  Ap  sin  p  cos  p 


All  of  the  elements  on  the  right  side  of  this  new  equation  are  available 
from  the  raysets,  and  so  the  derivative  can  be  calculated  without  further 
ray  tracing. 

Returning  to  the  monostatic  backscatter  radar,  it  is  seen  tliat  an  R 
subscript  must  be  put  on  the  R's,  'I^'s,  and  P's  in  the  above  equations 
to  signify  that  the  values  are  selected  from  the  rayset  that  describes 
the  returning  raypath  (which  may  be  different  from  the  outgoing  raypath). 
Thus,  the  area  is  calculated  according  to  the  relation 


=  (A^tp  A-t^; 


'^c^R^  sin  cos 

R  ^R  ^R  .  R 

r  sin  — 

op„  sin  p  cos  P  e  r 

V  R  R  R  e, 


It  is  not  necessary  to  worry  about  the  algebraic  sign  in  such  expressions 
as  these  because  they  will  all  bo  factors  in  an  expression  for  the 
absolute  value  of  the  received  power,  which  is  known  to  be  positive. 
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Since  the  amount  of  energy  striking  area  is  known,  the  portion 

of  this  energy  which  returns  toward  the  radar  can  be  calculated  by 
application  of  the  ground  backscatter ing  coefficient,  K(\lr).  Measure¬ 
ments  such  as  those  of  SHI  were  taken  by  a  single  monostatic  radar  in 
an  aircraft  and  therefore  they  apply  only  to  the  case  where  the  angle  of 
incidence  is  equal  to  the  angle  of  scatter.  In  the  problem  being 
analyzed  here,  such  a  i-elation  will  not  always  hold  true.  For  example, 
compare  parts  (a)  and  (b)  of  Fig.  12  where  it  can  be  seen  that  the 
angles  and  \l/j^  are  not  equal.  In  fact,  they  must  be  unequal  for 

all  cases  except  the  commonly  assumed  one  in  which  the  energy  travels 
to  and  from  the  backseat tering  ground  via  the  same  raypath.  However, 
Ihei'e  is  no  experimental  information  known  to  this  authoi'  which  indicates 
the  beliavior  of  the  backscattering  coefficient  when  and  are 

different,  so  it  seems  reasonable  to  calculate  their  average  value  and 
to  use  that  as  the  argument  of  the  backscattering  coefficient. 

Each  time  any  ray  penetrates  the  absorbing  region,  there  will  be  an 

energy  loss  which  must  be  calculated  by  the  program.  This  loss  will  be 

accounted  for  by  multiplication  of  the  total  energy  by  a  constant  whose 

value  is  between  0  and  1.  This  constant,  calculated  by  the  method 

Q'j’ 

previously  described,  will  be  symbolized  e  for  the  outgoing  (trans¬ 
mitted)  laypatli.  Similarly,  the  attenuation  of  tlie  incoming  (received) 
ray  pa  til  will  bo  symbolized  e 

The  factors  which  account  for  the  incomplete  reflection  off  the 
ground  of  the  transmitted  and  received  rays  will  be  symbolized  and 

Lot  E  represent  the  amount  of  energy  which  returns  to  area  A  . 
Then,  the  above  arguments  can  be  combined  to  show  that 


^2  = 


eP.j,G.j,  AC.J,  cos  fl.j, 

4ji 


t-T  ® 


^T  ^ 


%  ^"^R  4 

2k 


‘■r  ® 


Qr 
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In  this  expression  the  2:1  appears  in  the  last  denominator  because  the 

definition  of  the  function  K(\lf)  involves  a  power  comparison  against 

an  ideal  radiator  which  transmits  isotropically  into  a  hemi sphere.  If 

a  complete  sphere  had  been  used,  the  factor  would  have  been  4n  and  the 

parameter  K  would  have  been  more  closely  analogous  to  a  gain  function. 

Let  the  receiving  antenna  have  an  equivalent  area  symbolized 

The  use  of  an  equivalent  area  such  as  this  implies  that  the  area  is 

oriented  in  space  perpendicular  to  the  direction  of  the  incoming  rays  to 

be  received.  Thus,  when  comparisons  are  made  between  and  the  area 

A  which  lies  flat  in  the  earth,  a  correction  must  be  made  because  A 
^  2 

is  not  perpendicular  to  the  raypaths;  the  raypaths  make  an  angle 

with  area  A  . 

2 

In  a  practical  situation,  the  antenna  used  will  transmit  and  receive 

only  one  of  the  two  possible  polarizations.  Suppose,  for  example,  that 

the  antenna  transmits  and  receives  horizontal  polarization  only.  During 

both  the  transmission  and  reception  parts  of  the  ray  propagation,  the 

plane  of  polarization  will  rotate.  On  the  average,  only  half  of  the 

returning  energy  will  exist  in  the  form  of  horizontal  polarization.  Thus, 

a  correction  factor  of  1/2  must  be  inserted  to  account  for  the  fact  that 

half  the  energy  is  undetectable  due  to  polarization  rotation. 

Returning  to  the  analysis,  let  E  represent  the  total  received 

R  "  ■ 

energy.  The  above  reasoning  then  indicates  that 

E 

R  2  A^  sin 
2  2  2 

Now  c  =  FA  and  G  =  AnA/A  ,  so  A^^  =  c  C^/4-nF  .  Also  G.^,  G(P,j,) 

and  G  =  g(P  ).  With  this  notation  it  is  finally  possible  to  express 
1\  K 

the  received  energy  in  terms  of  known  parameters  as  follows: 
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E„  =  - 


4nF 


'  ’’’'-4^)  j;  1%  "r  ♦r)  '’’t  “(Pt) 


/g%  »tn  ♦„  °°=  ♦„  _g_ 

^  R^  \dp„  sin  cos  ®  "  r 


A^T  ^T 

4:i 


c„  e  •  e 
“^R 


B.  POWER;  THE  TIME  DISTRIBUTION  OF  ENERGY 

At  this  point  it  is  necessary  to  include  time  in  the  analysis.  The 
above  equation  gives  the  total  amount  of  energy  received  by  a  selected 
combination  of  modes  for  a  very  short  transmitted  pulse  length.  Now  the 
energy  must  be  traced  in  both  space  and  time,  a  process  which  is  most 
easily  visualized  on  a  three-dimensional  graph  whose  coordinates  are 
range,  time,  and  power.  Such  a  plot  is  shown  in  Fig.  15  which  sketches 
the  distribution  of  energy  on  the  ground  after  an  ionospheric  reflection 
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METER^  TRANSMITTED  PULSE  SPACE-TIME  DISTRIBUTION 

OF  ONE- HOP  ENERGY 
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FIG.  15.  THE  DISTRIBUTION  OF  ENERGY  IN  SPACE  AND  TIME  AFTER  ONE  HOP. 
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similar  to  the  one  described  by  the  family  of  raysets  printed  on  Fig.  6. 
Notice  that  volume  on  this  particular  three-dimensional  graph  will 
represent  energy  per  meter,  but  it  will  not  be  necessary  to  calculate 
the  energy  quantitatively  because  the  total  received  energy  has  already 
been  expressed  in  terms  of  known  parameters.  Consequently,  this  dis¬ 
cussion  seeks  only  to  clarify  qualitatively  the  time  distribution  of 
received  power. 

The  space  and  time  distribution  of  energy  shown  on  Fig.  15  will  be 
affected  by  such  things  as  the  antenna  gain  and  the  ionospheric  absorp¬ 
tion.  Ionospheric  focusing  and  defocusing  also  play  a  strong  part  in 
many  circumstances,  but  this  simple  Chapman  layer  shows  only  the  defocusing 
of  upper  rays,  which  leads  to  the  low  power  level. 

Now  consider  a  single  scatterer  taken  along  the  lower-ray  portion  of 
the  energy  in  Fig.  15.  This  scatterer  is  shown  in  Fig.  16  where  it  now 
assumes  the  role  of  a  transmitter.  The  energy  scattered  back  toward  the 
radar  assumes  a  space-time  distribution  like  that  in  Fig.  15. 


LOWER- RAY  ENERGY  RETRANSMITTED 
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FIG.  16,  THE  DISTRIBUTION  IN  SPACE  AND  TIME  OF  THE  ENERGY  FROM  A  SINGLE 
SCATTERER. 
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strictly  speaking,  the  energy  would  go  both  to  the  right  and  left 
on  these  figures,  but  this  fact  will  be  ignored  since  it  would  serve  no 
purpose  here.  Also,  no  attempt  is  being  made  to  depict  actual  relative 

power  levels  from  one  hop  to  the  next. 

Since  the  radar  is  located  only  at  range  zero,  it  would  receive 
energy  only  at  the  two  time  delays  where  the  zero  range  line  intersects 
the  scattered  energy  distribution  on  Fig.  16.  Actually,  in  this  case 
the  single  scatterer  would  receive  energy  by  two  modes  and  would  scatter 
it  back  by  two  modes,  with  a  consequent  reception  of  four  different 
elements  of  energy  at  three  different  time  delays. 

Figure  17  is  similar  to  Fig.  16  except  that  it  describes  the  action 
of  a  group  consisting  of  all  scatterers  in  a  strip  of  terrain.  If  only 
lower  rays  are  considered,  this  energy  scatters  back  toward  the  radar 
where  it  lands  in  a  space-time  distribution  which  appears  to  be  a  solid 
bar  on  the  figure.  The  intersection  of  the  zero-range  plane  with  this 
solid  bar,  as  shown  in  red,  is  therefore  a  pulse  of  energy  which  is  very 
nearly  rectangular.  The  power  level  in  this  pulse  will  be  symbolized 
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FIG.  17.  THE  TIME  DISTRIBUTION  OF  ENERGY  RECEIVED  FROM  ONE  RANGE 
INCREMENT  VIA  SELECTED  MODES:  A  MODESET . 
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P, ,  and  its  beginning  and  ending  times  will  be  called  T  and  T 
A  L  H 

respectively.  These  three  parameters,  taken  as  a  group,  fully  describe 
the  energy  which  is  received  for  a  particular  combination  of  modes,  and 
so  the  threesome  will  be  termed  a  "modesei." 

In  order  to  correlate  the  discussion  which  led  to  Fig.  17  with  the 
discussion  of  energy  centered  around  Fig.  12,  a  number  of  the  parametert. 
shared  by  these  two  figures  are  drawn  in  on  Fig.  18. 
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FIG.  18.  KNOW  MEASURES  OF  THE  ENERGY  DISTRIBUTION. 


From  the  arguments  given,  it  can  clearly  be  seen  that 


\  = 


T  +  T 
IT  IR 


and 


2T 


+  T 


2R 


Thus,  the  energy  bundle  E  will  arrive  in  the  time  interval  (t  ,T  ). 

R  L  H 

Under  the  assumption  of  a  uniform  distribution  of  energy  within  a  flux 
tube,  it  follows  that  the  power  level  is 
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P.  =  - - —  ,  T  <  t  <  T„ 

A  -  Tl  ■  L  H 


P^  =  0 


otherwi se 


C.  SUMMATION  IN  AZIMUTH  AND  ELEVATION  ANGLES 


Next,  it  is  necessary  to  carry  out  a  summation  in  azimuth  and  ele¬ 
vation  to  account  for  all  possible  flux  tubes  to  and  from  the  back- 
scattering  ground.  The  summation  in  azimuth  is  relatively  easy 

because  the  only  expression  which  varies  in  azimuth  is  the  antenna  gain 
so  that  the  summation  can  bo  changed  to  an  integration.  The  influence 
of  antenna  gain  and  its  method  of  variation  constitute  a  rather  extensive 
subject  and  so  the  entire  matter,  including  the  integration  in  azimuth, 
is  discussed  in  Chapter  IV. 

Although  the  summation  of  results  with  respect  to  azimuth  can  be 

changed  into  an  integral,  the  same  is  not  true  of  the  summation  with 

lespect  to  elevation  angle  p,  because  the  energy  is  a  multiple-valued 

function  of  p.  The  complication  arises  for  the  following  reason: 

Suppose  three  different  flux  tubes  reach  the  same  spot  on  the  eartli's 

sui-face,  each  with  a  different  p.  Then,  energy  contributions  will 

exist  for  each  possible  combination  of  these  three  paths,  and  it  will  be 

,  2  ^ 

necessary  to  calculate  (.3  =9^  modesets. 

The  summation  with  respect  to  elevation  angle  will  actually  be  carried 
ou 1  in  terms  of  the  parameter  range  R.  Tlie  surface  of  the  earth  will  be 
considered  divided  into  a  series  of  concentric  rings  as  defined  by  the 
intersections  with  the  ground  of  the  individual  raysets.  All  energy 
which  strikes  each  ring  will  be  accounted  for  in  such  a  manner  that  every 
possible  propagation  mode  is  included. 

In  the  section  which  discusses  gain,  it  will  be  shown  that  the  term 
A*!*  can  be  eliminated  from  the  expression  for  provided  gain  is 

leplaced  by  a  scalar  multiplier  g^  operating  on  the  function  g(p). 

The  factor  g^  is  a  result  of  the  integration  in  azimuth. 

The  computer  will  calculate  the  value  of  E^  and  P^  for  every 
possible  combination  of  modes  that  can  lead  to  an  element  of  the  received 
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backscattered  power.  It  will  add  each  of  the  resulting  power  levels  along 

a  time  base  between  the  values  T  and  T  .  When  this  process  is  finished 

L  H 

there  will  be  some  distribution  of  power  as  a  function  of  time,  which 
would  represent  the  backscatter  signal  if  the  transmitted  pulse  were  of 
short  duration  e  and  if  errors  were  permitted  of  the  order  of  e.  How¬ 
ever,  the  actual  transmitted  pulse  has  a  measurable  duration  and  a 
definite  shape.  Also,  after  this  pulse  has  traveled  out  and  been  reflected 
off  a  scatterer,  it  will  be  received  with  a  shape  that  is  further  modified 
by  the  influence  of  the  receiver  bandpass  characteristics.  These  processes 
are  simulated  by  the  use  of  a  weighted -summation  scheme  within  the  com¬ 
puter.  The  analogous  procedure  in  the  domain  of  continuous  variables  can 
be  expressed  by  use  of  the  convolution  integral  as  indicated  below. 

Let  the  transmitted  pulse  shape,  as  modified  by  the  receiver  bandpass 
function,  be  symbolized  f(t)  whore  t  is  time.  Similarly,  let  the 
power  level  which  was  calculated  for  pulse  duration  e  be  symbolized 
Pg(t).  The  actual  received  power  will  be  symbolized  and  is 

related  to  known  quantities  by  the  convolution  integral, 

p^(0  =  f  p.(0  f(t  -  0  d| 

-^0 

In  the  discrete  number  system,  such  as  is  used  in  the  computer,  the 
pulse  shape  f(t)  is  given  as  a  sequence  of  numbers.  Similarly,  the 
power  distribution  ( L )  is  given  as  a  sequence  of  numbers  which  repre¬ 
sent  the  powei'  level  sampled  at  equal  increments  of  time,  each  of  wliich 
equals  the  "infinitesimal"  duration  .  Anntlier  way  of  visualizing  this 
is  to  consider  tliat  (he  transmitter  power  level  changes  its  value  aljruptly 
at  each  increment  of  time  .  As  with  any  computer  solution,  the  method 
leads  to  errors  in  the  Lime  distribution  of  energy  on  the  order  of  e. 

Thus  llie  value  of  e  must  bo  selected  small  for  good  results  and  large 
for  quick  computation;  there  are  a  number  of  other  parameters  which  must 
similarly  bo  cliosen  as  a  compromise  between  conflicting  requirements 
( o ,  g .  ,  Ab)  . 

The  underlying  assumption  hero  is  that  the  elements  of  the  backscat- 
Lered  signal  add  incoherently.  If  each  element  added  coherently,  it 
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would  be  possible  to  carry  out  the  analysis  using  a  convolution  integral 
on  the  received  voltage.  If  the  actual  situation  lay  somewhere  between 
coherence  and  incoherence,  the  convolution  approach  probably  would  not 
work*  because  it  is  a  mathematical  model  based  on  linearity  [Ref.  44]. 

This  essentially  ends  the  analysis  of  ground  backscatter,  because 
all  the  factors  in  the  expression  can  bo  found  by  appropriate  cal¬ 

culations  starting  from  the  rayset  information.  Now,  it  will  bo  neces¬ 
sary  to  carry  out  the  process  of  synthesis  wliich  in  effect  will  close  a 
loop  wlien  it  arrives  at  the  expression  for  P^. 


-y- 

The  possibility  exists  that  satisfactory  results  might  bo  obtained  using 
voltage  raised  to  some  power  between  1  and  2. 
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IV.  THE  EFFECT  OF  ANTENNA  GAIN  PATTERNS 


Since  the  objective  of  this  study  is  synthesis  of  backscatter  as 
it  might  be  received  by  some  particular  radar,  it  is  appropriate  to 
consider  some  sources  of  information  concerning  the  gain  pattern  of  the 
antenna  G(p,<t!).  The  most  commonly  available  information  consists  of 
graphs  showing  the  gain  as  a  function  of  elevation  angle  g(P)  and  as 
a  function  of  azimuth  Cj(<t').  The  actual  measurement  of  the  gain  of  an 
antenna  is  an  expensive  and  time-consuming  project,  so  it  is  seldom 
carried  out  unless  the  antenna  is  part  of  a  system  where  calibration  is 
of  sufficient  importance  to  justify  the  effort  [kef.  45],  In  most  cases, 
the  functions  g(p)  and  q('t)  are  calculated  on  the  basis  of  theory, 
empirical  formulas,  or  measurements  made  on  a  small  model.  Such  results 
are  usually  of  acceptable  accuracy  as  a  description  of  the  performance 
of  high-frequency  antennas. 

In  the  equation  which  gives  the  energy  contribution  of  a  single 
modeset,  the  gain  is  entered  in  the  calculation  in  an  expression  of  the 
following  form: 


whore  G.^,  is  the  gain  of  the  transmitting  antenna  and  G  is  the  gain 
of  the  receiving  antenna.  Throughout  this  report  it  is  assumed  that 
there  are  no  ti-ansverse  lilts  in  the  ionosphere,  so  it  is  permissible  to 
drop  the  subscripts  on  it  because  0.^,  =  0  ;  that  is,  all  propagation 

is  by  great-circle  routes,  and  a  ray  which  starts  out  in  one  azimuth  will 
return  along  the  same  azimuth.  Also,  most  backscatter  radars  make  use  of 
a  single  antenna  so  that  the  subscripts  on  the  G  can  bo  eliminated. 
However,  it  is  not  possible  to  ignore  the  subscripts  on  P  because  a 
large  amount  of  energy  actually  leaves  the  radar  at  takeoff  angle 
and  travels  to  and  from  the  backseat toring  ground  by  different  propagation 
modes.  It  will  thus  arrive  back  at  the  radar  at  an  angle  of  arrival  B 
which  is  not  equal  to  P.^, . 


_  V  \' 


R  ^ 


^  .j,  .j, 


45 


SEL-65-002 


Now  suppose  that  gain  G  can  be  represented  by  the  product: 


Jr 


G(p,t)  =  ug(p)  q(«>) 

where  g  and  q  are  functions  of  arbitrary  size  but  appropriate  shape 
and  u  is  a  scalar  multiplier,  This  equation  implies  that  gain  is  a 
function  which  is  separable  in  terms  of  the  azimuth  and  elevation  vari¬ 
ables,  a  condition  which  is  not  always  true.  Nevertheless,  it  is  an 
assumption  whlcli  is  made  so  commonly  that  it  is  not  widely  recognized 
to  be  an  approximation.  There  are  probably  two  reasons  foi'  this: 

(l)  Tile  actual  gain  is  relatively  difficult  to  detei’mine,  as  was  men¬ 
tioned  above,  and  (2)  it  is  difficult  to  make  any  useful  calculations 
using  information  which  is  bettei'  than  that  which  can  bo  represented  by 
the  assumption  of  separable  variables.  Most  users  of  gain  information 
are  interested  in  either  g  or  q  but  seldom  have  enough  information 
about  other  system  variables  to  make  use  of  both  functions  simultaneously. 
Thus,  if  given  G,  they  would  find  it  necessary  to  calculate  either  g 
or  q  for  use  in  solving  the  pi’oblems  at  hand. 

One  of  the  most  gross  departures  of  gain  from  a  function  which  can 
be  represented  by  separable  variables  exists  in  the  form  of  side  lobes 
which  arc  centered  about  elevation  angles  different  from  the  elevation 
angle  of  the  main  antenna  lolie.  In  backscatter  synthesis,  it  is  fortu¬ 
nately  practical  to  simulate  this  situation  closely  while  still  making 
use  of  the  assumed  separable  variables.  Because  the  elements  of  the 
backseat tered  signal  add  in  a  linear  fashion,  the  superposition  principle 
will  allow  separate  calculations  to  be  made  for  each  lobe  which  cannot 
bo  simulated  by  a  product  such  as  G  =  ugq.  Separate  backscatter  can  be 
calculated  for  eacli  such  lobe  and  the  resulting  power  levels  can  be 
added  graphically  in  a  simple  manner.  To  state  this  mathematically,  it 
would  be  said  that 

G(f5,<t)  ^  u^g^(p)  q^{c)  +  • 


SEL-65-002 


46 


where  the  functions  and  constants  are  selected  so  that  each  term  repre¬ 
sents  one  or  more  of  the  actual  antenna  lobes.  The  number  of  terms  v/hich 
can  thus  be  used  are  limited  only  by  the  patience  of  the  user. 

From  the  above  arguments,  it  can  then  be  said  that  gain  enters  into 
the  backscatter  calculation  in  the  following  manner; 


E 


R 


\  \  2  , 

_  _  f(PT..Pp)  »  g(PT.)  g(Pp)  q“(«)  AS  At 


Notice  that  the  summation  in  C  involv'cs  only  the  single-valued 
2/  \ 

function  q  (Cj.  Thus,  this  summation  can  bo  transformed  into  an  inte¬ 
gral  and  taken  outside  the  summation  in  p,  as  follows; 


E 


R 


g(P.j.)  g(P^) 


This  operation  permits  the  integration  to  bo  done  analytically.  The 
bracketed  quantity  is  thus  reduced  to  a  scalar  multiplier. 

Gain  is  defined  as  the  ratio  of  (the  power  transmitted  in  the 
diin-ction  p ,  <5  )  to  (the  power  which  would  be  transmitted  in  that  direc¬ 
tion  if  the  antenna  wore  lossless  and  isotropic).  In  this  development 
t  lie  transmitted  power  will  be  considered  as  the  total  radiated  power  and 
llius  the  antenna  may  bo  assumed  lossless. 

From  the  above  definition,  it  can  bo  seen  tliat  conservation  of 
energy  requires 


cos  p  cIB  dl 


4n 
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which  can  be  written 


(u) 

■  /-t/2 

1  g{P)  cos  P  dp 

(p  q( 4 )  d4 

/p=o 

and  this  relation  will  hold  true  provided  ii  is  calculated  according 
to  the  relation; 


r 

/  b(P)  cos  P  dp 

®  q(4)  d4 

yp=o  j 

Looking  at  tlie  way  in  which  gain  enters  the  backscatter  calculation, 
it  can  t)c  seen  that  the  mathematics  will  take  on  its  simplest  form  if  w’o 
choose  to  simulate  gain  according  to  the  relation; 


G{P.t) 


'(4  )  d<t 


g(P) 


The  quantity  in  brackets  is  merely  a  calculated  constant.  II  tlie  value 
of  u  is  inserted,  the  factors  can  bo  rearranged  in  the  following  form 


G  =  g.B(f5); 


q“{4  )  d4 

^i  " 

r 

.tt/2 

(p  q(c)  d4 

1  g(P)  cos  p  dp 

yp=o 

Hero  it  can  be  seen  that  the  first  factor  in  the  expression  foi  g^  is 
a  constant  calculated  onl>  from  the  azimuthal  variation  of  gain.  The 
second  factor  is  calculated  only  from  the  variation  of  gain  with  eleva¬ 
tion  angle. 

It  is  interesting  to  note  that  the  first  factor  represents  the  only 
influence  of  q(4)  on  the  entire  backscatter  calculation.  Thus  the 
behavior  of  ttiis  factor  is  of  exceptional  interest  and  it  will  be  discussed 
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at  some  length,  bolli  to  illusti-ato  liow  azimuthal  gain  affects  back- 
scatter  and  to  derive  a  formula  for  mathematical  simulation. 

A.  METHODS  OF  PROGRAMMING  GAIN  FUNCTIONS 

In  the  computer  program  which  makes  synthetic  backscatter  records, 
it  is  sufficient  to  have  a  scalar  multiplier  g^  together  with  the 
function  g(p).  Ttie  function  is  inserted  into  the  program  bv  two  liif- 
ferent  methods,  eacli  of  whicli  is  well  suited  to  a  particular  form  of 
input  information. 

In  one  approach  tire  computer  is  actually  given  a  rattier  complex 
analytic  equation  which  has  been  designed  so  tliat  a  largo  number  of  dif¬ 
ferent  shapes  of  functions  can  lie  syntiiesized .  Then,  I  lie  user  can  insert 
various  constants  into  the  equation  liy  moans  of  punched  cards. 

Anotlier  separate  computer  program  was  designed  to  support  this 
operation  liv  carrying  out  i  lie  integration  wliich  loads  to  the  second 
factor  in  g^  and  at  tlu'  same  time  making  a  graphical  plot  of  tlie  func¬ 
tion  g(t'').  This  form  of  computation  i  .s  well  suited  to  studies  where 
one  is  most  intorosted  in  other  parameters  and  only  cares  that  tlio  gain 
is  easily  pi-ograminod ,  fairly  realistic,  and  consistent  with  the  known 
variation  of  a  gain  function. 

The  formulas  irsed  for  t  lie  gain  model  can  be  solectofi  from  the  follow¬ 
ing  array  which  is  built  into  tlie  computer  program: 


G(p) 


r 

cos”  a(p  t-  r^) 

1 .0 

■ 

sin  [a  sin  b  (f  -i  )  ] 

a  sin  b  ( p  ^  ) 

sin  [k  sin  (p  i-  p^) 

sin  a  (p  +  p^  )” 

I^(P  +  Pj  ) 

_ 

y 

In  this  formula  g.  is  tlie  constant  calculated  according  to  the  criteiia 
mentioned  previously.  The  paiametcr  is  used  to  raise  or  lower  the 


r 
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beam. 


The  fii-st  array  represents  the  free-spaco  pattern  of  the  antenna. 
The  second  array  represents  the  effect  of  a  perfectly  conducting  eai th 
on  the  free-space  pattern  provided  that  the  lower  element  in  the  brackets 
is  selected.  If  unity  is  selected,  then  the  shape  of  the  antenna  gain 
pattern  is  determined  only  by  the  square  of  the  element  selected  from 
the  first  array.  The  parameter  k  allows  for  antenna  height  h  through 
the  relation; 


which  is  derived  for  horizontal  polarization  and  for  vortical  polariza¬ 
tion  under  the  nrewster  angle  in  Ref.  46. 

The  elements  in  the  first  factor  were  selected  somewhat  subjectively 
by  the  author  in  order  to  make  maximum  use  of  tlie  constants  n.  a,  b, 
while  at  the  same  time  including  as  special  cases  some  of  the  theoretically 
derived  oquatiotis  for  froo-spaco  antenna  patterns. 

Perhaps  of  more  interest  is  the  programming  technique  which  has  been 
termed  "itemized  gain.”  Tliis  approach  is  better  suited  to  the  use  of 
specific  experimental  or  theoretical  gain  functions.  As  the  name  implies, 
the  ftinction  g{t:')  is  fed  into  the  computer  by  moans  of  a  large  numlter 
of  sample  points  along  tlie  function.  Tlie  computer  ilien  performs  inter- 
polatioti  between  tlie  data  points  to  determine  tlie  value  of  tlie  function 
at  intermediate  values  of  the  takeoff  angle.  The  integration  of  the 
function  g(tt)  must  bo  done  so  often  that  it  lias  Iteen  found  convenient 
to  invent  a  graphical  technique  wliich  permits  use  of  a  planimeter. 

Figure  19  shows  the  graph  devised  for  this  purpose.  The  horizontal  axis 
is  foreshortened  by  the  factor  cos  p  and  the  vertital  axis  is  labeled 
so  that  gain  can  be  plotted  in  decibels,  but  the  vertical  axi.s  scale  i.s 
linear  in  power.  Thus,  if  the  itemized  gain  points  are  plotted  on  this 
giaph  and  the  planimeter  is  run  along  the  curve  defined  by  the  points, 
the  area  measured  will  be  proportional  to 

at/2 

/  g(p)  cos  p  dp 
-/pzrO 
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A  normalization  area  can  be  measured  on  the  graph  and  the  ratio  of  the 
two  readings  then  used  to  calculate  g^ . 


B.  AZIMUTHAL  GAIN 

The  first  factor  in  the  expression  for  g^  warrants  further  dis¬ 
cussion  because  it  represents  the  only  effect  of  azimuthal  gain  on  the 


FIG.  19.  A  CHART  DEVISED  TO  PERMIT  USE  OF  A  PLANIMETER  FOR  THE 
INTEGRATION  OF  GAIN. 


51 


SEL-65-002 


POWER  RATIO 


>■ 


backscatter,  provided  the  underlying  assumptions  are  valid.  This 

excludes,  for  example,  the  effect  of  side  lobes  which  can  be  calculated 

by  a  separate  run  of  the  same  computer  program.  For  convenience,  define 

the  numerator  and  denominator  of  the  expression  to  be  q  and  q  , 

n  d 

that  is. 


^  q^($)  dt  :  ^  n('^) 


d4 


The  ratio  of  those  factors  will  be  calculated  for  various  widths  and 

shapes  of  the  azimuthal  gain  pattern  illustrated  o!i  Fig.  20.  Here  are 

shown  three  functions  q('l-)  of  arbitrary  width  and  varving  shape.  The 

first  example,  function  q  ,  has  a  constant  value  between  two  azimuths 

a 

and  is  zero  elsewhere:  it  might  be  termed  rectangular.  Shown  on  the 

drawing  is  the  3-db  beamwidth,  symbolized  i,,,,  ,  and  the  null-to-null 

3db 

beamwidtli,  symbolized  C^.  Each  of  the  three  examples  on  the  figure 
was  selected  so  that  the  first  null  falls  at  azimuth  ±r./2n,  and  thus 
tile  total  width  of  these  patterns  between  the  nulls  has  the  valtic  of 
4  =  Tt/n.  Similarly,  lot  =  ji/l.  Going  back  to  Fig.  20a,  it  is 

seen  that 


m  n 

mn 

”  n 

~  n 

\  /  n  \ 

=  (m  /- 

■■=  /  -  •  n  =  0 . 564 

’  V  n 

/  \  mn  f 

V  Jt 

In  this  case,  L  =  n  and  so  q  /q,  =  0.564  L  . 

n  d  v 

Now  considei'  case  (b)  where  qj^(4)  is  a  sinusoid.  Here, 


2  m  It 
“^n  ^  ^rT 


2m 

'’d  "  T 


\  ^d 
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(c) 


32683 

FIG.  20.  HYPOTHETICAL  SHAPES  USED  FOR 
CALCULATING  THE  EFFECT  OF  THE 
AZIMUTHAL  VARIATION  OF  GAIN. 


Since  the  3-db  points  occur  at  <$> 


±r>/3,  L 


3n/2  and  therefore 


q  /q,  =  0.513  /h 
n  d  '' 
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In  case  (c)  the  function  is  a  squared  sinusoid,  probably  the  most 
realistic  of  the  three  patterns  shown.  i^ere, 


q  /q^  =  0.690  Jn  =  0.488  J~h 
n  d  '' 

These  i 
Null  ■ 

3-db 


esults  can  be  summarized  as  follows: 
to-null  beamwidth  =  Ai 
beamwidth  =  t\/i, 


Shape  of  Gain 

Function,  q 

Value  of  q  /q., 
n  d 

Multiplier  of  ^ r\ 

, 

Value  of  q  q , , 
n  d 

Multiplier  of 

Square 

0.564 

0.564 

Sine 

0.626 

0.513 

Sine^ 

0.690 

0.488 

Variation 

(^) 

22 

15 

Two  interesting  conclusions  follow  from  the  argument  given:  (l)  A 
high  directivity  in  azimuth  leads  to  a  stronger  backscatter  signal  in 
approximate  proportion  to  the  inverse  of  the  square  root  of  llie  beam- 
width,  and  (2)  the  shape  of  the  azimuthal  gain  function  does  not  have  a 
very  strong  effect  on  tlie  backscatter  signal  strongtli  for  a  given  beam- 
width.  In  fact,  we  liave  shown  that  the  shape  can  undergo  the  drastic 
changes  shown  on  Fig.  20  and  Lliat  the  backscatter  signal  strength  will 
vary  by  only  30  percent  provided  the  beamwidth  remains  constant  in  terms 
of  the  pai-ametei'  L.  [The  variation  would  be  reduced  to  8  percent  if 
(  remained  constant.] 

Finally,  this  result  is  useful  because  it  opens  up  a  number  of  pos¬ 
sibilities  for  the  mathematical  sin.ulalion  of  the  effect  of  azimuthal 
gain  patterns.  Uncertainties  in  tie  experimental  measurement  of  back¬ 
scatter  signal  levels  are  usually  far  in  excess  of  30  percent  and  prob¬ 
ably  should  be  measured  in  orders  of  magnitude.  Thus,  if  the  only 
information  available  is  a  number  giving  the  beamwidth  of  the  azimuthal 
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pattern,  a  good  approximation  can  bo  made  to  the  value  of  the  parameter 
by  using  either  0.69  or  0.5  y/h,  depending  on  which  of 

these  factors  is  known.  Of  course,  if  the  actual  azimuthal  pattern  is 
known,  the  integration  can  be  carried  out  and  this  approximation  may  be 
eliminated . 

C.  SUMMAuy  OF  GAIN  DISCUSSION 

The  effect  of  antenna  gain  on  the  backscatter  signal-strength  dis¬ 
tribution  has  been  analyzed  in  some  detail.  A  single  transmitting- 
receiving  antenna  is  assumed,  although  this  argument  could  bo  readily 
extended  to  cover  use  of  separate  antennas.  The  gain  is  mathema I ically 
broken  down  into  the  sum  of  several  functions,  each  of  which  can  bo 
represented  by  a  function  which  is  separable  in  terms  of  the  variables, 
azimuth  and  elevation  angle.  This  procedure  results  in  a  marked  sim¬ 
plification  of  the  mathematics  and  yet  utilizes  all  the  information 
which  is  normally  available  concerning  antenna  gain  patterns. 

The  effect  of  the  elevation-angle  variation  in  gain  is  thus  separated 
out  from  the  main  problem  and  two  different  techniques  are  shown  for 
simulating  this  pattern  during  computer  operations.  A  graphical  method 
of  carrying  out  the  needed  integration  i.s  also  given. 

Backscatter  signal  stif'ngth  is  shown  to  vary  approximately  in  inverse 
proportion  to  the  square  root  of  azimuthal  beamwidth.  The  shape  of  the 
azimuthal  gain  pattern  does  not  have  a  marked  effect  on  the  backscatter 
signal  strength  for  a  given  Ijeamwidth,  assuming  that  the  ionosphere  and 
terrain  do  not  vary  in  the  direction  transverse  to  the  propagation  path. 
Such  variations  can  bo  simulated  by  moans  of  separate  computer  calcula¬ 
tions  foi-  each  azimuth  and  subsequent  addition  of  the  results,  by  virtue 
of  the  fact  tliat  signals  add  incoherently  in  such  a  manner  that  super¬ 
position  can  bo  applied. 
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V.  THE  TECHNIQUE  OF  DIGITAL  SYNTHESIS 


This  chapter  describes  the  operation  of  the  computer  program  which 
synthesizes  backscatter,  together  with  comments  which  illustrate  the 
connection  between  some  of  the  processes  carried  out  digitally  and  the 
corresponding  backscatter  mechanisms.  Input  data  to  this  program  will 
include  a  deck  of  cards  containing  a  sufficient  number  of  raysets  to 
describe  a  family  of  calculated  raypaths  through  the  ionosphere  of 
interest.  Tlie  number  of  hops  can  be  varied  between  one  a>'d  three, 
although  it  would  be  easy  to  modify  this  procedure  to  accommodate  four 
or  more  hops. 

The  assumption  is  made  that  the  raypa'hs  represent  streamlines  in 
a  flux  of  energy,  and  that  the  streamlines  are  sufficiently  closely 
spaced  that  the  energy  can  be  considered  to  be  uniformly  spread  between 
them  in  space  and  time.  The  energy  tubes  defined  by  these  streamlines 
will  begin  at  the  transmitter  where  the  energy  density  in  each  tube  is 
determined  by  the  antenna  pattern.  Then,  as  the  tubes  travel  out  into 
space,  mix,  .spread,  and  cross  in  various  ways,  the  energy  within  each 
tube  is  conserved.  The  flux-tube  concept  permits  calculation  of  the 
density  distribution  of  the  energy  in  space  and  time  as  it  strikes  the 
ground  after  ionospheric  refraction  of  some  number  of  hops.  Notice  that 
two  or  more  flux  tubes  may  strike  the  same  area  on  the  ground. 

As  the  energy  travels  outward  within  these  flux  tubes,  corrections 
are  made  when  it  passes  the  D  region  to  account  for  absorption  at  the 
particular  range  and  angle  of  each  such  encounter.  Also,  corrections 
are  made  if  the  ravpaths  reflect  off  the  earth,  to  account  for  the 
incomplete  reflection  coefficient  at  the  particular  terrain  whore  the 
ray  landed.  For  example,  there  may  be  synthesized  a  mountain  range  of 
high  scatter  coefficient  but  low  reflection  coefficient,  which  has  a 
limited  extent  in  radial  range. 

When  the  tube  of  flux  arrives  at  the  ground  where  it  is  to  bo 
scattered,  the  scattering  coefficient  K(t)  is  used  to  determine  how 
much  of  the  energy  is  scattered  back  toward  the  transmitter,  taking  into 
account  the  angle  of  incidence  and  the  angle  of  scatter  necessary  to 
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permit  the  energy  to  return  to  the  transmitter.  Geometrical  relations 
applied  to  the  rayset  information  provide  full  allowance  for  the  effect 
of  earth  curvature  and  ionospheric  curvature,  both  along  and  transverse 
to  the  propagation  path. 

On  the  return  trip,  the  energy  is  again  attenuated  due  to  absorp¬ 
tion  in  the  D  region  and  incomplete  reflections  off  the  earth.  When 
it  reaches  the  receiver  location,  it  is  corrected  to  account  for  the 
antenna  gain  at  the  elevation  angle  of  arrival.  Also,  half  the  energy 
is  lost  because  of  polarization  rotation  in  the  ionosphere  which  causes 
it  to  bo  undetectable  by  the  recoiv'or  antenna. 

Tyiiically  there  will  bo  so\'oral  thousand  elements  of  energy  calculated 
in  this  manner.  Throughout  this  process,  the  time  of  flight  of  each 
element  is  accounted  for  and,  finally,  the  convolution  principle  is 
applied  to  coi-roct  for  the  actual  shape  of  the  individual  ti-ansmitted 
pulses  as  modified  by  the  receiver  bandpass  characteristics. 

The  preceding  discussion  is  a  rough  sketch  of  the  digital  computer 
process  which  is  about  to  bo  described.  During  the  actual  computer 
operation,  a  number  of  complicating  problems  arise,  such  as  the  presence 
of  raypaths  that  do  not  return  to  the  earth  because  of  a  tilted  ionosphere 
or  because  they  penetrate  and  escape  from  the  earth  completely.  Since 
.such  raypaths  may  form  one  of  the  boundaries  of  a  flux  tube,  it  is 
necessary  to  decide  how  to  process  the  semiboundod  energy  increments 
to  maximize  the  similarity  between  tlio  synthetic  backseat  tor  and  that 
wliich  would  be  obtained  experimentally  in  such  a  situation. 

.A  number  of  other  complications  arise  because  synthesis  is  done  in 
the  domain  of  discrete  rather  than  continuous  numbers.  The  particular 
computer  in  use,  an  IBM  7090,  provides  slightly  more  than  eight  decimal 
digits;  consequently,  round-off  orroi’  is  frequently  a  source  of  difficulty, 
particularly  when  it  describes  the  height  parameter  w’hich  is  mea.sured  from 
the  center  of  the  earth.  In  this  case,  the  round-off  error  is  approxi¬ 
mately  a  meter,  whicli  is  j'ust  large  enough  to  cause  occasional  trouble. 

The  authoi'  hopes  that  otJiers  may  wish  to  write  siinilar  computer 
programs  and  it  is  unlikely  that  the  Fortran  II  version  of  this  p.rogram 
would  bo  directly  useful.  Consequently,  this  chapter  will  bo  given  in 
sufficient  detail  to  enable  a  scientific  programmer  to  write  a  successful 
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backscatter  synthesis  program  in  another  computer  language  if  given  some 
assistance  by  a  person  familiar  with  the  physics  of  backscatter.  A  con¬ 
sistently  successful  program  can  be  written  in  from  two  to  four  weeks 
starting  from  this  description. 


A.  LIST  OF  SYMBOLS  AND  DEFINITIONS 

At  this  point  it  is  convenient  to  list  and  define  the  symbols  used 
in  this  study,  many  of  which  have  already  been  introduced. 


Symbol 

A 


a ,  b 


Description 

the  effective  receiving  area  of  an  antenna 
areas  on  the  surface  of  the  earth 

constants  used  to  calculate  gain 


B 


c 


bandwidth  of  a  receiver,  in  cycles  per  second; 
also,  in  raytracing,  the  number  of  times  a  ray 
hit  the  earth 

the  speed  of  light,  299.7925  km/msec 


D.E,F^,F2 


E,  E  ,  e.,  ,  e„,  e„ 
a  1  2  3 


E 


R 


conventional  designation  of  ionospheric  layers 

parameters  which  control  the  calculation  of 
absorption  according  to  the  formula, 

E  =  E^.(l  +  e^R  ^  e^R"  t-  e^) 

the  received  energy  in  a  single  modeset 


e 


2.7182818,  Napier's  constant 


F 


f 


frequency,  usually  given  in  kilocycles 
general  notation  for  a  function 


G 


antenna  gain  as  a  function  of  azimuth  and 
elevation 
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Symbol 


Description 


antenna  gain  separated  to  be  a  function  of 
elevation  alone 


g 


i 


H 


h 


I 


K 


a  constant  used  in  gain  calculations 

a  height  punched  into  a  rayset,  usually  the 
height  of  the  last  preceding  apogee 

height  in  general,  always  given  in  kilometers 

the  height  of  the  D  region 

an  indicator  number  in  raysets  giving  general 
ray  trajectory  shape 

a  function  of  angle,  terrain  type,  and  frequency 
which  determines  backscatter  coofficiont 

when  used  with  subscripts  of  presupersc r ipts , 
this  is  a  constant  used  to  calculate  K 


Boltzmann's  constant;  an  arbitrary  constant  used 
to  calculate  gain;  also,  the  absorption  coofficiont 

constants  usod  in  tho  descriptions  of  ionospheric 
tilts 


constants  usod  in  analysis  of  tho  gain  function, 

q 

a  constant  usod  to  describe  gain;  also  an  identi¬ 
fication  number  for  timosots 

specific  values  of  the  timoset  identification 

ionospheric  electron  density;  electrons  per  cubic 
cen  t imeter 
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Symbol 


Description 


N 

m 

P 


P 

n 


P 

e 


q 


R 


R.  ,R 


H 


1  2 
R,  R 


R 

P 


Hi.H^ 


R 


D 


maximum  electron,  density  in  an  ionosphere 
power,  measured  in  watts  or  milliwatts 
noise  power 

the  power  level  of  a  modeset 

received  power,  provided  transmitted  pulse  is  of 
"infinitesimal"  duration,  e 

gain  as  a  function  of  azimuth,  according  to 

the  relation,  G  =  gq 

range  in  kilometers.  However,  if  used  as  a  sub¬ 
script,  this  signifies  association  with  the 
receiver . 

the  lowest  and  highest  of  two  ranges  compared  by 
the  computer 

arbitrary  langes  used  to  establish  the  boundary 
between  differing  types  of  terrain  for  purposes 
of  calculating  K  and  (_ 

the  maximum  range  to  be  calculated 

ranges  calculated  by  the  computer  in  establishing 
the  identity  of  a  modeset 

the  center  of  the  range  interval  i-eached  by  a 
particulai'  energy  flux  tube 

for  a  particular  raypath,  the  ground  range  between 
encounters  with  the  earth  and  the  ionospheric  D 
layer 
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Symbol 


Description 


r/r 


r 


e 


T 


T  ■ 

P 


T  T 
^1’  2 


t/r 


t 


u 

V,  W 


X,  y 


z 


the  "received  rayset,"  that. is,  the  rayset  used 
to  describe  the  energy  path  on  its  return  from 
the  scattering  ground 

the  radius  of  the  earth,  6370  km 

time,  in  milliseconds.  However,  when  used  as  a 
subscript  it  signifies  association  with  the 
transmitter. 

group  and  phase  time  delays,  respectively,  as 
given  by  the  raysets 

the  maximum  time  used  in  calculating  backscatter 

the  time  limits  which  define  tlie  interval  during 
which  a  particular  flux  tube  encounters  the  earth 

the  time  limits  which  define  the  interval  spanned 
by  a  modeset 

the  "transmit  rayset,"  that  is,  the  rayset  used 
to  define  the  outgoing  flux  tube  which  strikes 
the  ground 

time,  as  a  variable  used  in  analysis 

constant  used  in  gain  analysis 

dimensions  used  in  a  reciprocity  calculation 

distance,  usually  in  Cartesian  coordinates 

distance  measured  along  a  raypath;  also  used  as  a 
parameter  in  the  formula  for  a  Chapman  ionospheric 
layer 


/• 
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Symbol 


Description 


a 


y 


y^.y^. 


.  .  y 
'  ll 


A 


e 


A 

5 


l.ll-  Ip.  I.p 


an  angle,  usually  measured  from  the  center  of  the 
ear  th 

the  ionospheric  absorption  index 

the  angle  of  a  raypath,  relative  to  the  ground, 
measured  at  the  radar  end  of  the  path 

the  average  beta  of  a  flux  tube 

a  constant  used  to  raise  or  lower  a  mathematically 
synthesized  antenna  beam 

an  angle  used  in  absorption  calculations 

various  intermediate  products  calculated  during 
backseat  tor  synthesis 

used  as  a  prefix  to  indicate  the  increment  of 
some  parameter 

in  continuous  numbers,  an  infinitesimal  time 
increment;  in  discrete  numbers,  a  short  time 
increment  on  the  order  of  0.1  msec 

an  angle  used  in  absorption  calculations 

the  geocentric  angle  subtended  by  the  end  points 
of  a  line;  for  example,  0  =  R  ;  r^ 

the  wavelength  of  the  radio  energy 

in  continuous  numbers,  a  dummy  variable 

intermediate  products  calculated  during  back- 
scatter  synthesis 
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Symbol 


P 


1  2 
P.  P. 


3 

P 


Pt-  Pr 


T 


<t> 


Description 

the  percentage  of  energy  specularly  reflected  off 
the  earth;  also,  the  radius  in  circular  or  polar 
coordinates 

the  reflection  coefficient  of  the  three  types  of 
terrain  which  are  separated  at  ranges  '^R  and 

the  product  of  all  the  reflection  coefficients 
encountered  by  a  multiple-hop  flux  tube  as  it 
travels  to  or  from  the  backscattering  ground 

the  duration  of  the  transmitted  pulse,  if  it  is 
rec  tangular 

azimuth  angle  in  polar  coordinates;  angle  in 
circular  coordinates 


4-  . 

W’  3db 


■if 


'  (prime) 


Ravset 


the  azimuthal  beamwidth  of  the  antenna,  measured 
between  nulls  or  3-db  points 

the  angle  of  a  raypath  relative  to  horizontal  at 
the  ground  end  of  the  path.  (The  analogous  angle 
at  the  radai-  end  is  p.) 

the  average  landing  (or  takeoff)  angle  of  an 
energy  flux  tube  as  it  arrives  at  (or  leaves) 
the  scattei-ing  ground 

The  primed  superscript  on  a  rayset  parameter 
means  that  it  is  selected  from  the  next  punched 
cai-d  in  a  deck  of  raysets. 

A  set  of  numbers  that  describe  a  particular  hop 
of  a  particular  ray.  It  can  also  contain  numbers 
that  characterize  the  energy  flux  tube  bounded 
by  the  ray  and  some  other  ray  with  a  higher  p. 


63 


SEL-65-002 


Symbol 


Description 


Modeset 


Tlmeset 


Each  propagation  mode  to  and  from  the  backseat terlng 

ground  leads  to  a  small  portion  of  the  total 

received  signal  of  power  level  diii-lng  the 

time  Interval  bounded  bv  T  and  T  .  The  three 

L  H 

parameters  thus  defined  are  termed  a  modeset. 

A  set  of  approximately  one  thousand  numbers  whlcli 
represent  the  received  power  distribution  as  a 
function  of  time  sampled  at  Intervals  of  e. 


B.  REARRANGEMENT  OF  THE  POTOR  EQUATION 

Chapter  III  led  to  a  long  expression  for  P^  which  may  be  rearranged 
Into  eleven  factors,  7^^  through  follows; 


cos 


V _ _  , _ J 


SP,,C 

T 


e 


5 


> 


6 


Here  a  prime  superscript  (e.g., 
from  the  neighboring  raysel  with 
a  symbol,  as  In  p^^,  means  that 
flux  tube;  that  Is, 


R')  means  that  the  parameter  Is  selected 
a  next  higher  value  of  p.  The  Itar  above 
It  Is  the  average  angle  from  a  particular 


P 


R 


P'  ) 
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Notice  that  the  original  raysets  described  single  raypaths  which  could 
not  carry  energy  because  they  had  no  cross-sectional  area.  However,  in 
these  calculations,  each  rayset  is  compared  to  another  which  pertains 
to  the  neighboring  ray  with  the  next  higher  p.  The  paired  rays  are 
considered  to  bound  a  flux  tube  so  that  p  is  the  elevation  angle  of 
the  flux-tube  centerline.  Then,  it  can  be  seen  that 


> 

>v  = 


% 


\  -  % 


For  purposes  of  summation,  the  ground  will  bo  considered  to  consist 
of  concentric  rings,  each  of  which  will  bo  used  in  sequence  as  a  source 
of  backscatter.  Each  flux  tube  defined  by  a  pair  of  rays  will  arrive 
at  the  backscattering  ground  to  strike  a  ring  which  includes  all  the 
area  between  range  R,^,  and  R,j,.  This  annular  section  of  ground  will 
be  denoted  by  the  expression  (R,j,,R,j,).  The  terrain  which  is  used  for  a 
single  summation  may,  or  may  not,  bo  identical  to  (R.j,,Rp  but  the 
computer  logic  will  be  such  that  the  ground  annulus  under  consideration 
always  lies  within  the  boundaries  of  (R.j,,Rj),  Suppose  the  ring  being 
calculated  is  (Rj^.Rg)-  Then,  is  the  increment  of  takeoff  angle 

which  describes  a  flux  tube  that  leads  only  to  (Rj^,R  ).  Consequently, 
^11  calculated  by  the  following  relation; 


Pf  - 

^11  R.,  -  R'  ^^2  ■  ’'l- 


Tu  <“2  ■  "P 


Lsing  this  notation,  it  is  seen  that  the  range  R  which  appears  in 
is  related  to  in  the  following  manner: 


r 


6 


R 


"1  "  "2 


Since  R^  and  R,^  are  known  prior  to  the  selection  of 

raypath,  the  synthesis  is  speeded  if  the  second  factor 

f erred  to  y^. 

6 


the  transmitted 
in  y^^  is  trans- 
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will  be  consLanL 


Except  for  frequency,  all  of  the  factors  in  y 

5 

in  every  backscatter  calculation.  The  power  is  chosen  as  1  kw 

and  the  infinitesimal**  pulse  direction  €  is  chosen  as  0.1  msec. 

Calculation,  in  an  appropriate  system  of  units,  then  shows  that  if  y 

equals  710/F  ,  the  power  will  come  out  in  milliwatts.  Rewritinc  P 

A 

incorporating  the  above  operations  yields 


y  ' 

10  n 


There  are  some  interesting  symmetries  in  the  foregoing  equation.  Notice 

that  y^  and  are  the  same  except  that  one  is  selected  from  the 

transmitted  rayset  and  one  is  selected  from  the  received  ravset-  y 

*  9 

^10  similarly  related.  Also,  y^  contains  one  factor  from 

each  of  the  two  raysets,  and  the  cosine  terms  in  y ^  and  y ^  follow 
the  pattern.  Thus,  each  of  these  can  be  calculated  as  soon  as  the  ray- 
sets  are  read  in  by  the  computei',  since  their  values  depend  on  no  other 
parameters.  Therefore,  they  will  be  combined  to  speed  the  process. 

Finally  the  power  can  be  written  in  an  expression  which  is  the  most 
convenient  form  for  digital  synthesis. 
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This  equation  looks  wrong  because  the  presence  of  Implies  that 

depends  on  which  of  the  flux  tubes  is  chosen  to  be  the  transmission 

path.  Time-reversal  invariance  tells  that  the  choice  should  not  matter; 

the  flow  of  energy  should  look  the  same  when  time  runs  backwards,  as  it 

will  appear  to  do  when  the  transmit  and  receive  raysets  are  interchanged. 

The  source  of  this  asymmetry  is  the  absorption  of  a  esc  term 

into  because  of  the  definition  of  backscatter  coefficient  which 

was  used  in  Ref.  33.  There  it  is  noted  on  page  43  that  K(\|f)  = 

sin  'If.j,,  where  cr^  is  the  standard  radar  cross  section  per 

unit  area.  This  seems  to  imply  that  cr  is  a  more  fundamental  measure 

o 

than  K(\ir). 

In  tlio  last  equation  for  P^,  absolute  values  are  taken  because  some 
of  the  factors  have  been  allow'ed  to  become  negative  in  the  expression. 
Notice  that  |  E  ,  7  and  7  can  be  calculated  as  soon  as  the 

1  j\  ^  *3 

raysets  are  read  into  the  computer.  Also,  7^  is  evaluated  only  once 

in  the  entire  program  since  it  depends  only  on  frequency.  When  the 

annulus  of  terrain  has  been  selected  for  calculation,  it  is  possible  to 

calculate  7^. 

6 

The  operation  which  is  carried  out  most  often  is  the  selection  of 
specific  transmit  and  receive  raypaths.  When  this  is  done,  it  is  only 
necessary  to  calculate  7  ,  as  all  the  other  factors  in  P.  are  already 
evaluated  and  stored  in  registei’s  of  the  computer.  The  three  parameters, 
P. ,  T  and  T  ,  then  constitute  one  modeset  which  might  be  thought 

ZA  11  Xj 


67 


SEL-65-002 


of  as  one  element  in  an  integration.  Because  of  the  streamlining  of 
this  operation  the  most  frequent  calculation  is  very  short.  Synthetic 
backscatter  calculations  in  which  there  are  6000  or  more  modesets  take 
less  than  60  sec  on  an  IBM  7090,  despite  the  formidable  appearance  of 
the  equations.  Thus,  careful  tailoring  leads  to  a  procedure  which 

is  very  fast  and  consequently  practical  for  synthesis  of  ground  back¬ 
scatter. 

C.  DETAILS  OF  THE  SY.VTHETIC  PROCESS 

When  backscatter  is  to  bo  synthesized,  the  digital  computer  is  given 
a  deck  of  cards  which  consists  of  three  groups;  (a)  the  program  deck, 
wi-itten  in  accordance  with  the  description  which  is  now  to  be  given; 

(b)  a  number  of  raysets  which  describe  propagation  through  a  specified 
ionosphere  at  a  specified  frequency;  (c)  a  set  of  propagation  parameters, 
which  consists  of  the  constants  which  the  operator  can  choose  at  will 
to  describe  the  ground  backscatter  and  refloctioti  characteristics,  the 
absorption,  antenna  gain,  pulse  length,  etc. 

A  good  deal  of  the  logic  inside  the  program  will  depend  on  the 
seciuence  in  which  the  raysets  are  road.  This  sequence  is  established 
prior  to  the  operation  of  the  computer  when  the  rayset  cards  are  shuffled 
in  an  automatic  sorting  machine.  All  of  the  cards  must  have  the  same  F 
and  the  same  IID.  Next,  the  cards  appear  in  major  groups  determined  in 
the  order  of  increasing  B;  that  is,  the  first  cards  ali  describe  the 
first  hop,  the  second  group  of  cai-ds  describe  the  second  hop,  and  so 
forth.  Within  those  groups,  the  cards  must  appear  in  order  of  increasing 
p.  In  operation,  the  program  will  read  the  first  rayset  card  to  obtain 
the  values  of  IID  and  F  to  bo  used  throughout  the  calculation.  The 
rayset  parameters  T^  and  H  are  never  used  in  this  particular  appli- 
ca  t ion . 

The  characteristics  of  the  terrain  are  described  to  the  computer  by 
the  following  set  of  14  constants; 
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Ill  this  list,  R  signifies  ground  range,  symbolizes  reflection 

coefficient,  and  K  symbolizes  ground-backscatter  coefficient.  As  the 
layout  of  those  constants  would  imply,  they  are  interpreted  according  to 
the  following  logic;  the  ground  range  radially  measured  from  the  radar 
is  divided  into  throe  sections.  At  all  ranges  loss  than  ^R,  the  con¬ 
stants  with  a  corresponding  prosuporscript  of  value  1  will  be  selected. 

1  2 

At  ground  ranges  botwoon  R  and  R,  the  reflection  and  backscatter 
constants  with  a  prosuperscript  of  2  will  be  selected.  Beyond  range 

n 

R,  the  constants  with  a  prosuporscript  of  3  will  bo  selected.  This 
system  could  easily  bo  extended  to  include  higher  numbers,  but  so  far 
added  versatility  has  not  been  desired. 

The  constants  are  applied  in  the  following  manner.  If  radio  energy 
strikes  the  earth  at  oblique  incidence,  some  of  it  reflects  as  it  might 
off  a  mirror,  and  some  of  it  is  scattered  in  all  directions.  The  amount 
reflected  is  then  ^  times  the  incident  power.  The  scattering  calcula¬ 
tion  is  more  complex  and  has  been  described  previously.  It  is  charac- 
toi’ized  by  the  tjackscatlor  coefficient  k(\1/),  which  is  calculated 
according  to  the  formula 

4  2 

K(tlr)  =  K  sin  tit  +  K  sin  '1/  +  K 

J.  ^  vi 

Generally,  the  user  of  this  program  will  have  some  maximum  time 
delay  beyond  which  results  are  not  desired.  In  order  to  take  advantage 
of  this  limit,  the  maximum  time  delay  is  given  to,  the  computer  so 

that  it  will  not  wastefully  calculate  power  which  would  appear  at  longer 
time  delays.  Associated  with  this  parameter  is  a  corresponding  maximum 


69 


SEL-65-002 


ground  range  which  is  also  given  to  the  computer.  This  is  symbolized 
R  ,  and  is  generally  calculated  as  the  round-trip  range  which  could  be 
reached  by  light  in  time  T^. 

In  this  particular  program,  the  decision  was  taken  to  make  e  =  0.1 

msec.  Any  function  of  time  is  sampled  at  increments  of  e  and  may  be 

described  either  as  a  function  or  as  a  sequence  of  numbers  which  represent 

sampled-data  points.  The  particular  value  chosen  for  e  could  easily 

be  changed,  but  0.1  msec  appears  to  bo  reasonable  for  two  reasons.  First, 

4 

the  bandwidth  of  a  backscatter  I’adar  seldom  exceeds  10  cps,  and  conse¬ 
quently  the  information  which  can  bo  received  is  fully  represented  with- 

-4 

out  any  approximation  by  samples  taken  at  increments  of  10  sec.  This 
result  follows  directly  from  the  sampling  theorem  of  information  theory 
[uof.  47].  Thus,  e  can  bo  taken  as  0.1  msec  and  the  resulting  record 
will  contain  as  much  or  more  detail  than  actual  experimental  records 
obtained  by  radars  of  bandwidth  10  kc  or  less. 

Another  reason  for  this  selection  is  that  the  final  result,  a 
drawing  of  power  as  a  function  of  time,  is  made  with  a  pen  in  an  auto¬ 
matic  plotting  machine.  With  the  largest  practical  scale  and  smallest 
pen  available,  the  line  width  is  about  0.1  msec. 

It  is  considered  that  the  maximum  transmitted  power  is  1  kw,  but  this 
is  not  very  significant  because  the  results  can  bo  multiplied  or  divided 
by  any  factor  to  represent  a  cori-esponding  higher  or  lower  transmitted 
power.  The  nature  of  the  transmitted  pulse  is  described  in  two  ways. 
First,  the  pulse  shape  is  considered  in  some  cases  to  be  square.  This 
approximation  permits  easier  running  of  llie  program  if  one  is  willing  to 
accept  the  corresponding  assumption  of  excessive  receiver  bandwidth,  as 
is  often  the  case  when  interest  is  centered  on  other  aspects  of  the 
backscatter  phenomenon.  The  square  pulse  is  fully  described  by  the 
selection  of  a  single  parameter,  pulse  duration,  which  is  measured  in 
milliseconds  and  symbolized  T. 

However,  if  the  pulse  shape  as  modified  by  the  receiver  bandpass 
characteristic  is  known,  it  is  sampled  at  intervals  e  and  the  resulting 
array  of  numbers  is  given  to  the  computer.  In  subsequent  operation  these 
data  points  are  considered  to  be  sample  points  of  the  function  f(T) 
previously  mentioned  in  the  description  of  the  effect  of  pulse  shape. 


SEL-65-002 


70 


The  necessary  convolution  is  carried  out  in  the  discrete  number 
system  where  it  appears  in  the  more  simple  form  of  a  weighted  summation. 

The  two  methods  by  which  the  gain  (as  a  function  of  p)  is  given  to 
the  computer  were  discussed  in  the  preceding  chapter.  When  using  the 
analytical  form  of  input,  the  computer  evaluates  the  expression  (pre¬ 
viously  given)  at  each  particular  value  of  p  where  it  needs  a  value  of 
gain.  Sometimes,  the  actual  gain  of  an  experimental  antenna  is  available 
in  graphical  or  tabular  form.  In  this  case  it  is  easier  to  use  the 
''itemized  gain"  approach  where  the  curve  showing  gain  plotted  against 
P  is  approximated  by  a  connected  group  of  straight-line  segments.  The 
connection  points  are  listed  as  a  series  of  pairs  of  values  of  gain  and 
p:  then  each  is  multiplied  by  g^  which  is  determined  by  planimeter 

integration.  These  sequences  of  paired  values  are  given  to  the  computer 
in  digital  form.  The  straight-line  segments  are  reconstructed  when  the 
computer  carries  out  a  linear  interpolation  to  find  a  value  of  gain  for 
some  value  of  p  between  two  of  the  given  pairs. 

The  amount  of  nondeviative  absorption  in  the  ionosphere  is  established 

by  the  selection  of  four  constants:  E  ,  e  e  and  e  .  These  are  used 

O  X  i(v  3 

in  the  cubic  equation, 


e(r)  =  E^|l  +  e^R  +  e^R"  + 

\diich  establishes  the  value  of  K  to  be  used  in  evaluating  the  absorp¬ 
tion  for  some  particular  pass  of  a  ray  through  the  D  region  at  range 
R.  Appropriate  corrections  are  made  for  the  angle  at  which  the  ray 
penetrates  the  D  region  and  also  for  the  effect  of  radio  frequency. 

D.  PROGRA.MMED  STEPS 

After  the  computer  reads  in  the  raysets,  it  changes  all  angles  to 
radians.  Then,  after  the  last  rayset,  it  makes  an  artificial  rayset  with 
the  follow'ing  values: 


71 


SEL-65-002 


B  =  previous  B  +  1 
T  =  T 

P 

R  =  R 

P 

P  =  previous  P  +  0.01 
>[r  =  previous  ilf  +  0.01 

This  procedure  facilitates  the  logic  which  will  follow,  and  serves  Lo 
shut  off  the  calculation  at  the  proper  point  without  special  tests. 

The  only  rays  which  reach  the  ground  will  have  I  =  0.  and  they  will 
be  the  only  ones  which  account  for  any  backscattered  enei g> .  To  simplify 
the  logic,  the  computer  takes  the  raysets  with  I  0  and  modifies  them 
as  follows.  In  the  first  rayset  of  each  hop,  it  replaces  T.  R,  and 
with  the  values  BT^.  BR^,  and  p.  In  all  other  cases,  it  replaces 
T  with  the  previous  T  plus  3  msec  and  it  replaces  R  with  the  pievlous 
R  plus  1000  km. 

Next,  the  computer  calculates  =  710/f^  This  is  merely  an  inter¬ 

mediate  product,  and  it  contains  the  implicit  assumption  that  peak  trans¬ 
mitted  power  is  1  kw.  Then,  the  computer  processes  the  ravsets,  adding 
and  deleting  values,  so  that  each  will  contain  information  which  describes 
the  flux  tube  bounded  by  itself  and  the  next  rayset. 

To  simplify  language,  the  following  description  will  make  use  of  the 
symbol  (•  )  which  means  "to  be  replaced  by."  Thus,  B  -  B  +  1  means 
the  same  as  the  command,  "Increase  B  by  unity. 

In  each  rayset,  compare  R  and  R' . 

If  ■R-R'l<1.0,  then  T'  *  T'  ( - )  and  R'  -  R  -  1.0 

This  step  is  necessary  because  the  array  of  R  represents  sample  points 
from  the  function  r(P)  and,  on  occasion,  two  successive  sample  points 
will  be  very  close  to  one  another.  The  result  would  imply  an  unreason¬ 
ably  high  energy  density  striking  the  ground  between  the  two  R's.  In 
some  cases,  further  ray tracing  would  resolve  the  problem  because  a  ray 
with  an  intermediate  value  of  P  might  fall  at  some  range  outside  of 
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the  interval  spanned  by  the  two  R's.  This  would  mean  that  the  flux- 
tube  approximation  broke  down,  and  that  the  energy  which  started  out 
between  the  two  rays  underwent  a  folding  operation  and  now  lay  partly 
outside  the  interval  spanned  by  the  rays. 

However,  raytracing  will  not  always  correct  this  situation  because 
the  corresponding  continuous  function,  R(P),  may  have  an  infinite 
derivative,  dp/dR.  This  is  a  very  common  phenomenon  and  occurs  at  the 
skip  distance  of  every  layer.  Ray  theory  in  the  continuous  number  sys¬ 
tem  indicates  that  energy  density  is  infinite  at  the  skip  distance. 

This  obviously  false  result  is  resolved  by  wave  theory  which  shows  that 
the  energy  distribution  at  the  skip  distance  is  somewhat  complicated 
but  not  infinite  [Ref.  48].  The  computer  operation  resolved  the  problem 
by  the  artificial  spreading  of  this  energy  over  a  1 -km  range  with  a 
corresponding  time  spread. 

If  R  lies  in  an  interval  whose  end  points  are  the  R's  of  two 
successive  raysets,  and  if  each  rayset  has  1=0  and  the  same  B,  then 
manufacture  a  new  rayset  with  R  •  ^R.  The  new  rayset  should  be  placed 
between  the  two  original  raysets  and  should  have  values  of  T,  p,  and 
wliicli  are  calculated  by  interpolating  between  the  values  in  the 
existing  raysets  in  a  manner  which  is  linear  with  respect  to  the  param¬ 
eter  range. 

2 

If  R  lies  in  the  range  interval  between  any  two  successive  raysets, 
make  new  raysets  again  according  to  the  criteria  described  above  for  ^R. 

This  operation  makes  an  artificial  i-ay  land  precisely  at  the  range 
which  separates  various  typos  of  terrain.  It  results  in  a  simplification 
in  the  logic  because  only  a  single  scatter  or  reflection  coefficient  will 
now  be  needed  in  the  calculation  of  a  single  modeset. 

In  eacli  rayset,  store  R^^  and  R^  which  arc,  respectively,  the 
lowest  and  the  highest  of  R  and  R'.  If  B  is  not  equal  to  B',  sot 
R^  •  Rp  +  1  anfl  "  Rp  +  2.  If  B  =  B'  and  1=0,  set  |  •  G(P). 
Otherwise,  set  |  =  0.  In  each  rayset,  store  p  =  (p  +  p ' ) / 2 ; 

\lf  =  ( +  \lf')/2;  erase  if  desired.  The  following  definitions  are 

made; 
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cos  p 


%  =  ■'c  (C  - 


R  = 


In  each  rayset,  multiply  ^  by  as  follows; 


If 

Q  =  1. 

go 

Lo 

Instruc  tion 

1 

If 

B  =  2, 

go 

lo 

Instruction 

2 

If 

B  =  3, 

go 

lo 

Instruc  tion 

3 

Instruction  1 


'  t'xp 


CS( 

I' 


CSC  r 


[e(Rj^)  +  E(r  -  Rp)] 


Go  to  Instruction  4. 


Instruction  2 


^1- 

^  exp 

CSC 

L  J 

[E(Rp)  +  e( 

'K 

2 

Elf  *  kJ  *  E(n 

If 

r/2  <  ^R, 

f;  • 

1 

p 

If 

’r  <  r/2 

i-  • 

2 

p 

If 

^R  <  r/2, 

• 

3 

p 

Go  to  Instruction  4. 
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Instruction  3 


*-  ^1^2 


[e(R„)  +  e(|  -  R^)  *  e(|  *  1^) 

*  =(f  -  »d)  *  '(f  *  "d)  *  -  “d)] 


If  R/3  ^ 


If  R^  <  R/3  R^, 


If  R^  <  r/3, 


Pi  P 


If  2r/3  s 

If  R^  <  2r/3  <  Rg, 

If  Rg  <  2r/3, 


Go  to  Instruction  4. 


In  each  rayset, 


Instruction  4 


2^ 


cos  P 


store  in  each  rayset  =  esc  \1/. 

Notice  that  each  rayset  has  a  value  for  each  of  the  following: 


B  p  P  \1/  T  R  R^  Rjj  e  >3 


Now  begins  the  creation  of  perhaps  several  thousand  modesets.  Examine 

all  the  raysets  and  find  the  lowest  R  . 

L 

Call  it  R^ ,  and  go  to  Instruction  5. 

Instruction  5 

Find  the  lowest  R  >  R^^  and  call  it  R^. 

If  R„  >  R  ,  go  to  Instruction  6.  Otherwise  calculate 

2  p 


r 


6 


/ 
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according  to 


Now  select  the  constants  for  use  in  calculating 
the  following  criteria; 


*^1  ^  1 
If  - "  R. 


1  "l  +  S  2 

If  R  <  - - -  ^  R. 


2  ""l  +  •'2 

If  R  <  r 


■^2' 

(\- 

‘^3) 

’^2- 

’^2- 

"'‘3) 

Select  for  processing  all  raysets  which  satisfy  both  of  the  following 
relations: 


0  ■;  R  -  R,  <  200  km 
■L  J_i 


•^2  ^ 


This  procedure  selects  for  processing  only  those  raysets  which  describe 

energy  that  propagates  to  or  from  the  range  interval  of  interest, 

(Ri.R^)'  Semibounded  flux  tubes  are  included  only  for  the  first  200  km 

beyond  R  .  (in  the  computer,  it  is  necessary  to  establish  a  maximum 
i_» 

number  of  such  raysets,  and  thus  far,  after  a  year  and  a  half  of  opera¬ 
tion,  the  number  has  never  exceeded  20.  Usually,  there  are  2  to  4.) 

If  no  such  raysets  are  found,  go  to  Instruction  5. 

In  each  of  the  selected  raysets,  store  and  calculated  as 

follows: 


R  -  R 

=  T  +  (t  -  T  ) 


R  -  R 


1 


1  *^2  ■ 

=  T  +  (t  -  T  )  - - 

R  -  R 
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It  automatically  follows  that  <  T^.  To  undoi-st and  this  opera¬ 

tion,  it  must  be  appreciated  that  we  are  now  processing  the  energy  that 
lands  within  the  range  interval,  The  raysets  will  commonly 

be  describing  energy  flux  tubes  that  land  on  the  ground  in  a  range 
interval  which  is  larger  than  and  Includes  (Rj^,R2).  It  is  then  neces¬ 
sary  to  calculate  the  time  of  flight  to  the  near  and  far  edges  of  the 
range  interval  under  consideration.  These  times  of  flight  are  T^^  and 
T^  under  the  assumption  that  the  energy  is  distributed  uniformly  in 
space  and  time  within  a  flux  tube. 

For  this  particular  range  increment,  the  processing  of  the  raysets 
is  now  complete.  Each  rayset  under  consideration  identifies  a  propaga¬ 
tion  mode  which  leads  to  the  strip  of  terrain  under  consideration,  and 
thus  each  such  mode  can  be  cithei*  a  route  fi'om  the  transmitter  to  the 
ground  or  a  return  route  from  the  ground  to  the  receiver.  If  there  is 
just  one  rayset,  the  situation  is  trivial  and  the  energy  goes  out  and 
back  along  the  same  path.  Suppose  that  there  are  two  raysets:  let  us 
call  them  A  and  B.  The  energy  can  go  out  on  A  and  return  on  A. 

It  can  go  out  on  A  and  return  on  B.  It  can  go  out  on  B  and  return 
on  A.  It  can  go  out  on  B  and  return  on  B.  Thus,  there  are  four 
propagation  modes  and  the  computer  must  calculate  four  modosets  for  this 

particular  range  interval.  If  there  were  three  raysets,  there  would  bo 
,  2  , 

(3  =  9)  modesets.  Similarly,  if  there  were  n  raysets,  there  would 

2 

be  n“  modesets  for  tliis  I'ange  increment.  One  of  the  novel  features  of 
tills  form  of  analysis  is  that  the  resulting  complicated  energy  distribu¬ 
tion  in  space  and  time  can  be  simulated  by  a  simple  accounting  procedure 
within  the  computer.  This  procedure  is  described  as  follows. 

From  among  the  selected  i-aysets,  pick  one  and  call  it  the  transmit 
ra\set  ( ablirei  iated  T^r).  Elements  selected  ft  an  this  rayset  will  be 
given  an  additional  T  subscript.  Select  another  rayset,  which  can  be 
the  same  as  the  T  K;  it  will  be  the  receive  rayset,  symbolized  R/R, 
and  its  elements  will  receive  an  R  subscript.  This  procedure  must  bo 
carried  out  for  each  of  the  possible  combinations  of  T'^R  and  r/r 
which  can  be  found  from  among  the  raysets  available. 
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1  2 

Define  T  =  T  +  T  and  T  =  T  +  T  .  Pick  the  lowest  of 
11  IK  2T  2R 

these,  round  off  to  the  nearest  multiple  of  e,  and  call  it  T  . 

L 

Similarly,  round  off  the  highest  of  the  numbers  and  call  it  T  .  If 

H 

the  roundoff  process  has  caused  them  to  be  equal,  set  T  =  T  +  e. 

2  1  H  L 

Also,  if  T  -  T  <  €,  set  T  =  T  +  €.  Now  calculate  t  according 

2  1  o 

to  the  following: 


(Recall  that  K  is  a  function.)  Now 

^  '^6  ^3R  ^R  ^  ^5^ 

This  corresponds  to  the  last  equation  in  Sec.  B,  thus  "closing  the 

loop"  with  the  analysis  previously  given. 

Now  create  1100  storage  registers,  each  identified  by  a  whole  number 

n  where  0  <_  n  <  1101.  These  registers  will  be  filled  with  the  numbers 

wliich  constitute  a  timeset.  Define  a  parameter  j  whicii  equals  the 

width  of  the  transmitted  pulse  measured  in  tenths  oi  milliseconds  (for 

6  =0.1  msec).  In  tiie  case  of  a  square  pulse,  j  =  lOr .  Then,  add 

to  each  element  of  the  timeset  in  which  lOT  +j«n<]0T  +j.  Now 

L  H 

go  back  and  select  a  now  t/r  and  a  new  r/r  from  among  the  ray sets 
which  roach  this  range  increment,  and  start  a  now  calculation.  If 

all  combinations  have  been  selected,  set  R^  =  R^.  Repeat  Instruction  5. 

Instruction  6 

The  timeset  registers  now  contain  P^(t). 

If  the  transmitted  pulse  is  square,  such  that  j  =  lOT,  find  the 
sum  of  the  first  j  timeset  elements  and  substitute  this  sum  for  the 
first  element.  Starting  at  the  second  element,  sum  the  first  J  ele¬ 
ments  and  substitute  the  sum  for  the  second  element.  Starting  at  the 
third  element,  sum  the  first  j  elements  and  substitute  the  sum  for 
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the  third  element.  Continue  this  operation  until  it  has  been  carried 

out  20T  times. 

P 

If  a  shaped  pulse  is  used,  substitute  the  following  procedure.  Pair 
up  the  first  j  timeset  elements  with  the  j  elements  in  the  trans¬ 
mitted  pulse  description,  in  reverse  order.  Multiply  each  pair  together 
and  take  the  sum.  Substitute  this  sum  for  the  value  in  the  first  timeset 
element.  Starting  at  the  second  element,  pair  up  the  first  j  elements 
with  the  j  factors  in  the  (reversed)  pulse  description.  Again,  sum 
the  product  of  the  pairs  and  this  time  put  the  result  in  the  second  time- 

set  element.  Repeat  this  process  20T  times.  The  timeset  register 

P 

now  contains  the  final  answer,  ready  for  plotting. 

Add  a  very  small  number  to  each  element.  This  number  should  be  small 
enough  to  be  negligible  and  yet  large  enough  to  permit  calculation  of  the 
logarithm  of  tlie  number  if  the  prior  entry  in  the  register  was  0.  Now 
find  the  each  clement  and  substitute  it  for  the  element. 

Plot  the  logs  on  a  linear  graph  as  a  function  of  n  for  all  values  of 
n  up  to  20T^.  The  result  is  then  a  semi logari thmic  representation  of 
tlie  received  power  as  a  function  of  time. 

Figure  21  shows  a  tviiical  result.  Although  100  db  of  amplitude  and 
100  msec  of  time  are  plotted,  the  useful  part  is  usually  about  30  db  X 
50  msec.  The  rest  of  tlie  plot,  particularly  the  bottom  portion,  servos 
mainly  to  provide  insight  in  case  of  computer,  plotter,  or  program  error. 
The  labels  showTi  were  also  automatically  drawn.  The  first  shape  shows 
f(t),  the  transmitted  pulse  as  modified  by  the  receiver  l^andpass  charac¬ 
teristic.  The  next  shape  plotted  is  the  resulting  backscatter  from  three 

-14 

hops.  The  backgiound  noise  level  would  be  near  10  mw. 


E.  ADDITION  OF  THE  THIRD  DIMENSION--FREQUENCY 

So  far  this  analysis  has  taken  place  only 
the  results  from  a  number  of  frequencies  will 
single  synthetic  "sweep-frequency  backscatter 
doing  this  is  a  straightforwai-d  simulation  of 
is  carried  out  manually  using  a  rather  simple 
technique  will  be  illustrated  by  an  example. 


at  a  single  frequency, 
be  combined  to  produce  a 
record."  The  method  of 
the  natural  process  and 
graphic  technique.  The 


Now, 
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P'lG.  21.  THE  FINAL  RESULT:  A  COMPUTER  PLOT  OF  THE  TIME  VARIATION  OF 
AVERAGE  BACKSCATTER  AMPLITUDE. 
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Figure  22  shows  a  group  of  f ixed-f requency  backscatter  records  obtained 
for  a  single-ionosphere  -under  a  set  of  fixed  circumstances.  The  only  thing 
that  changes  from  record  to  record  is  the  antenna  gain  pattern,  which  is 
varied  in  accordance  with  its  known  functional  relationship  to  frequency. 


DIFFERENT  FREQUENCIES. 
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One  other  parameter  varies  with  frequency  and  that  is  the  received 

noise,  which  is  known  to  change  in  the  upper  high-frequency  band  in 

pi-oportion  to  the  radio  wavelength  raised  to  a  power  which  is  near  2.3. 

In  a  previous  report,  the  author  consulted  a  number  of  sources  in  the 

literature  and  concluded  that  noise  power  can  realistically  be  simulated 

2  3 

by  the  expression  P  =  90kBA  ‘  [kef.  49].  In  this  expression,  P 

n  ri 

is  received  noise  power,  k  is  Boltzmann's  constant,  B  is  t  lie  receiver 

bandwidth,  and  A  is  the  radio  wavelength;  all  units  are  mks.  This 

expression  gives  the  average  noise  power  with  the  understanding  that 

actual  measurements  can  vary  from  this  value  by  quite  large  amounts.  In 

the  lower  half  of  the  3-  to  30-Mc  band,  where  man-made  noise  is  the 

primary  soui'ce,  the  day-to-day  and  band-lo-band  variation  can  be  orders 

of  magnitude  [kef.  50,.  Consequent! v,  the  expression  for  P  should  be 

n 

considered  as  an  indication  of  order  of  magnitude  and,  in  this  report, 
it  serves  the  very  useful  function  of  providing  a  reference  level  for 
power  measurements  which  behaves  in  a  manner  consistent  with  the  knowi 
variation  of  the  received  noise  levels  when  frequency  or  bandwidth  are 
changed . 

Since  it  will  be  desirable  to  assign  a  noise  power  level  to  each 
synthetic  backseat  ter  record,  the  receiver  bandwidth  must  bt*  specified. 

It  has  been  found  that  for  the  puipose  of  obtaining  sweep-frequency 
backscatter  records  experimentally,  the  optimum  bandwidth  is  very  near 
the  theoretical  minimum  value  of  1.2  divided  by  the  pulse  duration, 

T  [kef.  23].  Therefore,  unless  otherwise  stated  in  reference  to  a 
specific  syntlietic  record,  the  assumption  will  be  made  that  B  =  1.2/t. 
Using  this  relation,  the  noise  can  be  expressed  in  terms  of  frequency  and 
pulse  duration  only.  Wlien  evaluated  throughout  the  range  of  values  which 
are  useful  here,  the  received  noise  power  may  thus  bo  represented  by  the 
two-dimensional  graph  given  on  Fig.  23. 


Now  at  last  it  is  possible  to  synthesize  a  single  sweep-frequency 
backscatter  record.  Power  will  be  measured  in  decibels  above  the  noise, 
using  data  from  Fig.  23.  The  method  of  putting  together  several  fixed- 
fiequenct  synthetic  records  to  produce  a  single  sweep -f requency  record 
is  illustrated  by  the  three-dimensional  model  which  is  shown  on  Fig.  24. 
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FIG.  23.  THE  CALCULATED  NOISE  LEVEL  FOR  GIVEN  FREQUENCIES  AND 
PULSE  DURATIONS. 

Hero,  each  of  the  shaded  vertical  planes  represents  a  single  fixed- 
frequency  synthesis,  taken  from  Fig.  22.  The  gradually  descending  plane 
which  cuts  through  these  vertical  ones  then  represents  noise  power  taken 
from  Fig.  23.  The  gray  shades  painted  on  the  vertical  planes  represent 
bands  10  db  wide  measured  upward  from  the  noise  level. 

Of  course,  in  actual  practice  the  cardboard  model  is  not  marie,  but 
rather  a  tracing  process  is  used  to  produce  what  might  be  called  a 
■paint-by-number ”  pattern.  This  pattern  is  then  painted  in  shades  of 
gray  watercolor  to  produce  the  result  given  on  Fig.  24.  When  accompanied 
by  proper  frequency  and  time  scales  and  a  description  of  the  ionosphere, 
terrain,  and  radar  parameters,  such  a  painting  is  a  synthetic  sweep- 
frequency  backscatter  record. 

In  some  work,  however,  the  10-db  increment  of  power  level  is  too 
coarse,  so  that  o  her  representations  of  the  sweep-frequency  records 
are  occasionally  necessary.  For  example,  when  a  strip  of  mountainous 
terrain  is  simulated,  it  is  more  convenient  to  show  the  difference  that 
tliis  simulation  causes  in  the  sweep-frequency  record  as  compared  to  a 
similar  record  synthesized  without  the  mountain  range. 
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FIG.  24.  A  MODEL  SHOWING  THE  .ADDITION  OF  THE  FREQUENCY  DIMENSION. 


The  parlicular  backscatter  i-ecord  illustrated  on  Fig.  24  shows  little 
similarity  to  actual  backscatter  because,  among  other  things,  there  is 
no  absoi-ption.  Other  unrealistic  parameters  are  used  in  this  particular 
model  and  so  it  is  seen  that  the  third-hop  echo  is  almost  as  strong  as 
the  first-hop  echo,  whereas  the  natural  backscatter  dies  off  rapidly 
with  an  increase  in  the  hop  number.  Figure  24  was  generated  as  part  of 
a  study  in  which  the  effect  of  the  ionospheric  structure  was  isolated, 

F.  A  COMPARISON  OF  SYNTHETIC  AND  EXPERIMENTAL  RECORDS 

The  acid  test  of  the  complicated  technique  which  is  presented  here  is 
comparison  between  synthetic  and  experimental  recoi'ds  taken  under 
identical  conditions.  Previously  it  has  not  boon  possible  to  do  this 
under  well -controlled  circumstances  because  a  properly  instrumented 
backscatter  I’adar  did  not  exist.  Thus,  the  diiect  comparison  wa.s  pre¬ 
vented  by  limitations  in  the  experimental  technique  used  in  the  past. 

The  limitation  is  understandable  when  it  is  considered  that  it  has  not 
ooen  possible  before  to  make  use  of  all  the  information  which  is  now- 
found  to  be  necessary. 

In  order  to  permit  such  a  study,  a  wel 1 -cal i bra t ed  swoep-f requency 
backsca^  tor  radar  has  boon  assembled  by  the  Radioscience  Laboratory  for 
tlio  specific  purpose  of  permitting  the  comparison  to  be  made.  The 
radar  operates  into  an  antenna  whoso  pattern  has  boon  flown  and  the 
resulting  reduced  data  published  [Ref.  45].  The  use  of  this  antenna, 
kindly  made  available  by  the  Stanford  Itesoarch  Institute,  eliminates  one 
of  the  major  uncertainties  wliich  would  usually  limit  the  validity  of  a 
dii'oct  comparison  between  synthetic  and  expei-imental  records.  Also,  this 
I adar  has  good  calibration  in  power,  frequency,  and  group  time  delay. 
Vertical  ionograms,  also  provided  by  SRI,  are  obtained  to  maximize  the 
information  available  concerning  the  ionosphere.  The  effort  appears 
worthwhile  because  it  shows  a  promise  of  producing  information  about 
such  things  as  the  scattering  characteristics  of  the  ground  as  a  function 
of  frequency,  terrain  type,  and  incidence  angle.  An  example  of  such  a 
comparison  is  given  hei*e  to  show  that  the  computer  technique  is  capable 
of  producing  a  match. 
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Figure  25  shows  a  synthetic  record  and  a  corresponding  expei'imental 
record  which  was  obtained  by  a  Granger  Associates  sounder  located  at 
Stanford.  The  striking  resemblance  between  the  structures  is  readily 
apparent.  Notice  in  particular  that  the  shape  and  location  of  the  two- 
hop  echo  is  closely  simulated. 

Figure  26  shows  a  corresponding  experimental  and  synthetic  one-hop 
record.  For  synthesis,  it  was  necess.ary  to  use  horizontal  electron- 
density  gradients  in  order  to  simulate  the  sharply  upturned  tail  at  the 
high-frequency  end.  The  experimental  record  is  one  of  the  early  practice 
scans  taken  with  the  calibrated  radar  when  it  was  being  prepared.  Here 
again  it  can  be  seen  that  the  correspondence  between  svnthesis  and 
experiment  is  quite  close.  The  experimental  record  does  show  an  addi¬ 
tional  echo  on  the  near  side  of  the  main  trace;  this  is  probably  due  to 
an  antenna  side  lobe  which  looks  in  a  different  azimuth  at  a  slightly 
different  ionosphere;  however,  these  two  records  have  not  been  fully 
Investigated,  and  they  are  presented  here  only  to  show  that  there 
actually  is  some  correlation  between  synthetic  and  experimental  back- 
sca  1 1  er . 
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FIG.  25.  EXPERIMENTAL  AND  SYNTHETIC  RECORDS  OF  TWO-HOP  BACKSCATTER. 
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FREQUENCY  (Me) 

a.  Experimental  record 


FREQUENCY  (Me) 


b.  Synthetic  record 

FIG.  26.  EXPERIMENTAL  ANT)  SYNTHETIC  BACKSCATTER  WITH  A  TILTED 
IONOSPHERE. 
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VI .  ILLUSTRATION  OF  THE  SYNTHESIS  OF  A  SPECIFIC 
EXPERIMENTAL  RECORD 


There  are  two  primary  uses  for  the  digital  computer  technique 
described  in  the  preceding  chapters.  First,  it  can  be  used  to  investi¬ 
gate  the  effect  on  backscattor  of  changes  in  only  one  of  the  various 
factors  which  influence  the  record.  For  example,  it  is  possible  to 
change  only  the  antenna  pattern  and  hold  all  other  factors  unchanged 
to  determine  the  effect  whicli  the  antenna  has  on  backscatter.  Similarly, 
the  absorption  may  be  changed  in  some  predetermined  way,  while  all  other 
things  are  held  constant;  an  isolated  irregularity  could  be  inserted  in 
the  ionosphere  at  a  distance,  and  the  resulting  backscatter  perturbation 
could  clearly  be  illustrated  in  a  quantitative  manner. 

The  second  potential  use  for  the  digital  technique  is  in  the  syn¬ 
thesis  of  specific  experimental  records  which  ai-e  acquired  by  a  well- 
calibrated  backscattor  radar.  Prior  to  this  time,  there  has  been  to  the 
author's  knowledge  no  backscattor  radar  which  was  adequately  calibrated 
to  support  compai-ison  studies  of  this  typo.  The  lack  is  understandable, 
because  there  has  never  l)oon  any  prior  need  for  such  preci.so  quantita¬ 
tive  data.  However,  now  it  will  be  necessary  to  know  the  transmitted 
power,  tlic  pulse  shape,  the  effect  of  the  receiver  on  the  pulse  shape, 
th<'  antenna  pattern,  the  alisolute  power  level  of  the  received  signal, 
and  the  precise  fro(|Uonc\  scale  of  any  sweep-frequency  records,  to  men¬ 
tion  a  few  examples.  Also,  it  will  Ije  very  convenient  to  have  nearly 
simultaneous  "A"  and  records  of  the  same  Ijackscat  ter.  (An  A 

record  shows  backscatter  amplitude  as  a  function  of  group  time  delay  at 
a  fixed  frequency,  while  a  Z  record  shows  frequency  and  group  time 
delay  on  the  X  and  Y  axes  of  an  oscilloscope,  while  tjackscatter 
signal  amplitude  modulates  the  Z  intensity  axis.) 

One  of  the  major  proltlems  has  been  the  need  for  a  well-calibrated 
antenna  suitable  for  operation  in  the  frequency  interval  of  approximately 
7  to  30  Me.  Fortunately,  the  Stanford  Research  Institute  had  recently 
measured  the  radiation  pattern  of  a  rhombic  antenna  [Ref.  45],  and  they 
kindly  made  the  antenna  available  for  this  experiment. 
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In  mld-1964,  Iladioscience  Laboratory  personnel  began  the  assembly 
and  calibration  of  a  sweep - f requency  backscatter  radar  to  operate  into 
the  SHI  antenna.  It  has  proven  possible  to  build  a  radar  which  is 
satisfactory  for  these  purposes  by  using  equipment  already  on  hand. 

The  basic  item  was  a  near-prototype  version  of  a  sweep-f I’oquency  pulse 
radai-,  built  by  Applied  Technology,  Inc.  of  Palo  Alto,  Calif.,  and 
designated  the  ES-IA.  This  device  consists  of  a  transmitter  which  is 
sei'vo-tuned  so  tliat  it  tracks  the  frequency  tuning  of  a  modified  SP-600 
receiver.  Tlie  ES-IA  transmits  approximately  1  kw  in  pulses  of  variable 
duration  and,  in  the  most  common  situation,  this  energy  is  used  to  drive 
a  PA-IA  power  amplifier  made  by  tlie  same  manufacturer  and  similarly 
sei'vo  -  tuned .  Tlie  peak  transmitted  power  is  on  the  order  of  5  kw, 
although  the  exact  value  is  a  function  of  frequency  due  to  the  inter¬ 
action  of  a  large  number  of  circuit  elements  which  have  frequency- 
sensitive  impedances. 

The  first  calibrated  record  suitable  for  syntliesis  was  made  between 
1032  and  1042  Pacific  Daylight  Time  on  20  October  1964  ( PDT  is  GMT  minus 
7  hr).  This  record  is  shown  as  Fig.  27.  The  array  of  dots  is  placed  on 
the  record  b\-  a  system  which  provides  precise  quantitative  indicators 
of  the  frequency  and  the  group  time  delay  of  the  record.  It  can  be  seen 
that  there  was  some  Jitter  in  the  lime  marks  toward  the  high-frequency 
end  of  the  record,  but  stability  at  the  low  frequencies  pi’ovides  an 
unambiguous  indication  of  the  group  time  delay. 

Figure  28  shows  calibiated  A  records  which  were  acquired  shortly 
before  the  Z  record.  These  will  provide  the  quantitative  measure  of 
backscatter  signal  strength  for  use  in  the  synthesis.  However,  it  will 
be  seen  that  the  high-frequency  records  were  acquired  as  much  as  1  hr 
20  min  prior  to  the  Z  record,  and,  since  this  was  in  the  morning, 
adjustments  have  to  bo  made  to  allow  for  the  changes  in  the  ionosphere 
which  took  place  between  acquisition  of  the  A  and  Z  records. 

The  remainder  of  this  chapter  deals  with  the  use  of  the  digital- 
synthesis  technique  for  the  simulation  of  frequency-time  data  of  Fig. 

27  together  with  the  amplitude  information  from  Fig.  28. 


SEL-65-002 


90 


91 


SEL-65-002 


FIG.  27.  THE  FIRST  CALIBRATED  SWEEP-FREQUENCY  RECORD  FOR  SYNTHESIS. 


A.  DESCRIPTION  OF  THE  CALIBRATED  RADAR 


Figure  29  shows  the  equipment  which  has  been  assembled  to  gather 
calibrated  backscatter  records.  In  the  rack  nearest  the  camera  are  the 
exciter-sounder  and  the  power  amplifier  previously  mentioned.  In  the 
other  two  racks  are  the  modified  SP-600  receiver,  linear  detectors, 
logarithmic  detectors,  oscilloscopes  for  data  monitoring  and  recording. 
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a  polaioid  camera  for  A-scope  records,  and  a  strip-film  camera  for 
Z-scope  records. 

The  equipment  is  located  in  a  bus  which  is  temporarily  stationed  at 
one  end  of  the  rhombic  anteiina.  This  antenna  has  a  gain  pattern  which 
has  been  measured  by  aircraft.  The  azimuthal  gain  integration,  described 
in  Sec.  IVB,  can  be  done  quite  accurately  on  this  antenna  because  the 
patterns  are  precisely  known.  The  result  of  the  calculation  is  given 
on  Fig.  30  which  shows  the  scalar  mulliplioi-  of  the  backscatter  amplitude 
as  a  function  of  frequency. 


AZIMUTHAL  GAIN. 

The  variation  of  gain  with  respect  to  elevation  angle  involves  a  good 
deal  more  work.  Reforonce  45  shows  the  shape  of  a  number  of  gain  functions 
in  a  typical  polar  plot.  Those  are  transferred  by  a  simple  coordinate 
transformation  to  the  chart  which  has  been  devised  for  the  integration  of 
itemized  gain  (Fig.  19).  The  result  is  shown  on  Fig.  31  for  the  frequencies 
which  were  plotted  by  Barnes.  It  should  be  mentioned  that  those  curves, 
shown  on  Fig.  31,  are  the  theoretical  curves  given  in  Ref.  45,  and  that 
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Ihere  are  some  data  points  along  some  of  these  curves  which  differ 
measurably  in  a  few  cases.  When  these  data  are  used  for  the  synthesis 
of  backscatter,  it  will  be  shown  that  the  gain  patterns  lead  to  a  better 
match  between  synthetic  and  experimental  backscatter  records  when  they 
are  corrected  for  the  actual  gain  measurements  acquired  by  Barnes' 
aircraf  t . 
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It  is  necessary  to  be  able  to  interpolate  between  the  curves  given 
on  Fig.  31  in  order  to  arrive  at  a  gain  pattern  for  any  frequency 
desired.  Since  almost  nothing  involved  is  linear,  a  graphical  inter¬ 
polation  technique  is  used,  as  shown  on  Fig.  32.  If  one  visualizes  a 
three-dimensional  array  of  gain  functions,  with  the  absolute  value  of 
gain  being  displayed  normal  to  the  plane  of  the  paper,  then  the  various 
curves  on  this  chart  represent  the  intersections  of  the  gain  surface 
with  a  series  of  ton  planes  which  are  equally  spaced  in  amplitude  between 
zero  and  maximum  gain.  Notice  that  maximum  gain  is  normalized  to  be  1.0, 
since  tlie  gain  function  g(p)  can  have  any  amplitude  so  long  as  it  has 
the  correct  shape.  It  is  thus  possible  to  select  any  frequency  and  to 
road  from  this  chart  a  series  of  21  points  along  the  curve  of  gtp). 

Next  a  series  of  measurements  was  conducted  to  determine  the  shape 
of  the  transmitted  pulse  and  also  to  determine  the  shape  of  this  pulse 
after  it  is  received  and  detected.  To  this  end,  three  measurements  were 
made,  as  shown  on  Fig.  33.  The  top  part  shows  a  direct  measurement  of 
the  transmitted  pulse  wliich  was  obtained  simply  by  holding  an  oscilloscope 
probe  in  the  vicinity  of  the  antenna  lead  wires.  The  vertical  axis  is 
thus  uncalibraiod,  but  the  relative  proportions  are  correct.  It  can  be 
seen  that  the  transmitted  power  level  decreased  during  the  half¬ 
millisecond  pulse.  This  situation  will  bo  corrected,  but  it  was  occurring 
at  the  time  when  the  calibrated  record  of  20  October  was  made.  Careful 
measurements  show  that  the  transmitted  power  level  actually  decreased 
17  percent  from  the  leading  to  the  trailing  edge  of  the  pulse. 

The  second  picture  on  Fig.  33  shows  the  received  pulse  which  was 
acquired  by  removing  the  receiver  and  detector  about  5  miles  to  a  pai tially 
shielded  site  located  in  rolling  terrain.  It  is  seen  that  the  pulse  shape 
det ei'iora tes  even  further  after  it  is  received  by  this  method.  The  third 
picture  on  Fig.  33  shows  the  received  pulse  displayed  after  detection  by 
the  same  logari thmic -response  detector  which  was  used  for  the  measurement 
of  the  backscatter.  Careful  measurements  show  an  extremely  close  agree¬ 
ment  between  the  linear  and  logarithmic  records  of  the  received  pulse 
shown  here.  As  a  consequence  of  these  measurements,  the  pulse  shape  is 
represented  in  the  computer  by  a  series  of  numbers  between  0  and  1  which 
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FIG.  32.  C1L\RT  USED  FOR  DETERMINATION  OF  GAIN  AT  ALL  FREQUENCIES 


simulate  the  shape.  They  have  been  determined  to  be  the  following; 

0.177,  1.000,  0.762,  0.691,  0.623,  0.177,  0.036,  0.001.  When  synthetic 
backscatter  is  shown  later  in  this  chapter,  the  reader  will  be  able  to 
see  a  plot  of  this  synthetic  pulse  shape  and  compare  it  to  the  bottom 
picture  on  Fig.  33. 

The  receiver  passband  is  3  kc  wide.  Since  the  sweep  rate  is  approxi¬ 
mately  3  kc  per  100  msec,  it  can  be  seen  that  detuning  is  negligible  for 
the  under-30-msec  lime  delays  of  one-hop  Ijackscat t er .  The  A  records 
wore  generated  at  fixed  frequencies,  so  no  detuning  was  involved. 

The  exact  shape  of  the  frequency  scale  on  the  figures  is  a  function 
of  the  tuning  scale  on  an  SP-600  receiver.  The  swoop  is  accomplished 
by  turning  the  receiver  knob  at  a  constant  angular  rate  through  the 
7-15  and  15-30  Me  bands  of  the  receiver.  Within  those  bands,  the  fre- 
cgiency  variation  is  neither  linear  nor  logar  i  tlimic ,  but  is  actually 
somewhere  in  between.  .At  the  junction  of  the  bands,  a  series  of 
switchings  is  accomplished  manually  while  the  strip-film  camera  is  off, 
and  tlio  junction  is  not  generally  noticeable  on  the  records.  However, 
on  Fig.  27,  the  junction  can  clearly  be  seen  at  14.5  Me  because  an  auto¬ 
matic  camera  switch  had  not  yet  boon  installed;  also  the  gain  changed 
considerably  during  the  Ijand  switch.  Figure  34  shows  the  variation  of 
froquoncy  as  a  function  of  time  during  tlio  sweep.  Since  tlie  strip-film 
camera  moves  at  a  constant  rate,  the  horizontal  axis  may  l)o  interpreted 
as  a  direct  indication  of  the  horizontal  axis  on  the  liackscattor  records. 
Tlio  plotted  points  in  the  figure  show  the  variation  of  the  radar,  and 
th('  straiglit  lino  through  those  points  shows  what  the  ideal  track  would 
bo  if  the  frequency  sweep  wore  logarithmic.  It  can  bo  seen  that  this 
radar  is  nearly  logarithmic,  and  for  most  purposes  in  this  study  the 
straight  line  is  used  as  a  substitute  for  the  measured  values.  This 
loads  at  most  to  an  crrcir  of  perhaps  l/3  Me,  and  oven  this  small  error 
can  bo  taken  out  through  the  use  of  the  calibration  marks  on  the  films. 
However,  the  reader  will  see  that  great  advantages  follow  from  the  use 
of  a  logarithmic  scale  because  it  enables  the  use  of  transparent  overlays 
which  can  bo  slid  along  the  horizontal  axis  to  accomplish  what  is,  in 
effect,  a  multiplication  of  plasma  froquoncy.  This  will  beccimo  clear 
later. 
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FIG.  34.  SCALE  CO.MPRESSION  FACTOR  OF  THE  CALIBIUTED  ILADAR. 


B.  EXPERIMENTAL  CONTIITIOXS  FOR  THE  20  OCTOBER  Rl'N 

Figure  35  is  an  equidistant  azimuthal  projection  centered  on  San 
Francisco.  Con.sequen  t  ly ,  all  straight  linos  through  the  origin  are 
great  circles,  and  all  distances  measured  from  San  Francisco  arc  groat - 
circle  distances.  The  concentric  rings  drawn  around  San  Francisco  define 
ground  ranges  of  1000,  2000,  3000,  and  1000  km.  Since  backseat  ter 
time  delay  is  approximately  6.7  msec  per  1000  km.  it  can  be  soon  that 
tlio  backseat  tor  on  Fig.  27  came  from  ranges  not  exceeding  4000  km, 
approximately.  The  antenna  azimuth  is  122  deg  true,  as  shown  on  tlie 
figure.  Also  sliown  are  the  subsolar  points  at  various  times  of  day. 

These  are  computed  from  data  given  in  Ref.  51,  whore  the  declination  of 
the  sun  is  given  as  slightly  more  than  10  dog  south.  The  subsolar  point 
had  recently  crossed  the  west  coast  of  Peru  at  the  time  of  the  caliljratod 
radar  record.  Since  the  antenna  was  looking  almost  directly  toward  the 
sun,  the  electron  density  along  tlio  path  probably  increased  with  radar 
range.  SimilaiTy,  the  absorption  should  increase  with  range  away  from 
the  radar.  Tliroughout  the  path,  as  time  increases,  the  sun  approaches 
nearer  the  radar  and  consequently  the  electron  density  and  absorption 
will  bo  expected  to  increase  along  the  entire  length  of  the  path  during 
the  progress  of  this  particular  experiment.  The  A-scopo  records  which 
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IG.  35.  MAP  OF  THE  AMERICAS  SHOWING  RADAR  CONDITIONS  ON  20  OCTOBER. 
(Subsolar  points  reflect  various  times  of  the  morning.) 
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were  taken  earlier  than  the  Z-scope  record  should  reflect  fewer  electrons 
and  less  absorption. 

The  timing  of  the  various  events  is  given  on  Fig.  36.  As  shown,  two 
other  Z-scope  records  were  obtained  prior  to  the  one  sho\ni  on  Fig.  27. 
These  are  shown  on  Fig.  37;  both  were  obtained  prior  to  the  final  adjust¬ 
ments  on  the  radar,  and  consequently  they  have  a  flattened  shape  due  to 
a  different  camera  film  speed.  They  will  be  useful  hei’e  in  the  inter¬ 
pretation  of  the  A-scope  records  given  on  Fig.  28. 

Also  shown  on  Fig.  36  ai'o  the  times  at  which  a  large  number  of 
vertical-incidence  ionosonde  sweeps  were  obtained  by  a  r-2  sounder 
operated  by  the  Stanfoi-d  Research  Institute.  This  sounder  is  located 
near  the  calibrated  radai-,  and  consequently  its  measurements  can  be  used 
to  determine  the  state  of  the  ionosphere  over  the  radar.  A  sequence  of 
ten  sweeps  was  obtained  during  the  experiment.  These  traces  have  been 


THESE  Z  RECORDS  WERE 

THE  LAST  WITH  THE  THIS  Z  RECORD  IS 


PACIFIC  DAYLIGHT  TIME 

32689 


FIG.  36.  TIMING  OF  EVENTS  ON  20  OCTOBER. 
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copied  directly  from  the  experimental  films  by  the  autlior,  and  they  are 
presented  on  Fig.  38.  A  number  of  interesting  changes  can  be  seen,  but 
in  particular  notice  the  fortunate  circumstance  that  the  vertical- 
incidence  critical  frequency  remains  relatively  stable  over  an  hour- 
and-a-half  period  at  a  value  near  5  Me.  This  stability  proves  useful 
in  the  synthesis  of  the  ionosphere  as  it  varies  along  the  path.  Lower 
layers  were  present,  and  it  appears  that  they  increased  in  density 
during  Llie  advance  of  the  morning  hours.  Three  of  the  most  interesting 
freciuencies  were  transferred  to  Fig.  36  whicli  shows  the  time  variation 
in  the  vertical-incidence  ciitical  frequency,  the  cusp  in  the  upper  trace, 
and  llie  E  layer. 

The  staff  of  the  Stanford  Research  Institute  has  reduced  by  computer 
Lechnic|ues,  the  ionogi'am  for  0850  PDT,  anci  the  results  of  Llieir  calcula- 
li^ons  are  given  on  Fig.  30.  'fhis  I'oduction  indicates  the  presence  of  a 
main  layer  of  critic'al  frec|uency  of  5  Me  at  an  altitude  neai  2.>0  km. 

Al.so,  there  appear  to  l>o  hints  of  layei’s  at  120  and  160  km  with  critical 
frequencies  of  approximately  2.8  and  3.4  Me.  The  ionogram  itself  is 
sliown  in  Fig.  40,  and  it  can  be  seen  that  the  traces  are  not  sufficiently 
complete  to  tte  unamijiguously  interpreted.  Consequently,  the  layer 
structure  given  in  Fig.  30  can  only  l)o  considered  as  a  good  estimate  of 
fne  state  of  the  upper  atmosphere. 

.Since  the  experimental  recorti  w.as  oljtained  at  a  time  when  tiie  radar 
ant enra  was  directed  towaifl  the  subsolar  point,  there  was  a  considerable 
range  dependence  of  the  sun's  zenith  angle,  )(.  In  Sec.  IIC  it  was 
explained  how  absorption  i.s  accounted  for  with  a  function  E(r)  whoso 
magnitude  and  shape  are  dotorminecl  by  the  selection  of  the  four  parameters 
in  the  ociuation 

E(r)  =  ''3'^^) 

with  power  loss  represented  l)y  the  factors  of  the  form 
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FIG.  38.  VERTICAI.- INCIDENCE  TRACES 
DURING  THE  MORNING  OF  20  OCTOBER. 
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FIG.  40.  PHOTOGRAPH  OF  THE  TRACE  USED  IN  FIG.  39. 


The  mnsni  tilde  of  tlie  numbers  whicli  can  bo  inserted  in  these  equations 
can  bo  obtained  from  a  comprehensive  study  of  skywave  field  intensities 
published  by  Laitincn  and  Haydon  in  Ref.  52.  There,  analysis  of  a 
^reat  deal  of  experimental  data  led  to  the  conclusion  tliat  absorption 
could  Ito  simulated  by  a  function  of  two  variables;  the  solar  zenith 
aiiftle  and  the  /.uiich  sunspot  number  S.  In  Ref.  53  it  was  pi’odicted 

lliat  Iht'  value  of  S  at  llie  time  of  llie  20  0ctobf.‘r  e.xperimenl  should 
have  been  appro.ximatcly  17.  The  accuracy  of  sucli  a  prediction  is  cer¬ 
tainly  adequate  for  these  purposes,  because  tlie  expression  for  E  will 
be  seen  to  depend  only  sliRlitly  on  the  exact  value  of  S  during  those 
limes  of  weak  solar  activity. 

From  Fip;.  88  of  Ref.  52  it  is  found  that  for  S  =  0,  ^  =  0,  and 

F  7  Me,  tlie  \ erticul -incidence  absorption  is  approximately  10.5  db. 

At  vertical  incidence,  esc  .  =  1.0,  and  calculation  yields  a  value 

of  E  5,93  ■  10^.  This  establishes  the  maximum  value  of  E  anywhere 
on  the  surface  of  the  earth  for  S  =  0.  Then,  use  of  Eq .  15  from  Ref. 

52  gives  a  value  of  E  as  follows: 
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E(X.S)  =  (5.93  X  10^)  (cos^'^  0.881X)  (l  +  0.0037s) 

Inserting  S  =  17, 

E(x)  =  6.3  X  10^  ros^'^  (0.881x) 

The  final  step  is  the  expression  of  this  function  in  terms  of  the 
variable  R  through  use  of  the  cubic  equation  which  is  built  into  the 
program.  This  task  is  carried  out  on  Fig.  41  which  shows  two  families 
of  curves,  one  of  which  gives  the  zenith  angle  as  a  function  of  range. 

The  families  of  curves  are  unfortunately  necessary  because  of  the  long 
time  interval  during  which  the  records  were  gathered  on  20  October.  In 
the  ideal  situation,  the  experimental  A  and  Z  records  might  all  be 
obtaitied  within  a  time  span  of  perhaps  15  min,  and  in  such  a  case  it 
would  be  necessary  to  calculate  only  one  function  E(r). 

The  zenith-angle  curves  on  Fig.  41  wore  calculated  simply  by  laying 
out  the  radar  azimuth  atid  the  subsolar  point  on  a  globe.  Measurements 
wore  then  directly  made  of  the  geocentric  angle  between  the  subsolar 
point  and  points  along  the  antenna  azimuth  at  the  desired  ranges.  This 
is  a  quick  procedure,  and  yet  it  yields  data  of  adequate  accuracy.  The 
values  of  X  (-  geocentric  angle)  thus  derived  were  inserted  in  the 
equation  for  E(x),  yielding  points  along  tlie  curve  that  defines  E(r). 
The  final  stop  was  the  selection  of  the  parameter.s  in  the  culiic  equation. 
Tliis  is  done  as  follows:  the  value  of  E^  is  simply  the  value  of  E(o) 
and  thus  it  can  be  directly  road  off  the  graph.  Next,  a  straightedge  is 
laid  over  each  curve  of  E(r)  so  that  it  passes  through  E(o)  and  lies 
along  the  approximate  place  whore  the  user  desires  to  have  the  straight- 
lino  function,  E  (l  +  e  r).  Using  this  technique,  it  is  assumed  at  the 

’  a  1 

outset  that  o^  will  bo  zero,  so  that  only  the  second-order  term  will 

be  used  as  a  correction  to  the  straight  line.  It  is  then  possible  to 

2 

lay  the  straightedge  in  such  a  position  that  some  multiple  of  R  will 
provide  a  good  final  correction.  The  selected  multiplier  is  then  o^, 
and  the  value  of  e^  is  found  from  the  slope  of  the  straight  line. 
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FIG.  41.  GRAPIIICAI,  UEHIVATION  OF  THE  CUBir  EXPRESSIO.V 
FOR  ABSORPTION, 


Using  this  technique,  the  constants  were  selected,  and  they  are  given  in 
Table  2. 

In  order  to  show  the  quality  of  the  simulation.  Fig.  42  shows  E(r) 
calculated  from  the  values  given  in  Table  2,  and  it  should  bo  compared 
to  Fig.  41.  In  this  regard,  note  that  very  little  of  the  energy  involved 
in  the  calculation  of  backseat  ter  will  undergo  absorption  at  ranges  in 
excess  of  3000  km,  so  a  slight  mismatch  at  the  upper  ranges  is  tolerable. 
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TABLE  2.  CONSTANTS  SELECTED  FOR  SIMULATION 
OF  ABSORPTION  VARIATION  VS  RANGE  AND 
TIME  OF  DAY 


Time  of 
Morning 

(pdt) 

E 

a 

— 

‘^1 

— 

^2 

0900 

2.10  ■  lo' 

-4 

4.07  \  10 

-1.52  •  lo'^ 

0930 

2.55  lo' 

3.26  ■  lo'"^ 

-1.25  ■  lo'^ 

1000 

2.95  10^ 

-4 

2.63  10 

X 

o 

X 

o 

1 

1030 

3.32  X  lo' 

2.23  \  lo’"* 

X 

1 

o 

/ 

o 

o 

1 

CURVES  OF  FIG.  41. 
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C.  LOGARITHMIC  SLIDERS  SHOWING  BACKSCATTER  SHAPE  FOR  TYPICAL  ION 

LAYERS 

A  major  factor  in  the  synthesis  of  backscatter  is  the  generation  of 
a  correct  ionosphere  so  that  the  sweep-frequency  backscatter  echo  will 
have  the  correct  shape.  The  electron  density  will  vary  as  a  function 
of  altitude,  radial  range,  transverse  distance,  and  time  of  day.  It  will 
probably  never  be  possible  to  simulate  the  structure  of  an  entire  iono¬ 
sphere  over  a  backscatter  path,  although  some  success  has  been  attained 
in  measuring  the  variation  of  the  ionosphere  as  a  function  of  height  and 
radial  range  through  the  use  of  a  series  of  spaced  sounders  along  a 
great  circle  ’^Ref.  54]. 

One  of  the  best  ways  for  determination  of  the  variations  in  the  iono¬ 
sphere  is  through  the  use  of  sweep-frequency  backscatter  records  such  as 
the  one  being  synthesized.  The  shape  and  appearance  of  such  backscatter 
can  yield  a  great  deal  of  information  about  the  ionosphere  if  it  can  be 
properly  interpreted.  It  does  appear  that  the  use  of  a  digi tal -synthesis 
technique  will  lead  to  a  better  understanding  of  the  ways  in  which  iono¬ 
spheric  variations  affect  backscatter  so  that  eventually  records  them¬ 
selves  can  he  used  for  a  quantitative  study  of  upper  atmospheric 
cond i lions . 

In  order  to  simulate  the  20  October  record,  it  was  necessary  to 
generate  a  large  number  of  synthetic  backscatter  records  with  various 
ionospheres  in  order  to  determine  the  nature  of  the  horizontal  gr.'idients 
and  the  vertical  layer  structures  which  were  present.  The  first  require¬ 
ment  was  a  quantitative  guide  to  the  shape  of  the  leading  edge  of  back¬ 
scatter  which  results  when  the  ionosphere  is  composed  of  a  single  Chapman 
layer  without  horizontal  gradients.  In  order  to  determine  this  shape, 

24  Chapman  layers  with  various  heights  and  thicknesses  were  synthesized, 
and  rays  were  calculated  through  them  at  a  large  number  of  f reqtiencies . 

A  subsequent  computer  program  operated  on  the  raysets  to  determine  the 
frequency  variation  of  the  minimum  group  time  delay  for  any  given  iono¬ 
sphere.  Each  set  of  such  data  yields  a  curve  which  is  one  of  the  most 
important  characteristic  parameters  for  the  particular  ionosphere.  The 
24  resulting  curves  are  presented  on  Fig.  43  which  is  a  reduced  photograph 
of  transparent  plastic  ovei'lays  generated  from  the  data. 
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FIG.  43.  LOGARITHMIC  SLIDERS 
CHAPMAN  IONOSPHERES . 
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The  three  overlays  on  Fig.  43,  which  are  termed  "sliders,"  contain 
a  sufficient  amount  of  information  to  completely  describe  the  leading 
edge  of  backscatter  which  will  result  from  Chapman  layers  of  any  realis¬ 
tic  height,  thickness,  or  critical  frequency.  Because  of  this,  the 
sliders  represent  a  very  powerful  tool  for  the  interpretation  of  sweep- 
frequency  backscatter  records.  They  can  be  placed  over  the  experimental 
records  and  slid  parallel  to  the  frequency  axis  in  an  attempt  to  find  a 
match  between  tlie  theoretical  and  experimental  backscatter  edges.  If 
a  vertical-incidence  trace  is  available,  as  it  is  here,  then  the  critical 
frequency  f^  is  a  known  parameter,  and  the  slider  can  be  placed  in  a 
specific  location.  This  general  procedui’o  is  possible  only  because  the 
frequency  scale  is  logarithmic,  and  it  should  bo  emphasized  that  such  a 
convenient  technique  could  not  bo  developed  without  a  log  scale. 

For  calculation.  Chapman  layers  were  formulated  as  follows: 

-3 

N  =  electron  density,  cm 

N  =  maximum  electron  density  ^  10^  always,  altliough  it  does  not 
m 

matter  (f  -  9  Me ) 
c 

H  =  height  at  which  N  occurs 
m  m 

H  =  scale  height 
s 

N  =  N  exp  (l  -  z  -  e  whore  z  =  (height  -  H  )/lI 

m  'ms 

Note  that  tliere  are  three  variables:  N  ,  II  ,  and  H  .  On  each  slider, 

mm  s 

only  one  of  these  can  be  varied  to  generate  a  family  of  curves.  On  Fig. 

43  it  was  decided  to  var^■  the  parameter  II  because  it  has  a  very  strong 

m 

effect  on  the  shape  of  the  curves.  For  different  sliders,  different 

values  can  be  used  for  anothei'  parameter.  It  was  decided  to  make  one 

slider  for  each  value  of  H  because  this  has  only  a  weak  effect  on 

s 

backscatter,  as  examination  of  the  curves  will  show.  The  third  variable, 

N  ,  is  taken  care  of  without  further  calculation  by  virtue  of  the  log 
m 

scale.  Notice  that  the  index  of  refraction  is  calculated  from  the 
relat ion 


n  = 
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so  that  all  work  is  parametric  in  the  ratio  F  /f,  not  in  F.  Also, 

P 

F  is  related  to  N  by  the  relation 
P 

F  =  v/so.  6N 
P 

The  fortunate  result  is  that  any  calculation  for  a  Chapman  layer  with  a 

particular  N  will  yield  results  which  are  directly  applicable  for 
in 

other  yalues  of  N  provided  the  operating  frequencies  are  adjusted  to 
m 

maintain  a  constant  value  of  the  parametei'  slidei's,  f 

is  used  to  denote  frequency.) 

Usually,  the  ionosphere  has  horizontal  gradients.  In  particular 
during  the  20  October  run  when  the  antenna  was  directed  toward  the  sub¬ 
solar  point,  the  electron  density  must  have  increased  witli  increasing 
range.  As  an  initial  study  to  determine  the  effect  of  horizontal 
gradients  on  the  shape  of  the  slider  curves,  four  tilted  ionospheres 
were  genei'aied  and  their  corresponding  minimum- 1  Ime-delay  curves  have 
been  calculated.  The  results  are  presented  on  Fig.  4A,  togetlier  witli 
a  curve  from  one  of  the  nontilted  ionospheres  whicli  is  most  similar  to 
the  lilted  examples  given.  For  this  set  of  figures,  the  layer  scale 
height,  a  measure  of  thickness,  was  fixed  at  100  km.  Throughout  the 
I'ange  of  heights,  the  electron  density  was  programmed  to  tie  a  linear 
function  of  ground  range,  either  increasing  or  decreasing.  In  the 
stronge.st  tilts,  the  vortical -incidence  critical  frequency  was  allowed 
lo  floublo  in  a  distance  of  4000  km;  in  the  weakest  tills,  the  critical 
frequency  tlouliled  in  8000  km.  Hays  were  traced  in  both  directions 
through  tliese  tilted  ionospheres,  so  it  is  possiitle  to  see  also  the 
effect  when  the  critical  frequency  decreases  by  a  factor  of  2  in  these 

distances.  In  each  case,  f  is  measured  over  the  radar. 

c 

Ideally,  as  a  method  of  searching  for  the  shape  of  a  given  ionosphere, 
one  should  have  a  complete  set  of  curves  such  as  is  given  on  Fig.  44  for 
each  of  the  24  individual  situations  depicted  on  Fig.  43.  This  would 
require  a  total  of  24  sliders  so  that  the  user  could  interpolate  between 
the  various  situations.  Although  the  work  has  not  yet  boon  done,  it  is 
possible  to  use  the  tilted  slider  curves  of  Fig.  44  as  a  guide  for  the 
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interpretation  of  the  nontilted  curves  so  that  one  can  guess  fairly 

closely  the  nature  of  the  ionosphere  from  an  observation  of  the  shape  of 

the  experimental  minimum-T  curve.  This  procedure  was  used  for  the 

g 

synthesis  which  follows,  and  it  will  be  illustrated  by  a  sequence  of 
studies  wherein  the  author  has  estimated  the  shape  of  the  ionosphere, 
pypthesized  and  rayi raced  through  it,  and  then  produced  the  corresponding 
slider  to  sec  if  further  corrections  are  needed  to  match  the  shape  of 
the  20  October  experimental  curve.  It  will  be  seen  that  an  extremely 
close  match  is  obtained  in  only  two  iterations  by  this  method. 

D.  .MATCHI.N'G  THE  LEADING  EDGE  OF  THE  20  OCTOBER  BACKSCATTER 

The  first  step  in  matching  an  ionosphere  to  the  experimental  record 

consisted  of  a  comparison  between  the  shape  of  the  leading  edge  of  the 

backscatter  and  the  shape  of  the  various  curves  on  the  transparent 

sliders  when  f  was  set  on  the  5-Mc  mark.  As  an  initial  step,  it  was 
c 

decided  to  try  to  match  the  shape  of  the  backscatter  by  a  linear  varia¬ 
tion  of  electron  density  vs  range  at  every  height  within  the  layer. 

Study  of  the  curves  indicated  that  the  ionosphere  contained  a  layer 
which  was  somewhat  like  one  of  those  described  in  Table  3. 


TABLE  3.  FIRST  ESTIMATE  OF  ION-LAYER 
PARAMETERS  FOR  20  OCTOBER 


Ionosphere 

I  den  tification 
Number  (IID) 

m 

«m 

Hs 

Critical- 

Frequency 

Doubles 

(km) 

789 

3.1  X  10^ 

250 

83.3 

6370 

790 

3.1  X  10^ 

230 

76.7 

6370 

791 

3.1  X  10^ 

270 

90.0 

6370 

792 

3.1  X  10^ 

250 

83.3 

8000 

1  793 

3.1  X  10^ 

230 

76.7 

8000 
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In  this  table  the  first  column  contains  the  three-digit  IID  which 

is  a  convenience  and  serves  as  a  name  for  the  layer  during  its  study. 

Notice  that  all  these  layers  have  an  N  of  3.1  x  10^,  so  that  all 

m 

f^  =  5  Me.  This  group  of  ionospheres  was  selected  to  permit  a  comparison 

in  the  following  manner.  Ionospheres  789,  790,  and  791  all  have  the 

common  property  that  H  =  1 ^3  H  .  Thus,  each  ionosphere  would  be 

s  m 

classed  as  "thick"  in  the  terminology  of  Fig.  43.  In  other  words,  the 
first  three  ionospheres  in  Table  3  are  thick  ionospheres  with  the  same 
amount  of  tilt,  and  the  only  thing  that  is  different  is  the  layer 
height.  The  results  of  ray  calculations  through  these  ienospheres  are 
shown  in  Fig.  45a.  The  effect  of  the  layer  heiglu  can  clearly  be  seen 
and,  when  the  curves  are  compared  to  the  experimental  data,  it  is 
evident  that  the  horizontal  gradient  is  too  strong. 

Also  shown  on  Fig.  45a  is  another  curve  which  might  be  thought  of 
as  the  trailing  edge  of  the  backscatter;  that  is,  it  is  a  group  time 
delay  of  the  ray  which  takes  off  at  zero  degrees  tangent  to  the  earth. 

If  one  considers  only  the  lower  rays,  then  this  curve  represents  the 
longest  time  delay  that  can  he  generated  with  the  particular  ionosphere. 
However,  it  will  be  shown  later  that  in  the  higher  frequency  portion  of 
the  curves,  the  upper  (Pederson)  rays  play  a  major  role  in  the  structure 
of  backseat  ter . 

In  Fig.  45b  are  shown  the  ionospheres  790  and  793.  Both  are  thick 

Chapman  layers  at  an  altitude  of  230  km,  and  they  differ  only  in  the 

» 

severity  of  the  horizontal  gradient  present.  Figure  45c  shows  a  similar 
compari.son  between  thick  layers  wliich  are  250  km  high.  ttTien  these  over¬ 
lays  were  compared  to  the  experimental  data,  it  could  be  seen  that  the 
later  height  was  slightly  in  excess  of  250  km,  and  that  the  horizontal 
gradient  must  have  had  a  value  somewhere  between  the  two  values  which 
were  computed.  Ionosphere  792  is  quite  close  to  the  right  answer,  and 
conseciuent Iv  another  iteration  was  tried  with  smaller  variations  in  the 
parameters.  The  new  ionospheres  selected  for  trial  are  described  in 
Table  4. 

The  results  of  this  calculation  are  given  in  Fig.  46.  In  part  (a) 
are  shown  the  comparisons  between  the  results  for  III)  797,  798,  and 
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a.  Thick  layers  at  various  heights 


b.  Layers  230  km  high 


c.  Layers  250  km  high 


FIG.  45.  BACKSCATTER  OUTLINE  FOR  TILTED 
IONOSPHERES  OF  VARYING  HEIGHTS  AND 
GRADIENTS. 


SEL-65-002 


118 


TABLE  4.  SECOND  ESTIMATE  OF  LAYER  PARAMETERS 


IID 

N 

m 

H 

m 

H 

s 

Distance  to 

Double  f 

c 

795 

3.1  X  10^ 

260 

60.0 

7650 

796 

3.1  X  10^ 

260 

40.0 

7650 

797 

3.1  X  10^ 

260 

86.7 

9560 

798 

3.1  X  10^ 

260 

86.7 

7650 

799 

3.1  X  10^ 
_ 

260 

86.7 

6370 

25[- THIN  LAYER:  Hs--40km 

MEDIUM  LAYER:  Hs=60km  v 

6  20F 


I 

cr 
o  0 


|•spTHICK  layer  06.7km 

^  .0 


110*796 

7% 

7qH 


10 


*  -  -i  -i-  .  -L  i  .  X 

12  14  16  18  20  22  24 

FREQUENCY  (Mcl 


b.  Layers  of  various  thicknesses 

( II  varies) 

^  s  ' 

FIG.  46.  BACKSCATTER  OUTLINE  FOR 
VARYING  THICKNESSES  AND  GR.ADIENTS. 
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799  which  are  all  thick  Chapman  layers.  Comparison  of  these  three  curves 
with  the  experimental  record  shows  that  the  center  one,  No.  798,  has  the 
same  shape  as  the  experimental  record.  However,  the  trailing  edge  of 
this  echo  does  not  have  a  sufficient  time  delay  to  permit  simulation  of 
the  record  because  the  low-takeoff -angle  rays  do  not  have  a  sufficiently 
great  time  delay.  This  effect  is  shown  more  clearly  in  part  (b)  in  which 
all  three  ionospheres  have  the  same  height  and  horizontal  gradient  as  the 
matching  one,  except  that  their  thickness  (^g)  varies. 

It  can  clearly  be  seen  that  the  thinner  the  layer,  the  "thicker”  the 
backsctit ter .  This  should  prove  to  bo  a  very  handy  rule  of  thumb  for 
futui-e  backscattor  simulation.  There  is  another  effect  which  must  be 
taken  into  account,  and  that  is  the  presence  at  high  frequencies  of  an 
uppoi’-ray  structure  in  the  backscattor  trace  only  when  the  layer  is 
relatively  thick.  This  effect  is  shown  very  clearly  in  the  next  illus¬ 
tration,  Fig.  47.  The  two  arrays  of  horizontal  lines  on  this  figure 
wore  drawn  by  a  plotter  which  is  associated  with  the  IBM  7090  used  for 
calculations.  Each  horizontal  lino  on  this  drawing  shows  the  group  time 
delay  and  frequency  of  a  raypath  which  has  boon  computed  through  the 
ionospheres.  Figure  47a  shows  the  results  for  ionosphere  798,  wliich  is 
thick.  The  second  part  shows  the  exact  same  calculation  when  t  lie  iono¬ 
sphere  is  thin,  that  is,  is  33.75  km  instead  of  86.7.  Notice  how- 

thick  the  backscattor  trace  would  be  at  the  lower  frequencies  for  a  thin 
layer,  while  at  the  high  frequencies  the  upper-ray  structure  gives  a 
rounded  shape  to  the  terminus  of  the  backscattor  from  a  thick  layer. 

The  thin  layer,  on  tlie  other  hand,  ends  in  a  sharply  pointed  structure. 

A  comparison  of  those  records  with  the  experimental  one  shown  on  Fig.  27 
gives  a  clear  indication  that  the  ionosphere  was  behaving  like  a  "thick" 
layer  at  the  time  of  the  20  October  experiment. 

In  an  effort  to  make  the  backscattor  "thicker"  without  creating  a 
sharply  pointed  structure  at  the  high-frequency  end,  it  was  decided  to 
insert  an  E  layer  beneath  the  layers  shown  on  Fig.  47.  This  effort 
was  only  partly  successful  in  accomplishing  the  desired  effect,  but  it 
unexpectedly  provided  an  explanation  for  a  number  of  commonly  encountered 
structures  of  sweep-frequency  backscatter. 
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FIG.  47.  COMPUTE R-DRAi\'N  PLOT  OF  R/\YPATH  TIME  DELAY  VS  FREQUENCY. 


121 


SEL-65-002 


A  study  of  the  experimental  records  using  the  slider  technique  of 
Sec.  VIC  indicated  the  possibility  that  there  might  exist  an  E  layer 
with  a  ci-itical  frequency  of  2.84  Me  at  an  altitude  of  120  km.  A  medium 
thickness  (Chapman)  layer  according  to  these  specifications  was  then 
inserted  underneath  each  of  the  two  ionospheres  shovtTi  on  Fig.  47.  In 
both  cases,  the  E  layer  has  a  horizontal  gradient  of  lesser  intensity 
than  the  gradient  in  the  F  layer;  the  critical  frequency  doubles  in  a 
distance  of  9.160  km  at  all  heights  below  120  km.  When  this  E  layer 
was  inserted  under  the  ionosphere  shown  in  Fig.  47a,  calculations  yielded 
the  results  which  are  shown  on  Fig.  48.  Careful  examination  of  part  (a) 
of  the  figure  reveals  the  leading  edge  of  tlie  E  and  F  echoes  as 
discrete  cui-ves  similar  to  those  which  would  be  expected  from  the  data 
previously  giv'en  on  sliders.  The  intersection  of  these  two  curves  appears 
as  a  slight  concav'ity  on  the  near  edge  of  the  echo  which  has  been  labeled 
the  "E-F  inflection."  The  high-frequency  terminus  of  the  echo  still 
has  the  rounded  appeai-ance  characteristic  of  the  F  layer  alone,  as 
-hown  on  Fig.  47a.  This  structure  has  been  termed  the  "upper-ray  flare." 
An  unexpected  i-esult  was  the  appearance  of  the  feature  called  a  "low- 
angle  flare."  Examination  of  the  numerical  data  which  accompanies  these 
calculations  shows  that  the  latter  feature  is  composed  primarilv  of 
energy  which  propagates  at  extremely  low  takeoff  angles.  The  peak  of 
this  flare  occurs  at  exactly  the  frequency  of  the  terminus  of  the 
E-layer  backscatter.  The  E-F  inflection  point  is  at  a  slightly  lower 
f  requency . 

Figure  48b  is  an  experimental  record  acquired  by  the  same  calibrated 
ladar  using  the  same  antenna  at  approximately  the  same  time  of  day  as 
the  i-eeord  which  is  being  synthesized.  The  only  significant  difference 
is  that  this  lecord  was  acquired  two  days  after  the  one  shown  on  Fig.  27. 
The  E-F  inflection  is  visible  and  the  upper-ray  flare  is  clear.  The 
low-angle  flare  is  faintly  visible  on  the  original  record  but  was  too 
weak  to  show  in  the  reproduced  version  of  Fig.  48.  It  was  doubtlessly 
verv  weak  because  of  the  combined  action  of  the  low  antenna  gain  and 
heavy  absorption  at  low  takeoff  angles. 
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GROUP  TIME  DELAY  ( msec) 


UPPER- RAY 


10  16  23 


FREQUENCY  (Me) 

□.  Experimental  record;  0855,  22  October  1964;  same  radar  and 
antenna  as  produced  Fig.  27 

FIG.  48.  EFFECTS  OF  AN  E  LAYER  AND  A  "THICK”  F  LAYER. 


123 


SEL-65-002 


The  broad,  nearly  horizontal  swath  across  the  lower  side  of  Fig.  “ISb 
is  caused  by  sporadic  E.  On  Fig.  49  is  shown  the  result  of  raytracing 
through  an  layer  and  it  can  clearly  be  seen  that  the  presence  of 

such  a  low  thin  layer,  which  has  a  generally  "patchy"  structure,  would 
account  for  the  simultaneous  presence  of  both  tlie  echoes  seen  on  Fig.  48b. 

When  the  same  normal  E  layer  as  described  above  wa.s  inserted  under 
the  thin  ionosphere  shown  in  Fig.  47b,  rayset  calculations  produced  the 
result  sliowii  on  Fig.  .'iOa.  Here  it  can  bo  soon  that  t  lie  sliarply  pointed 
higli-f roquoncy  terminus  of  tlio  thin  F  layer  is  preserved,  while  at  Llie 
same  Lime  the  E-F  inflection  and  the  low-angle  flare  liavo  appeared  on 
the  record.  Tliis  combination  of  features  produces  a  sweep-frequency 
backscattor  record  wliich  seems  to  have  an  abrupt  cutoff  at  a  fixed  ground 
range.  Such  backscattor  records  are  very  common  and  a  good  example  was 
acquired  5  days  prior  to  the  20tli  of  Octolter,  again  using  the  same  radar 
and  the  same  antenna  at  approximately  the  same  time  of  day.  This  record 
is  shown  in  Fig.  50b  where  all  of  the  features  can  be  seen.  .\gain. 
liowever,  tlie  low-angle  flare  is  almost  invisible,  probably  because  of 
tlio  combination  of  low  antenna  gain  and  high  absorption.  To  the  author's 
knowledge,  tlieso  commonly  encountered  structures  have  never  been  explained 
before . 

E.  THE  FINAL  STEP:  MATCHING  THE  AMPl,ITn)E  UECORDS 

Tlie  object  of  this  report  lias  been  the  development  of  a  method  for 
the  study  of  the  amplitude  of  backscattor  signals.  Consetiuen  tly ,  it  is 
of  interest  to  synthesize  the  A  recoiels  of  Fig.  28.  The  first  step 
in  such  a  matching  process  is  the  conversion  of  these  data  to  the  stan- 
diirdized  format  of  Fig.  51.  In  part  (a)  of  this  figure  are  shown  the 
estimated  time  histories  of  the  average  value  of  the  backscattored  signal 
derived  from  the  data  of  Fig.  28  after  all  corrections  have  been  made  for 
measurements  on  the  calibrated  radar  parameters.  Great  care  has  boon 
exorcised  in  the  acciuisition  of  those  c|uanti  tativo  data;  a  number  of  the 
calibration  instruments  wore  oven  returned  to  their  original  manufacturers 
to  permit  final  small  corrections  to  the  numbers  given  on  Fig.  51.  The 
peak  amplitudes  are  within  3  db  of  similar  measurements  published  by 
Shearman  [hof.  20;. 
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b.  Experimental  record:  0929,  15  October  1964;  same  radar  and 
antenna  as  produced  Fig.  27.  (Note  the  apparent  shape 
change  due  to  ordinate  compression.) 

FIG.  50.  EFFECTS  OF  AN  E  LAYER  AND  A  "THIN"  F  LAYER. 
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from  data  on  P  =  1  kw 

Fig.  28 

FIG.  51.  EXPERIMENTAL  "A"  RECORDS  IN  A  STANDARDIZED  FORMAT. 
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In  part  (b)  of  Fig.  51,  the  same  data  are  presented  but  the  vertical 
scale  is  slightly  adjusted  so  that  the  effective  transmitted  power  is 
1  kw  at  each  of  the  six  calibrated  frequencies.  A  seventh  record  was 
acquired,  at  21.5  Me,  but  it  was  so  far  removed  from  the  other  records 
in  time  that  it  was  not  possible  to  calibrate  the  group  time  delay. 
Consequently,  it  was  not  included  in  this  analysis. 

The  rest  of  this  chapter  is  devoted  to  the  synthesis  of  backseat  ter 
records  which  match  the  amplitudes  given  on  Fig.  51b.  As  a  basis  of 
calculation,  the  ionospheres  of  the  preceding  section  will  bo  used  oven 
though  they  wore  synthesized  to  match  the  record  of  1030  PDT.  Conse¬ 
quently,  there  will  be  small  errors  in  the  group  time  delays  of  the 
records  which  follow  but  these  could  be  removed  by  small  changes  in  the 
ionosphere  and  subsequent  re-ray trac ing.  However,  it  will  be  soon  that 
the  timing  error  inserted  through  the  use  of  the  same  ionosphere  for  all 
A  records  is  acceptable  and  the  added  expenses  of  raytracing  through 
six  new  and  different  ionospheres  would  doubtless  cause  an  insignificant 

cliange  in  the  conclusions  which  follow. 

Figure  52  illustrates  a  number  of  things:  in  the  lower  left  cornei 
of  part  (a),  note  the  synthetic  received  pulse  which  is  drawn  by  the 
computer  in  order  to  show  the  calculated  shape  of  the  echo  from  a  single 
scatlerer.  It  is  a  direct  plot  of  the  sequence  of  numbers  which  were 
cho.sen  to  describe  the  received  pulse.  (The  absolute  values  plotted  are 
not  significant.)  This  shape  is  identical  to  the  bottom  picture  on  Fig. 
33  except  for  the  obvious  compression  of  the  horizontal  scale  and  smooth¬ 
ing  over'  0.1 -msec  Lime  intervals. 

Also  shown  on  Fig.  52a  is  a  dotted  outline  showing  what  might  bo 
considered  as  the  "useful  portion”  of  the  plot.  Since  this  figure  is 
supposed  to  show  a  synthetic  version  of  the  record  which  might  be 
acquired  by  a  real  backscatter  radar,  there  should  be  some  noise  level 
below  which  echoes  cannot  be  received.  In  the  computer  there  is  no 
noise  level  but,  because  of  ihe  discrete  integration  stops  used,  the 
quality  of  this  synthetic  record  deteriorates  at  the  lower  signal  levels. 
In  fact,  it  was  through  trial  and  error  with  different  integration  steps 
that  the  author  determined  which  values  yield  good  results  on  signals 
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strong  enough  to  be  received  experimentally.  It  would  bo  wasteful  to 
use  smaller  integration  steps  since  this  would  increase  the  cost  without 
increasing  the  useful  quality  of  the  product.  For  the  analysis  which 
follows,  a  large  number  of  backseat  ter  records  have  boon  run  and  in  each 
case,  only  the  "useful  portion"  v^ill  bo  presented. 

The  backseat  tor  record  of  Fig.  52a  has  been  selected  because  it  is 
ail  exceptionally  good  illustration  of  the  strength  of  "mixed  mode"  echoes. 
Those  are  the  echoes  composed  of  energy  which  travels  out  to  the  ground 
and  then  returns  by  a  different  mode.  For  example,  the  energy  may  go 
out  by  two  hops  and  return  by  one  hop.  Such  energy  accounts  for  the  peak 
in  the  received  signal  between  the  one-  and  two-hop  peaks  on  Fig.  52a. 
Similarly,  energy  which  goes  out  by  two  hops  and  returns  bv  three  (or 
vice  versa)  accounts  for  the  peak  which  is  located  approximately  midway 
between  the  two-  and  three-hop  echoes  on  Fig.  52a. 

The  fact  that  these  peaks  are  indeed  due  to  energy  which  is  propagated 
by  mixed  modes  is  graphically  shoini  in  part  (b)  of  the  figure  where  the 
same  calculation  is  carried  out  by  the  synthesis  program  after  it  has 
been  modified  so  that  it  cannot  calculate  mixed  modes.  It  can  thus  be 
seen  that  in  some  circumstances  the  mixed-mode  echoes  may  account  for  a 
considerable  amount  of  the  receivetl  energy.  At  the  time  of  this  writing, 
the  author  had  not  yet  demonstrated  the  existence  of  such  echoes  on 
experimental  records,  but  it  is  hoped  that  this  will  be  possible  through 
use  of  the  calibrated  radar. 

The  first  attempt  to  match  experimental  A  records  resulted  in  Fig. 

52  which  shows  a  dozen  calculations  using  t  lie  antenna  patterns  derived 
from  Fig.  32.  The  absorption  values  are  approximately  ecpial  to  the  0900 
data  given  in  Table  2.  Part  (a)  of  the  figure  shows  the  situation  when 
the  backscatter  coefficient,  K,  is  a  multiplier  of  sin  ;.  Part  (b) 
of  the  figure  shows  the  identical  situation  when  the  backscatter  coef¬ 
ficient  is  independent  of  '.  Notice  that  both  sets  of  data  produce 
spikes  which  are  too  thin  (not  sufficiently  spread  in  group  time  delay) 
but  which  liavc  amplitudes  approximately  cori-ect.  However,  in  Fig.  53a 

4 

the  multiplier  of  sin  ,  had  to  be  increased  by  nearly  an  order  of 
magnitude  for  each  2-Mc  increase  in  frequency.  At  20  Me,  it  reached  a 
value  of  100  which  is  unrealistic  because  values  much  in  excess  of  1 
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violate  the  principle  of  conservation  of  energy  (e.g.,  see  the  definition 

of  K(i|r)  in  Sec.  IIB).  To  receive  echoes  that  large  would  require  a 

source  on  the  ground.  Not  all  synthetic  backscatter  is  shown  here,  but 

a  number  of  other  calculations  similar  to  Fig.  53a  produced  the  same 

indication  that  the  predominant  source  of  backscatter  could  not  be  pro- 
4 

portional  to  sin 

Figure  53b  is  an  excellent  match  for  the  experimental  data  except 
for  the  narrowness  of  the  spikes,  as  has  been  mentioned.  Notice  that 
the  backscatter  coefficient  K  is  a  function  of  neither  frequency  nor 
angle  of  incidence.  It  is  merely  a  constant.  Both  portions  of  Fig.  53 
were  calculated  using  the  "thick”  ionosphere  shown  on  Fig.  47a. 

Figure  54  shows  two  more  sots  of  calculations,  each  of  which  should 

be  compared  to  Fig.  53b.  Figure  54a  shows  the  situation  when  the  back- 

2 

scatter  coefficient  varies  as  sin  Here,  the  situation  is  similar 

4 

to  that  encountered  when  sin  was  used;  that  is,  the  multiplier  of 
the  sine  function  increases  with  frequency  and  gets  quite  large.  In 
fact,  at  20  Me,  the  multiplier  is  almost  largo  enough  to  violate  the 
energy-conservation  principle.  It  can  bo  concluded  that  a  sin^ijr 
dependence  of  the  backscattoring  coefficient,  although  not  impossible, 
appears  lo  bo  much  less  likely  than  a  simple  constant. 

The  major  objection  to  the  synthetic  backscatter  calculated  thus  far 
has  boon  its  insufficient  spread  in  time  delay.  An  examination  of  Fig.  47 
indicates  tliat  bettor  results  would  be  obtained  with  the  "thin"  ionosphere 
of  part  (b).  The  result  of  this  substitution  is  shown  in  Fig.  54b.  iiero, 
and  in  all  the  calculations  which  follow,  the  10  AM  absorption  values  are 
used . 

It  can  bo  seen  that  Fig.  54b  is  considerably  better  as  a  match  of 
the  experimental  data  than  were  any  of  the  previous  calculations.  Novor- 
theless,  the  returns  are  still  too  thin  and  a  study  of  the  numerical 
data  \thicli  accompany  these  plots  indicated  that  the  energy  along  the 
trailing  edge  of  the  echoes  was  being  severely  attenuated  by  the  very 
low  antenna  gain  at  low  angles.  The  experimental  antenna  measurements 
were  again  examined  and  it  appeared  advisable  to  increase  the  assumed 
gain  of  the  antenna  at  low  angles  by  a  small  amount.  Figure  55  shows  the 
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Thick  Chapman  layer  b.  Thin  Chapman  layer 

FIG.  54.  SECOND  ATTEMPT  AT  SYNTHESIS. 
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I'lG.  r>5.  SYNTHETIC  ANTENNA -PATTEKN  CHANGE  BASED  ON 
A I KCUAFT  MEASUREMENTS . 


synthetic  gain  pattern  at  14  Me  before  and  after  this  change  was  accom¬ 
plished.  At  other  frequencies,  the  situation  was  similar.  Barnes  [uef. 
45]  made  no  measurement  of  the  antenna  at  the  frequencies  used  for  this 
synthesis,  t)ut  his  measurements  at  nearby  frequencies  indicated  that  the 
gain  at  low  angles  was  in  many  cases  considerably  higher  than  indicated 
by  the  theoretical  curves  which  have  been  used  until  now. 
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Figure  56  shows  the  new  calculations  with  the  thin  ionosphere  (l78) 

and  the  modified  antenna  patterns.  Two  cases  are  shown:  a  constant 

2 

backscatter  coefficient  and  a  sin  \lf  backscatter  coefficient.  A  cal- 

4 

culation  was  also  carried  out  for  sin  'Ja  and  again  the  result  was  found 
to  be  unrealistic. 

It  can  be  seen  that  Fig.  56a  would  be  a  good  match  for  the  experimental 
record  except  for  the  presence  of  a  very  strong  loading  edge  on  the  echoes. 
Examination  of  the  printed  data  showed  that  all  of  this  energy  came  from 
rays  which  had  apogees  in  excess  of  253  km.  If  the  actual  ionospheric 
electron-density  profile  on  the  20th  of  October  followed  the  general  con¬ 
tour  of  the  assumed  ionosphere  only  up  to  253  km  and  above  that  altitude 
the  electron  density  decieased,  then  there  would  bo  no  rays  with  apogees 
above  253  km.  This  is  a  reasonable  assumption,  since  rocket  measurements 
of  the  electron-density  profiles  usually  show  data  curves  which  are  not 
smooth  like  the  Chapman  layers  used  for  synthesis  here. 

Using  this  assumption,  final  calculations  wore  carried  out  and  the 
results  are  shown  on  Fig.  57b.  It  can  bo  seen  that  the  match  between 
the  synthetic  and  experimental  records  is  quite  good  at  the  four  lower 
frequencies.  At  18  and  20  Me,  the  echoes  are  still  too  thin.  Examination 
of  Fig.  36  stiows  that  tliese  two  records  were  acquired  more  than  an  hour 
prior  to  the  time  of  the  Z  record  which  was  matched  by  use  of  the  thin 
ionosphere.  Consequently,  it  is  not  surprising  tliat  the  agreement  is 
poor  at  the  higli  frequencies  and  good  at  the  low  frequencies  because  the 
latter  were  measured  last  at  a  time  not  far  removed  from  the  Z  record. 
Figure  57a  sliows  the  experimen  .il  data  of  Fig.  511)  so  tliat  it  may  be 
directly  compared  to  tlie  final  syntlietic  data. 

From  tiie  amplitude-matching  exercise,  tlie  following  conclusions  and 
comments  can  be  made: 

1.  A  ground  reflection  coefficient  which  is  neith('r  frequency  nor 

angle  sensitive  can  account  for  the  experimental  backscatter  record. 
In  terms  of  the  more  familiar  radar  parametc-j  cr^  this  implies  a 
variation  like  sin  y.  However,  other  variations  are  not  ruled  out. 
It  can  only  bo  said  at  this  time--in  the  light  of  the  match  between 
experimental  and  computed  resul ts--tha^  if  the  coefficient  is  a 
function  of  angle,  then  it  must  be  a  function  of  frequency;  if  it 
is  not  a  function  of  angle  it  is  not  a  function  of  frequency  within 
tile  10-20  Me  band.  Future  work  with  different  ionospheres  and 
better  records  should  permit  a  decisive  resolution  of  this  matter. 
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FIG.  56.  FURTIER  SYNTHESIS. 
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GROUP  TIME  DELAY  (msec) 

(o)  (b) 

FIG.  57.  COMPARISON  OF  EXPERIMENTAL  AND  SV'NTHETIC  RESULTS. 
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2. 


Comparison  between  synthetic  and  experimental  backscatter  has 
been  severely  hindered  by  the  long  time  interval  for  acquisition 
of  records.  In  future  work,  emphasis  will  be  placed  on  the  mea¬ 
surement  of  A  and  Z  records  in  a  short  interval  of  time  so 
that  they  can  be  used  together  as  a  set  without  inaccuracy  due 
to  changes  in  the  ionosphere. 

.3.  The  antenna  pattern  has  a  major  influence  on  the  time  distribution 
of  received  backscatter  energy.  The  strength  of  this  influence  can 
be  seen  clearly  by  the  program  user  because  printed  data  which 
accompany  the  record  break  the  calculation  down  into  stages  where 
each  individual  factor  can  bo  seen.  On  occasion,  when  the  wrong 
synthetic  antenna  was  accidentally  used,  the  backscatter  echo  was 
severely  distorted. 

4.  Future  work  of  this  nature  will  be  greatly  aided  by  a  more 

thorough  study  of  the  effect  of  various  ionospheric  structures 
on  the  general  shape  of  backscatter.  A  number  of  features  are 
still  not  understood  and  the  metnod  given  hero  seems  to  show 
promise  of  being  able  to  answer  many  questions.  The  synthesis 
of  an  ionosphere  to  match  the  20  October  record  would  have  boon 
much  easier  if  it  could  have  followed  a  thorougli  study  of  the 
effect  of  ionospheric  layer  structures  on  swoop-frequency  back¬ 
scatter  shapes. 
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13.  ABSTRACT  ^  hlgh-f roqucncy  pul.so  radar  often  receives  energy  which  has  propugaied 

through  tlK.'  ionosphere  to  the  ground,  scattered,  and  then  returned  to  the  radar 

again  via  the  ionosphere.  When  this  energy  is  displayed  as  a  function  of  botli 

tinu.'  del-jy  and  operutln,g  frequency,  tlic  record  is  called  sweep- frequency  back- 

seta'.,  ter.  Tn  this  report  a  digital  computer  technique  is  devised  to  synthesize  such 

records  in  a  manner  which  is  quantitative  in  all  parameters. 

.'’he  synthesis  is  designed  .so  tliut  very  few  approximations  are  needed.  Full 

allowance  is  made  for  the  following;  u  spherical  earth  and  iono.sphore;  a  wide 

variet.y  t)  f  electron-density  variations  in  height  and  range;  Ionospheric  absorption 

which  v.Ti'ie.s  wi  tli  frequency,  azimuth,  and  elevation;  ground-scattering  cliarac  ter  1  s  i 

wliicli  vary  wi  '  li .  f  req..iency ,  range,  and  the  angle  of  incidence;  ground-reflection 

coe  f  f  i  0  I  onts  which  vary  with  frequency  and  range;  the  tran.smitted  pulse  shape, 

'.ind  power  level;  receiver  bandwidth  and  total  received  noise  power. 

Only  three  approximations  are  made:  the  discrete  number  systems  is  subs 1 1  tu ti.'d 

for  contlnuou.s  numbers;  ray  theory  is  used  with  the  assumption  that  energy  1:. 

distributed  uniformly  between  adjacent  closely  spaced  raypaths;  tiie  effect  of  Llie 

geomagnetic  field  on  rayputh  range  and  group  time  delay  is  neglected. 

The  digital  synthosi.s  has  two  major  uses.  First,  it  is  pcjssiblc  to  vary  a 

single  parameter  in  a  known  way  and  then  sec  the  effect  that  the  variation  should 

have  on  ground  backscatter .  For  example,  a  ripple  of  known  characteristics  can 

bo  inserted  in  the  ionosphere  at  a  known  range  and  the  resulting  perturbation  in 

the  ground  backscatter  can  be  calculated.  / r,  .  <  ,  n\ 
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The  second  major  use  of  the  program  is  the  synthesis  of  a  close  likeness 
of  a  specific  experimental  record.  This  .synthetic  record  serves  to  confirm 
knowledge  of  the  mechanisms  which  play  a  part  in  ground  backscatter.  For 
example,  the  ground-scatter  coefficient  varies  with  frequency,  terrain  typo, 
and  incidence  angle  in  a  manner  which  is  not  presently  well  undor.siood.  The 
process  of  making  a  synthetic  record  to  match  an  experimental  record  will 
involve  the  trial  of  several  hypothetical  scatter  functions  until  one  is 
found  which  duplicates  observed  behavior.  Since  most  of  the  other  con¬ 
tributing  mechanisms  exhibit  a  predictable  behavior,  it  will  thus  be  possible 
to  isolate  the  scatter  pattern  and  study  it  independently. 
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